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SUMMARY 

The problem of flow of a compressible fluid past a "body 
with subsonic flow at. infinity is formulated by the hodograph 
method. The solution in the hodograph plane is first con- 
structed about the origin by superposition of the particular 
integrals of the transformed equations of motion with a set 
of constants which would determine, in the limiting case, a 
known incompressible flow. This solution is then extended 
outside the circle of convergence by analytic continuation. 

The previous difficulty ?f the Chaplygin method of slow 
convergence of the series has been overcome by using the as- 
ymptotic properties of the hypergeometric functions so that 
numerical solutions can be obtained without difficulty. It 
is emphasized that, for a solution covering the whole domain 
of the field of flow, both fundamental solutions of the hyper- 
geometrical differential equation are required. 

Explicit formulas for numerical calculations are given 
for the flow about a body, -such as an elliptic cylinder, and 
for the periodic flow such as would exist over a wavy surface. 

Numerical examples based on the incompressible flow so- 
lution of an elliptic cylinder'of thickness ratio of 0.6 are 
computed for free-stream Mach numbers of 0.6 and 0.7« 

The results of this investigation indicate an appreciable 
distortion in the shape of the bodies in compressible flow 
from that of incompressible flow, which necessitates a series 
of computations with various values of the geometric parameter 
in order that the desired body shapes can be selected for a 
given Mach number. It also is shown that the breakdown of ir- 
rotational flow depends solely upon the occurrence of limiting 
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lines, which, in turn, are dependent on the boundary condi- 
tions. 

The numerical calculations show that at a free-stream 
Mach number of 0.6, irrotational supersonic flow exists up to 
a lfecal Mach number of 1.25; whereas breakdown occurs at 1.22 
for a Mach number of 0.7. 


INTRODUCTION 


When a flow of nonviscous incompressible fluid is irrota- 
tional, it is well known that the problem can be reduced to 
either the problem of Dirichlet or that of Neumann, and that 
there exists a unique solution for any given boundary condi- 
tions. When the fluid is nonviscous but compressible, the va- 
riation of density makes the mathematical problem very diffi- 
cult and complex. In this case, a pure pdtential flow through- 
out the region is not alv/ays possible for a given body; this 
depends very much upon the condition at infinity. If a certain 
speed of the flow at infinity is reached, regions within the 
field of flow will be created in which the irrotational flow 
does not exist owing to the appearance of "limiting lines." 

Such regions were picturesquely designated a6 "forbidden re- 
gions" by Th. von Kdrmdn (reference l), and they appear when 
the local speed of the flow considerably exceeds the local 
speed of sound. It has been shown that the occurrence of lim- 
iting lines is directly connected with the breakdown of irro- 
tational flow and with the resultant increase in drag of the 
body due to shock waves. In other words, if there is a limit- 
ing line in the field of flow, the isentropic irrotational 
flow must break down. However, the irrotational flow may break 
down before the appearance of limiting line due to the insta- 
bility. of the velocity field. On the other hand, shock waves 
can occur only in supersonic flovi. Therefore, there is no 
danger. of breakdown of isentropic flow if the whole field of 
flow i^ subsonic. C on sequent ly , the Much number correspond- 
ing to^the first appearance of local speed equal to that of 
sound can be designated as the "lower critical Mach number"; 
and thd Mach number corresponding to the first appearance of 
limiting lines can be designated as the "upper critical Mach 
number J 1 The actual critical Mach number for a given body 
will be! influenced by the boundary layer and hence the Reynolds 
number., However, it must lie between these two limiting crit- 
ical values. (See reference 2.) Thus, knowledge of these 
critical speedy of the flow are essential for the design of 
efficient aerodynamic bodies. 
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To determine the critical Mach numbers, the general prob- 
lem of flow of a compressible fluid about a given body must 
be solved. The often-used methods treating such a problem 
are Janzen-Rayleigh' s method of successive approximations and 
Glauert-Prandtl 1 s method of small perturbation. The latter 
method has been extended recently by both Hantzsche and Vfendt 
(reference 3) and C. Kaplan (reference 4). Indeed, both meth- 
ods yield valuable information regarding the effects of com- 
pressibility and are useful for many practical design problems, 
particularly the determination of the lower critical Mach num- 
ber of a given body. But, so far as the general problem of 
limiting line and upper critical number is concerned, none 
seems to be adequate, owing to the doubtful convergence of 
such successive approximations at the required high Mach num- 
bers . 

An entirely different approach first was made by Molenbroek 
(reference 5) and Chaplygin (reference 6) by introducing the 
velocity components instead of the usual space coordinates as 
independent variables. The advantage of the method is that, 
instead of a nonlinear differential equation as is the case in 
the physical plane, it leads to a linear one in the velocity or 
hodograph plane. The particular solutions of this linear equa- 
tion are found to be products of trigonometric functions of the 
angle of inclination of velocity vector and hypergeometric 
functions of the magnitude of the velocity vector. It is then 
possible to construct a general solution from the particular 
solutions of the differential equation. The difficulty, how- 
ever, is that the character of the field in the physical plane 
to which the solution in the hodograph plane corresponds cannot 
be determined beforehand. This difficulty prevents the exact 
formulation of the boundary value problem in the hodograph 
plane. Chaplygin has overcome this handicap by first choosing 
a “suitable solution” in the hodograph plane and then proceeding 
to find the corresponding flow in the physical plane. The 
euitable solution is one which, in the limiting case of zero 
Mach number at infinity, becomes identical with the incompress- 
ible flow over a body similar to the body concerned. This will 
ensure the satisfaction of the proper boundary conditions in 
the physical plane. Furthermore, such a solution would be ex- 
act both for the subsonic and for the supersonic regions, as 
no approximation is introduced. Therefore, it is particularly 
suitable for the problem of determining the upper critical Mach 
number for a given body, as limiting lines occur only in mixed 
subsonic and supersonic flows. This method is followed in the 
present report, except for the introduction of the transformed 
potential function for easy calculation' of the space 

coordinat e s . 



4 


NACA TN No. 955 


For the flow around a hody, Chaplygin’s procedure will 
lead to a solution in the form of an infinite series, each 
term of which is a product of a trigonometric function and a 
hypergeometric function. To put the method on a firm founda- 
tion, it is necessary to establish the convergence of the in- 
finite series. Chaplygin himself has done this for the sub- ■ 
sonic region. Thus, only the extension to include the super- 
sonic region remains to be completed. In part I of this re- 
port, the general properties of hypergeometric functions of 
large order are investigated in preparation for the proof of 
the convergence given in part II. The essential point in 
these parts is to establish the upper and lower bounds for 
the hypergeometric functions so that the sum of the infinite 
series can be discussed. It is appropriate to mention here 
that for the proper representation of the general solution 
in the holograph plane, both fundamental solutions of the hy- 
pergeometric differential equation are required. This fact 
has not been considered by many of the previous investigators 
in this field. In other cases (reference 7) the investigator 
has chosen to work with only the first solution. 

The general solution constructed by the Chaplygin method 
is really an existence theorem. The extremely slow, conver- 
gence of the series makes numerical calculation very difficult, 
if not impossible. This, in fact, constitutes the main diffi- 
culty of the method. In part III of the present report, this 
difficulty is svercome by using the asymptotic properties of 
the hypergeometric functions. The result is the separation of 
the solution in the holograph plane into tvo parts. One part 
is of clos.ed form and is the product of a universal function 
of the velocity and the same solution as for incompressible 
flow but with a velocity distortion, or velocity correction. 

For instance, the first part -of the stream function for the 
compressible flow is equal to the product of the universal 
function of velocity and the stream function for the incom- 
pressible flow with the magnitude of velocity modified by a 
given rule. The other part is an infinite series which con- 
verges rapidly everywhere except in a small region on both 
sides of a critical circle with a radius equal to q = c in 
the holograph plane. In practice, by using only a few terms 
of the infinite series, this zone of slow convergence ca&'be 
limited to such a small interval that, it is of no consequence. 
Thus the Chaplygin procedure is improved to a point where ac- 
tual numerical calculations can be made without difficulty. 

As a result of this part of the study it becomes cl ear 
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that h y the mere substitution of a different speed scale, or 
velocity distortion, in the solution for an incompressible 
fluid, an accurate enough solution for the compressible flow 
cannot be obtained. For if this were the case, then not only 
the second part of the solution (the rapidly convergent se- 
ries given by the present method) would be negligible, but 
also the value of the multiplying universal function of veloc- 
ity in the first part of the solution would be unity. How- 
ever, the value of the second part of the solution is not 
small compared with that of the first part for' a speed near 
that of sound, and the value of the multiplying function of 
velocity i s f ar from unity. In other words, the usual so- 
called hodograph method (reference 8) cannot, in general, 
yield satisfactory results, for mixed subsonic and supersonic 
flow. On the other hand, the present method doeB show that 
the second part of the solution is zero and the multiplying 
function in the first part takes the constant value of unity, 
if the isentropic exponent is equal to -1 . This means that 
for this particular case, a simple speed distortion is suffi- 
cient. This is, of course, in accordance with the previous 
investigation of von Karman (reference l) and Tsien (reference 
9) and L. Bers (reference 10). 

Furthermore, the present method also shows that the rules 
of speed distortion for the first part of the solution can be 
used . only for subsonic flow and that there is a singularity at 
the local sonic speed. For regions of supersonic flow, the 
first part of the solution involves both the incompressible 
stream function and the incompressible potential function. 

Thus even without considering the second part of the solution, 
there is no possibility of making the compressible stream 
lines coincide with those for incompressible flow in the hodo- 
graph plane by a simple stretching of the speed scale. The 
mathematical basis of this fact is the change in character of 
the differential equation from elliptic to hyperbolic in the 
transition from subsonic to supersonic flow. For the super- 
sonic regions, it is not possible to use a real transformation 
of the velocity variable to convert the differential equation 
of flow to the Laplace equation, and thus make a simple con- 
nection between the compressible and the incompressible flows. 
This is one of the difficulties of the previously proposed 
hodograph method. In fact, writers using this method must 
generally limit their calculation to subsonic speeds. (See 
references 9, 10.) Now this limit is removed, and the whole 
field of mixed subsonic and supersonic flows can be treated at 
once with ease . 


For the purely subsonic flow, the second part of the 
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solution is small compared with the first part and may he 
neglected. furthermore, if only the zero streamline repre- 
senting the body is considered, the universal multiplying 
function of velocity is of no importance. In other words, 
for this case, a simple speed distortion from the solution of 
incompressible flew is -sufficient to give accurate enough re- 
sults. However, the subject of the ,, best tt velocity distor- 
tion rule in subsonic regions ha3 been the subject of many 
discussions. (See references 1 and g.) The present analysis 
id considered t« settle this question. This is due to the 
fact that the present velocity distortion rule is obtained 
from the asymptotic properties of the hypergeometric functions, 
and that such properties are definite and unique. Therefore, 
the resultant velocity distortion rule is not the result of 
uncertain speculation. furthermore, it is also the best rule, 
because the analysis implies that this rule will make the 
second part of the solution, or the correction terms, the 
smallest. This distortion rule is found to coincide with that 
of Temple and Yarwood. (See reference 11.) 

for the purely supersonic flow, the second par't"sf-~the - 
solution is again small compared* vrith the first part and may 
be neglected. In fact, the solution then can be reduced to 
that of the simple wave equation with the inclination of the 
velocity vector and the distorted velocity as independent va- 
riables. This is, of course, the counterpart of the fact that 
by a simple distortion in velocity, the, differential equation 
for subsonic flows can be reduced to the Laplace equation. 

The usefulness of this new result fo-r purely supersonic flow 
has yet t* be exploited. 

Once the general problem of mixed subsonic and supersonic 
flow around a body is solved, the determination of the upper 
critical Mach number or the Mach number for the first appear- 
ance of the limiting lines is a simple matter. This problem 
is discussed in part IT of the report. A simple method is de- 
veloped, based on the properties of the limiting line as given 
by von Karman (reference 1), Ringleb preference 12), Tollmien 
(reference 13 ) , and Tsien (reference 2 ) ' 

To test the practicability of. the method developed, tw# 
numerical examples are worked out in detail. However, in 
order tq reduce the amount of computational work and in view 
of the limit ed 1 t ime available, a slightly different procedure 
actually is used. This procedure is only approximate but is 
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"believed to "be sufficiently accurate in the supersonic region 
to give a satisfactory description of the most interesting 
features of such flows. The examples chosen are derived from 
the incompressible solution of an elliptic cylinder of thick- 
ness ratio 0.6. The free-stream Mach numbers of the compress- 
ible flow are 0.6 and 0.7 for these two examples. The first 
case gives a smooth flow over an "elliptic 11 cylinder of thick- 
ness ratio 0.42. The maximum local Mach number is approxi- 
mately 1.25. Thus a considerable supersonic region exists. 

The second case gives a flow with limiting line. 

Finally, it must be said that owing to the limitation of 
time, only the case of flow without circulation is investi- 
gated in detail. The explicit formulas for numerical calcu- 
lation are given for two cases: (a) Flow around a body such 

as an ellipse, (b) peri-odi-e-fl-ow -pa^t4^r-n_.aucli__as that over a 
wavy surface. However, it is believed that nor e~‘gen-ers.l~ case a 
can be studied by a slight extension of the present results 
and use of the same method of approach. 

This investigation, conducted at the G-uggenheim 
Aeronautics Laboratory, California Institute of Technology, 
was sponsored by and conducted with the financial assistance 
of the National Advisory Committee for Aeronautics. 


hot at ions 

The symbols used in this report are classified according 
to the following groups: 

A. Physical Quantities 

x,y Cartesian coordinates 
u,v the velocity components 

q the absolute value of the velocity vector 
d the Inclination of the velocity vector with x-axis 
p density of the fluid 
P 0 density of the fluid at q = 0 
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p pressure within the fluid corresponding to p 
p 0 pressure at q =» 0 
y ratio of the specific heats 
c the local speed of sound 

c Q the speed of sound at q e 0 

U the value of q at infinity, assuming parallel to the 

x-axis. With subscript, however, it may he a function 
of T. 

B. Hydrodynamic Junctions in the Physical Plane 
z » x + iy 

W 0 (z) = c P 0 (x,y) + i\j/ 0 (x,y) complex potential for incompresB- 

ible flow in z 

✓ 

cp 0 velocity potential for incompre ssihle flow 
\|/ 0 stream function for incompr essihle flow 
cp velocity potential for compressible flow 
v|/ stream function for compressible flow 

0. Hydrodynamic Junctions in the Hodograph Plane 
w ss u - iv 

W 0 (w) ss cp Q (u,v) + i \{f 9 (u , v ) complex potential for incompress- 
ible flow in w 

cp 0 ("u,v) - velocity potential for incompressible flow 

'l' 0 (u,v) stream function for incompressible flow 

A 0 (w) - zw - W Q (w) = X 0 (u,v) - iC7 0 (u,v) transformed complex 

•ay potential function 

X 0 (u.v) B ux + vy - cp 0 (x,y); x = 

a u 

:w 

y s 0 transformed potential 
8v function 
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W(-w; T ) 
\J/(u,v) 
A(w; T ) 

X(u,v) 


the complex potential function for compressible flow 


= Im 


|w(w;T)| 


Stream function for compressible flow 


transformed complex potential function for compress- 
ible flow 


= ux + vy - cp(x,y) 


HI 


{■ 


A( w; T ) 


} 


transformed potential 
function for com- 
pressible flow 


® 0 (u,v) 


dXp 

3£ 




S£o 
3 •§ 


\J/( q , -d ) = ^(q.S) + \[/ 3 ( 1 ) ( q, £) ; ^(q.d) represents the con- 

tribution by the 
velocity distorsion; 

\|/ s ( * ^ (q, $) stands 
for the transformed 
infinite series, 
where the super- 
script 1 may ei- 
ther mean i the 
inner, or o the 
outer solution. In 
the case of coordi- 
nates, the notation 
is exactly the same. 

= J* (t)AB Jo) + E n. A lp( T J 
U ^ fCr^T^T,) 


Gp (a) (T) 


£ v ( T)AB n (a) 


B n A| i; (T) 

+ f (T 1 )T u ( T x ) 




= p - 1 

V + 1 


Ip 


l (T)AB a (») 


+ 


5nAf„.,CO 

f(T J )T V (T 1 ) 
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E. Parameters and Variables 
v positive rational numbers 
m t n positive integers 

a denotes 1 or 2 when used as superscript with a bracket 
or a = 



P 


A 


denotes the dependence on P when used as subscript 

or p * — 1 

7-1 


3(2p ) a ^ 3 1 1 

(1 + a) a T(t ] l ) 


the ratio of the distorted speed 
to that at infinity 



€ geometrical parameter of the body 

A laplacian or difference between exact and approximate 
values of a function or a constant 
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E. Hyperge ometr i c Functions 


a,b,c parameters of the hyperge ometr ic functions. In par- 
ticular, , "by, c v are defined by (29). 

Xy ( T ) = F(a^j, b.y ; C-yjT) fir st ' integral of the hypergeometric 

equation associated with the 
stream function 


x_ v) (' r ) a F.(l + a v - c V t 1 + - Oyj 2 - C v i T) 


v ~ w v l 

-V 


V T > - 


TT T 




.1 


T‘(a v )T (b u ) 


v 


(2pc 0 2 ) T( C y-l)T(cy) LT(i + a v ,Cy)T(l+V°v J 


T (T) ' 


V 


T(c v ) 


1 „(T) 


T(2~c v ) ~~ V 

2 v (r) - q au F v( T ^ 

l Vil (T) = F(l+a ut 1 +Td v ; 1 +c^j t) 

l v (r) ( T) = X 1 J (T)/X u (T l ) 
l v> 1 (r) (r) = F v> 1 (t)/f u (t 1 ) 

X^( t) = F ( t) + i Fy( t) 

B U (T) * F* v (t) 
cj>y(t) « arg I* (T) 


second integral of the same equation 


If any function or a constant is associated xvith xLl»^)» 
it will be marked on top by a symbol such as Xp( T )» 
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PART I 

DIFFERENTIAL EQUATIONS OF OOMPRESSIBLE FLOW AND 
' PROPERTIES OF THEIR PARTICULAR SOLUTIONS 
1. Equations of Motion 


It is proposed to study the irrotational steady motion of 
an inviscld nonconduot ing compressible fluid in an infinitely 
extended domain containing a cylindrical body with itB axis 
perpendicular to the constant . velocity at infinity. The flow 
is then two-dimensional. Let x and y be the Cartesian co- 
ordinates and u and v the velocity components parallel to 
the x- and the y~axis. The dynamical equations governing such 
a motion, in the absence of body force, are 


pu + p T If 


0v , dv 

pu a; + pT sF 


_ op 

(1) 

OI 


_ 5p 

(2) 

oy 



Here p is the pressure and p the density of the fluid, 
both being continuous functions of x and y. In addition, 
the following equation of continuity must be satisfied: 

— (pu) + — (pv) = 0 (3) 

3x ^y 


Furthermore, since the velocity is constant at infinity, the 
flow is irrotational there. Then, according to Thomson's 
theorem, if the pressure is a function of the density alone, 
the flow' will remain irrotational; that is. 


8 v 8u 
3x 3y 


0 


(4) 


In the case of flow of an inviscid nonconducting gas, the 
thermodynamic change of state of the gas is adiabatic. If 
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the flow is assumed to he continuous, excluding shock waves, 
then the relation between p and p must he that of an 
isentropic process: 

p es constant (5) 


where 7 is the ratio of the specific heats. 


As in the case of incompre ssihle flow, there are more 
equations than the number of the variables. However, by vir- 
tue of equations (4) and (5), the dynamical equations (l) and 
(2) reduce to a single differential equation and can be inte- 
grated easily to give a relation between the pressure and the 
magnitude q of the velocity; namely, 

7 

P = P?J -alj V-\ with q2 = u* + v3 (6) 


Here 

speed 

It is 
<1 t>y 


Po 

of 


and c G are respectively the pressure 
sound at the stagnation point q » 0 and 


and the 



possible to obtain a similar relation between p and 
means of equation (5): 


P - P 0 



(7) 


where p 0 denotes the value of p at q = 0. 


After integrating the dynamical equations, the velocities 
u and v can be determined from the kinematic conditions 
specified by equations (3) and ( 4 ). By eliminating p from 
equation (3;, the result is 


(1 - 


2uv 6u + ( 
c 3 3y ^ 


v 3 \ 8v 
c 3 ' 3y 


0 


( 8 ) 


where c 3 = 7p/p and thus can be calculated in terms of the 
speed by equations (6) and (?), It is of interest to note 
that the equation of continuity (8) now, unlike the case of 
incompressible flow, becomes dependent on the dynamical equa- 
tions and, consequently, is nonlinear. This change in the 
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character of the fundamental equation makes the direct solu- 
tion of the problem in space coordinates very difficult. 


2. Transformation of the Differential Equations 

The assximption of irrotati onality implies the existence 
of a velocity-potential for such a flow. If this function is 
introduced to eliminate u and v, equations (4) and (8) 
would give, rise 1 to a nonlinear partial differential equation 
of the second order. The problem is further complicated by 
the' possible appearance of supersonic regions, or regions 
where the speed of flow is larger than the local sonic speed. 
This means that for some part of the domain, the equation Is 
of the elliptic type; while in the other part, it is of the 
hyperbol i o type. Thus the equation not or, ly is nonlinear but also 
is of mixed typs, and there is as yet no successful method to 
deal with it directly in the physical plane . Molenbroek (ref- 
erence 6) and Chaplygin (reference 6) made sone progress in 
solving the problem by transforming the equations from the 
physical to the hodograph plane in which u and v are taken 
as the independent variables. If this is done, the differen- 
tial equations become linear and tlms can be solved by well- 
known methods. 

Let the transformation be defined by 


u = u(x,y) (9) 

v = v( x, y ) (10) 


If u and v are continuous functions of x and y 
continuous partial derivatives, and if the Jacobian 



is finite and nonvani shin'g, a unique inverse 


with 

trans- 


formation exists. Under these condi ti ons ,eauat i ons (8) and 
(4) are. easily transformed into 



( 11 ) 


3* _ _ o 

dv ou 


• I 

( 12 ) 
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Corresponding to cp(x,y) in the physical plane»there is in- 
troduced here a function xCu-i -7 ) defined by 

X * xu + yv - cp; ■ x = 4r^» • y s (13) 

- ou dv 

While equation (12) is satisfied identically, equation (ll) 
becomes 


(i - H*') + svu &I>l + A _ zlN i!* = 0 (i4) 

'• c 3 '' 3v 3 c®' 3v3u ^ c s/ 3u 3 


As c is a function of q alone, the equation fer X(u,v) 
iq then. linear. From equation (13) it is recognised -that if 
X(u,v) is known, a one-to-one correspondence between the 
space coordinates and the velocity components can be easily 
e stabli shed . 

However, it is also clear that this function is incon- 
venient for obtaining the streamlines and the flow in the 
physical plane. To solve this part of the problem, a plan 
may be adopted similar to Chaplygin'. 6 by introducing both the 
potential function ?(x,y) and the stream function 'Wx,y) 
defined by: 



ay 


(15) 


pu 


_ a^ „ 

-r-. pv 

o Sy 



(16) 


From these definitions are obtained immediately the following 
equivalent relations: 


dy = udx + vdy 


p^dty = -pvdx + pudy 


(17) 

(18) 


For the subsequent calculations, it was found convenient to 
introduce the polar coordinates in the hodograph plane de- 
fined by: 
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u » q cos 6, v = q sin 8 


(19) 


where 6 is the inclination of the velocity vector to the x- 
axis. Functions cLx and dy can he solved for from equa- 
tions (17) and (18). As dx. and dy are exact differentials 
the conditions of" integrability then give! ■ 


dq p V. . a 3 / q d9 


( 20 ) 


i JSB _ £a 

q 96 p Bq 


( 21 ) 


By eliminating cp between equations (20) and (21),' an equa- 
tion for \|/ is obtained: 


s d 3 ^ 
L 

3q 




aq, . 


5w 

S6 a 


= 0 


(fe2) 


Equation (14) can also be transformed in polar coordi- 
nates. The procedure is straightforward and yields 


a d x 

dq 3 




13 

ae 3 


(23) 


There is an additional relation between X and cp de- 
rived from equation (13): 


cp = qXq - X (24) 

Since cp is connected with 'k, this relation ensures that 
\j/ and X are properly connected and represent the same flow 
pattern in the physical plane. It can be thus considered as 
the equation of compatibility. Bquations (22), (23), and (24) 
are the three fundamental equations in the present problem 
dealing with the t wo-dimensi opal ’ f 1 ow of a compressible fluid. 
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3. The Particular Solutions of the Differential Equations 

As the differential equations for ^(q,8) and X(q.©) 
are linear, a general solution can certainly he built by 
superimposing the particular integrals of the equations. To - 
obtain the particular integrals, let 'Wq.Q) and X(q, 6) 
be of the following forms: 


\J/(q,6) = q u v^ v (q) e ilJ0 
X(q,8) ■ q u Xy(q) e iu6 


where u 
(22) and 


ar e : 


9 


s 


9 


a 


is any real number. By substituting in equations 
(23), the equations satisfied by ^(q) and Xy(q) 


I!** + 

dq S 

(j3u + 1 

c 3; 

|q~^. + 
d q 

u( u + 1) 

a 

9 

c 2 


0 (25) 

d % + 
dq 2 

(2v + 1 

3 X 
C S > 

1 9. “ — + 
dq 

1>( V - 1 ) 

a 

9 

c s 

*u = 

0 (26) 


Now each of these equations can be further reduced by changing 
the independent variable. The appropriate transformation is 
found to be 


T 



with 5 = 7— i 

Y _ 1 


By expanding the gas to zero pressure, or vacuum, the maximum 
velocity is obtained. Equation (6) shows that the maximum 


speed is q„ = / — -s — e . Therefore, the maximum value of 
e u max v Y - 1 0 

T is unity, Similarly, it is found that for the speed of 


the flow equal to the local sonic speed, t 
tions (25) and (26) then become 


25 + l’ 


equa- 


T ( 1 - T)^' 1 (t) + 



(a v + b 


u 


1)t 


NV 


( t) 


= 0 
(27) 



18 


NAOA TN No. 995 


T(l r T)X^"(T) + c p . - (ap + p + b v + p + 1 )t]xp'(t) 

L _j 

. ■ . - Up + P ) ( b p + P)Xp(T) * 0 (28) 

where 

a p + ^p * u - P . a y bp * - \ PuCu + l), and Op = v + 1 (29) 

'O 

These are the hyper geometric equations, of which equa- 
tion (27) was first obtained by Chaplygin in 1904. (See ref- 
erence 6.) The differential equation of this type has three 
regular singularities at 0, 1, and +co . If the differences of 
the two exponents at the respective singularities; namely, 
c - 1 , a - b , a + b - c , are not integers or zero, the two 
fundamental independent solutions are F(a,b; c; t) and 
T i-c + a *- e, 1 + b - c; 2 - c; T ) . They are single-val- 

ued and regular in the whole plane with a cut from +1 to +co . 
The function’ F(a,b; c; T) known as the hypergeometric func- 
tion of- general' para.met er s' a,- b, and c, is defined by the 
hyper ge omet ric series which is absolutely and uniformly con- 
vergent when I T | <1, provided Kl(c - a - b) > 0. For 
lT| > 1, analytic continuation has to be used. Furthermore, 
it- is normalized, so that at . T a 0 


. F(a,b{. c; 0) = 1 (30) 

Hence , the particular, solutions of equation (27) are 

FUp,bp..;Cp;T) , t 1 "' c u F ( 1 + a v -Cp, 1 + bp - c v ; 2 - c v ;t) 

(33.) 

The particular solutions of eauations (28) are 

F(ap + P; bp + P; c y ; t), T 1 ”^ F(l + a y > p - c y , l‘+-bp 

' ~ ‘ 1 ...+ P - Cpj 2 - Cpj- t) (32) 

Here a.,, b,,, and c,, are parameters defined by equation 
(29), U 
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When u Isa positive integer while and "by re- 

main as they are, the second integral will reduce to a con- 
stant multiple of the first one. This case was first studied 
"by Gauss (reference 14), who found a second integral involv- 
ing a logarithmic term "by considering the limiting value of 
the integrals given as v tends to an integral value. The 
method has been further developed by Tannery (reference 15) 
and Goursat (reference 16). However, the .form regarded as 
conventional nowadays was that obtained by Frobenius 1 general 
method. According to this method, the' "pair- of • f-uudajman-tal- 
solutions of a hypergeometric equation are 


FU.b; n + 1; t), K a T~ a -jV 1 F(a,b; n + 1; T)log T 

+ T n Q n ^(a,b; t) + (33) 

when c n = n + 1, n “being a positive integer; and 

00 

Q k)(a a bj t) = V -.P.(a + m) . m ) T m 

n r(a)T(b) 4-r(m. + l)r(n + 1 + m) . 




n-i 


(-l ) 13 * 1 ^( n + l) y (L n« r(a - n + m)r(b - a + 
r(a)T(b) L F(m + 1) 


m)F(u - m) m 

■— — T 


m - 1 

$(a,bj m) * y 


m 


-I- + -JL 


; L a + r b+r n + 1 + rJ 
r = o 


-y ± 


r = i 


Here a,b may be either a n , b n or a n + 5 , b n + p de- 
fined in equation (29) according to whether the system (33) 
is referred to as solutions of equation (27) or (28). And 
K a can be determined so that the product of the second inte- 
gral and q Sn satisfies the condition (30) . 


In view of the fact that the second integral in (33) 
does not constitute a family of solutions with the second in- 
tegral given in (3l) or (32), it. is very desirable to define 
a new function as second integral which will be continuous in 
u as well as in T. Let JL.(t) denote the first integral 
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I’Ca.'b; Cy ; t). As a second integral, take the- linear combi- 
nation of the solutions: 

^(t) = Ky |r(i - eu >r(a)r(b)iy(T) + rd - c v )rd + a 

- Cyjrci + b - c u )t 1 '*° U £- v (t)J (35) 

where 

X_ u (t) » F(1 + a - c v , 1 + b - c v ; 2 - Cy{ t) 

This is evidently a solution and valid for all values, of u. 
The constant Ky is determined subject to the following con- 
dition: 

q 2U J* V (T) = 1 for T * 0 _ (36) 

The value of Ky then is found to be 

K“ l = (2pc 0 ) U T( C y -l)F(l,+ a *. Cy)r(l + b - c y ) 

Using the relation 

P(c V )r(l - Op) = IT CSC C^TT 


equation (36), when multiplied by q 2U , will define a new 
function G v (t): a, b^-n 


1 U <T) 


TT 

Sin CyTT 


JP 


( c-y )r ( cp - 


ru)r(b) t v Xy( T ) 

i)rd + a - Cy)r(i + b - 



r(c v - i)r(2 - Cy). 


(37) 


When u takes integral values, the expression in the 
bracket vanishes; however, the limit of the ratio exists: 
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The usual definition of the limit of a quotient gives 


(T) = (-1) 


n+i 




r(a)r(h) T u 1»(T) 


- 9 l-u( t) 


r(u + i)p( v )r(a - v)ro> - v) ^p(«)r(i - u)J 


V an 


By considering separately the first n terms in X (t), as 

r(l - ) has poles at V a n, a straightforward reduction 

yields: 

£„< T > * ° n lo e T I n (T) + rU ^i s> < T > + p n-i s) <^ < 39 > 


where 

q, ^ (T) a 

n r(n)F( 


(-i ) n+1 v 1 r 

— rrr r ) v|r(a + m) + \{r(b + m) 

7 (-n + a)T(~n + b) — > L 

m=o 


- \Kc n + m) - M/(a + 1)1 r(a + m)3?(b + 

J r(c. + + 


m)_ T m 
1) 


p. - 

n- 1 


n-i 


r(n)r(a-n)r(b-n) <- 

mao 


V ( i) m P(a-n+m)P (b-n+in)P (n-m) T m 

r (m + 1) 


ULita)r(b) 


r(n)P(n + l)F(a - n)D(b - n) 


and \|/( | ) denotes the derivative of log F(|). It can "be 
seen that the difference between (33) and (39) lies only In a 
constant multiple of the first integral which has been absorbed 
in Q n < 3 )( T). 

In the following discussions, the two fundamental solu- 
tions of the hypergeometric differential equation will be 
taken as and q-au £ (t). The normalization condi- 

tions given by (30) and (36) are chosen for the continuous 
passage of a compressible to an incompressible flow. Ultimate- 
ly, these functions are again defined in terms of power series 
which are absolutely and uniformly convergent within the do- 
main t t j < 1. However, since the maximum value of T at- 
tainable by the fluid is unity, the continuation of the solu- 
tions beyond the unit ci-rcle will not be discussed here. 
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Thus £y(T) and 2 ( T ) denote the two independent 

integrals of equation (27) where p is any positive number. 
The particular solutions of equation (22) are then: 


q u Iy(T) 


M 


cos 


pd + A 


P 


(a) 


sin pd 


q~ v G (T) 

P 


By^^cos pd 


+ B„ (s ^ sin yd] (40) 


where By^ 1 ^, and B y ^ 8 ^ are constants. Sim 

ilarly, those of equation (23) are 


[r,d 


cos pd + Ay^ 2 ^ gin Pd 


. q" V 5p( T ) 


By^ l ^cos Pd 


+ B 2 ^ sin pd 


1 


(41) 


where J! p (t) and q~ au G^(T) are the two independent inte- 
grals of equation (28) and Ay(*), Ay(®), By( 1 ), and By(®) 
are constants. 

In addition to these solutions, there are two other inte- 
grals each of which is 4 function of one variable only. 
Assuming >4/ * ty(q) or ^(d), then equations (22) and (23) 
yield, respectively: 


c 1 d 


and 


and 


(1 - 

y J P d T 

(42) 

1 

(1 - 

T)- 9 il 

(43) 


which correspond to the fundamental solution of the Laplace 
equation. 

As c Q approaches infinity, all these particular solu- 
tions reduce to the familiar harmonic functions: namely. 


“p 


(O 


cos Pd + A 


V 


(a) 


sin pd 


> , q" v B, 


(l) 


cos 


pd 


By (2 > 


sin pd 


(44) 
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S3 


and 

ci&, c 8 log (i (45) 

This property which is the consequence of (30) and (36) is 
essential in the method presented in this report for connect- 
ing a compressible flow with the Incompressible flow of simi- 
lar configuration. 

In the subsequent calculations, another integral will be 
encountered for the function x(q»d) which corresponds to 
the imaginary part of w log w e*™ or q log q sin d 
-qC^T -* d) cos d of the ineompre ssihle-- flow.- -• Suppose the so- 
lution possesses the form: 

X(q,d) = X x ( q) sin d - X^qK-rr - d) cos d (46) 

By substituting the expression in equation (23), X a and X s 
.. are found to satisfy simultaneously the following differep- 
- t.ial equations: 

, s x 1 «(q) + (1 - si) (dV - X,) . a (1 - sf)x s i* 7 ) 

c 0 „ 

: . ^ q a x 4 n + (1 - ^) - x a ) . 0 us) 

Equation (48) can be easily integrated by putting X 3 = ak a (q) . 
The condition that X a -5> q as e 0 — > 00 requires k s (q) to 

be a constant,. The second integral of equation (48) is just 
' 'the second of (43) which, in the limit, tends to log a. Thus 
X 3 «= q- is the appropriate solution. With this solution 1 . It 
1 is possible to proceed to solve equation (47) by aeusimiag- 
X'i = qkj(q) . The equation for ki(q) is again integrable by 
quadrature, and the result is 

T 

ki (q) = ^ - -f (28 + 1) log T - I •+ Ki / (1 - T)" P H + E s (49) 

2(g + 1) L T J T a J 

where K l and. K a 'are the constants of integration. Hence, 
the desired particular integral is 
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X(q t £) = qki(T) sin d - q(n - d) cos $ (50) 


The correspondence between solutions for compressible flow 
and for incompressible flow is summarized in table 1, 


4. The Properties of the Hyper geometri c Junctions 

of Large Order 

The behavior of E^a^.b^} c^; t) for large positive 

values of v has been discussed by Chaplygin in connection 
with the question of convergence of his series solution for 
the flow of a gas jet . However, his discussions are limited 
to the subsonic flow and, for this reason, the value of T 

is restricted to the interval 0 < T < 1 , In the more 

~ - 28 + 1 . .. ... 

general problem where both subsonic and supersonic flow may 
exist, the whole interval 0 < T < 1 has to be considered. 
Furthermore, both integrals of the hypergeometric equation 
are involved, as will be shown in part II, As a preparation 
for the proof of the convergence of the solutions, the prop- 
erties of the hypergeometric functions of large order in the 
extended interval will be discussed presently. 

Chaplygin (reference 6) introduced a new function 

p 

— L,(t) defined as the logarithmic derivative of J (t): 
2T u P 

namely, 

p| (t) = 2t~ log T2 2 (i), v ^ 0 (51) 

p dT P 


where IT ( T ) ' denotes the first integral of the hypergeomet- 
ric equation (27) or (28) and p can be either an integer or 
not an integer. Then in the place of equation (27) or (28), 
the corresponding differential equation for is a Riecati 

equation: 


x(£p) = 




1 (25 t 1) 

1 - T 


0 (52) 


where the lower sign corresponds to eqpat ion. ( 28 ) . As shown 
by Chaplygin, 2y( T )» although an oscillatory function, - - 
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can have no root in 0 < T < — : — ± 

“ ” 2p + 1 

is finite and continuous in the same 
can he deduced also that fy(0) = !• 


and, consequently, 

interval. Moreover 
and |y'(0) = -p. 


, it 
Since 


|„'(t) does not change sign in 0 < T < i ; |„(0 is 

monotonic decreasing and eventually vanishes at t q < T*, T* 

being the first root of the hyperge ometri c function for 
V > 0. Since T* is a decreasing function of u, when v 
becomes large, t* and consequently T Q will differ from 

1 by a small quantity. 

2(3 + 1 


Chaplygin's theorem ,- In 0 < T < -- if> a monotonia 

continuous function Tty(T) satisfies (i) 7]y(o) = 1 and 

(ii) X(T|y) ^ 0, then 


V T > < v > 1 


(53) 


The proof is given in Chaplygin's paper (reference 6,) In 
the case of the second integral 3?y(T), the theorem is Btill 
true with the signs of inequalities reverBei because it can 
be verified that X ( = 0, where l.ytT) corresponds to 
the case of F u (t) instead of 2.\>( T ) in (5l), and |_y(0)*-l 
therefore £_ n ,(T) is negative in 0 < T < — 

° — — Oft j. 1 


Corollary ( 51 ) . - In 
ly(T) and Gy(T) fall 


0 < T < t , the functions 

2P + 1 

respectively between the limits! 


(i> 

(ii) 

where 

Tj, (t ) 

®s(f ) 


T l U (T) < Fy(T) < T/ (t) 

T x U (t) > Gy(T) > T 2 U (t), u > 1 




+ ' 1 ) T ~ dT 

T J 2TJ 



(54) 

(55) 


(56) 


( 57 ) 
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‘This can be verified easily by choosing Tl^ to be 

/l' - ( £6 + 1 )t or a . T) P. As, evidently, in 0 < T < — - — 
V 1 — T “ ‘26+1 

when 

„> i, /IXISZIE< t,< a - tv* (58) 

and 

> > - (1 - T) * U9) 

and furthermore , X(T|^j) ^ 0 are satisfied, consequently, it 

follows the results. 

Corollary (52').- In 0 < T < , the absolute value 

*. - 26+1 

of the logarithmic derivative of f(a v ,b-y; a v \ t) divided by 
u, is* bounded both above and below - that is, 

M (t) < .+ .i*-! — 1 1 < M a (T) (60) 

S'(ay,byj Cyl i*) 

where M 1 (t) and M 2 (t) are independent of u. This really 
is a consequence of (58) and (59). 

It shall be noted that t^e results established in the 
foregoing are applicable to JLy(T) * ^(a^ + + P* c pS T ^» 

provided v is large, because then the two equations (2?) and 
(28) tend to be the sape. 

Obviously, Chaplygin's theorem ceases to be true when 

T > — . For in the interval - < T < 1, the solu- 

26 + 1 26 +1 

tions of the hyper ge omet ri c equation are oscillatory and, 

hence, within any closed interval ifc — — < T < 1 the num- 

26+1 

ber of roots of J? v (t) will be proportional to p. (See ref~ 
erence 17.) When p is large, there will be a large number 
of roots in the interval. As a consequence, the function 
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^(t) will have there an ever increasing number of simple 

poles, and a finite interval in wh-ich ! u (' r ) remains finite 
for all u does not exist. 

To carry the investigation over into < T < 1, 

2{3 + 1 

the method is modified. Let 3 ^(t) and T^t) be two inde- 
pendent solutions of equation (3?) or (28); and let the lin- 
ear combination be denoted by 

3V< T ) = X V (T) + i F p (T) (61) 

The complex function is, of course, a solution of the same 
differential equation. In terms of its modulus E^(t) and 
argument 4 >^(t), the function may also be expressed as 

i { d) f t) 

^•(t) = R v (t) e u (62) 

where both R v ( t) and 4> p ( t) are continuous functions with 
continuous derivatives. By comparing with (61), it is neces- 
sary to have 

F v (t) = R v ( t) cos 4> v ( T ) (63) 

ff v (T) = R V (T) gin <p v (T ) (64) 

According to the Sturm separation theorem, F^t) and S’ u (t) 
never vanish simultaneously in any closed interval and R-^t) 

never vanishes in i-, — < T < 1 and remains positive in 

28 + 1 

the whole interval. Then corresponding to (51), a complex 
function I v *(t) can be defined as follows: 

u !u*(t) = 2tA- l0 § t1 5 V*( t > <65) 

which satisfies the same equation (52). On separating into 
real and imaginary parts, the Riccati equation for | u *(t) 
becomes 
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: (t (i) p (2) N \=f , (i) + ! — : | (* ) + JL 

,& u 7 “ S i ~ r l v 27 


:3(l) fes(s) 

’U b P 


(8P + I )t - 1 
1 - T 


a 0 (66) 


X 2 (lv U) ,t„< 1 )) S * --S— l v ' a > + * l„ (s) E v (l) = o (67) 


where i and i are real continuous functions of T 

* V 0 


defined as 


£*(t) « + i £y (s >(T) 


( 68 ) 


Their connection with E^t) and ^(t) separately are giv- 
en "by means of (65):. namely. 


v |,^ 1 ^(t) * 2T-f- log t'S R (t) 
u dT v 

v tj‘h T) ■ 2t± V 0 


(69) 

(70) 


(i) 


Now equation (67) can he integrated in terms of '(t) 

and whence |y( s )(T) can be eliminated from equation (66). 
Then the equations for . and ^ 


’V 


are 


*»0» (1) ) s 5'v Cl> * “ t„ (l) + £ [fl U) 


- to (1 


T , V dT 

- T) 8p e“ 2 V T + » 0 


■ 1 - T 


(71) 


4 < T) 


■M 1 - t) 


6 „ r t (l) 51 

^ Q “ U / *V T 


( T <>* 


( 72 ) 
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Equation (7l) together with the condition fp^^Co) = -1 de- 
termines uniquely the solution tp^ 1 ^^). The actual value 
of | (^(t) can he expressed, of course, in terms of the 
known hypergeometr.ic functions. But the problem on hand is 
to determine the properties of for large u which 

are given by the, following theorem. 

Theorem (53) .- If Tlp^^T) is continuous and monotonic 
in To < T < 1 and satisfies X^Tju^ 1 )) < 0, then for all 
v > N _ 

T] v . (l) (' r ) > l V (l) (T) (73) 


The proof is given in appendix A. 

Ocr ol larv ( 53) . - In T 0 < T < 1, the following inequal- 
ity holds for the modulus of ]?*(t): 

v/s 

M t )/M t 0 ) < (If) . V > N (74) 

where 

(2P + lVr 0 - 1 > 0 

Tor in T 0 < T < 1 f 0; and hence T|p ^ 1 ( T ) - 0 

satisfies the condition 0 > tp^ l ^(T), which gives (74) by 
int egrat i on . 

Now, since ) (T) •is bounded by zero for all V ^ 0 

in T 0 < T < 1, it is implied also that 

R V (T) < ? 3 U(t) (75) 

where T 3 (t) a ~~r/a • ^ ere the constant t 0 can be determined 

by joining T 3 at T = t 0 with T x or Tg defined by 
eqjF.'. J.o.on (56) and (57). Then from equations (63) and (64) 
it follows that for v> > N 


!p( T )' 

G- (t) 

—v 


< 2 3 u (t) 

< V u (t), 


Tq < T < 1 


(76) 


(77) 
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; PARE II 


COHSERUCE ICH OP EHE SOLUEIOHS FOR 
COMPRESSIBLE PLOW AROUND A BODY 
5. Chaplygin’s Procedure 


In the previous sections, the particular solutions of 
the differential equations in the hodograph plane are ob— 
tained. Since the differential equations in the hodograph 
plane are linear, superposition of solutions is allowed. 

In other words, if these particular solutions are multiplied 
by different constants and then added together, the sum is 
again a solution of the differential equations. By this 
procedure, general solutions can be constructed from the 
particular solutions. 

However, there is a d iff i cult y in such a method of 
constructing the general solution — the difficulty of making * 

a proper choice of the multiplying constants for the partic- 
ular solutions so that the resultant solution will give a 
flow satisfying the boundary conditions specified in the . . 

physical plane. This can be seen from the fact that the 
space coordinates x and y are obtained from x which is 
not explicitly connected with the stream function. In 

fact, to obtain the coordinate x and y directly from 
would involve an integration in the hodograph plane, Ehus 
the linearization of differential equations in the hodograph 
plane i3 obtained at the expense of the simplicity in boundary 
value problem. Eo guarantee that ^ and X do actually be- 
long to the same flow in the physical plane, an additional 
condition besides the differential equations for 'l' and X 
has to be satisfied. This condition will be discussed in sec- 
tion 11. 

Chaplygin (reference 6) suggested an ingenious method 
of solving this difficulty by using the well-known solutions 
of the incompressible flow, Ehe first step in this method is 
to find the incompressible flow around a body ’’similar 1 ’ to 
the body concerned. (Ehe meaning of the word "similar" will 
be made clear in the following paragraph.) 

Ehe stream function , for instance, is then expressed 
in terms of the speed q and the inclination 6. Ehe function 
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( q. , 0 ) can be expanded into an infinite series each term, 
of which is of the form q. n cos n£ or q. n sin nd, For the 
flow around a body with constant velocity U at infinity, 
the function %(<1, d) has a singularity at the point q * U, 

9 * 0 in the hodograph plane, since there all the stream- 
lines, or lines of constant W 0 originate. Ihus, there- are 
two forms of the series expansion of One is convergent 

within the circle q = U{ while the other is convergent out- 
side of the circle q = U. The first, or '•inside, 11 series 
must be regular at the origin of the hodograph plane. There- 
fore, only positive values of the integers n can occur. 

The second, or "outside," series can have both positive and 
negative v. Chaplygin's method is to use the inside series 
for as the starting point for obtaining the desired 

solution if for the compressible fluid. He suggested 
choosing the multiplying coefficient of the particular sol- 
utions for the compressible flow by the condition that for 
the limiting case of infinite sonic speed, or incompressible 
fluid, .the series will degenerate to the inside series of 
the incompressible flow already obtained. The series for 
the compressible stream function so constructed can be 

called as the inside series of. The outside series for 

i 1 then can be obtained by the method of analytical continu- 
ation with the aid of the "outside series" of the incompress- 
ible flow. 

These -s olutions so constructed for the compressible flow 
contain the Mach number of the undisturbed flow as a parameter. 
They constitute a family of singly infinite solutions. In- 
cluded in this family of solutions is the limiting case of 
zero Mach number of the free stream. This limiting case will 
give the incompressible flow around a body used as the starting 
point of this method. For other values of the free-stream Mach 
number, the body contour is generally different from that corre- 
sponding to zero Mach number. Thus, if the compressible flow 
around a given body is desired, the body shape for the initial 
incompressible flow must be slightly different from the given 
body shape. However, if a geometric parameter is included in 
the solution, such an adjustment is not difficult to make. 

It may be stated here that owing to the regularity of the 
solution at the origin of the hodograph plane, only the first 
solution of the hypergeometr ic differential equation appears 
in the inside series. For the outside series, both the first 
and the second solution of the hyper ge ome tr ic differential 
equation are necessary. This is in direct analogy with the 
appearance of both positive and negative exponents of q in 
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the 'incompressible outside series. This fact is particularly 
important, since previous investigators seem to be unaware of 
it. /Chaplygin himself d.id not use the second solution of the 
hypergeometric differential equation, but that is simply be- 
cause, for his problem, there is no singularity in the hodo— 
graph plane and hence only the inside series is needed, 

6, The Functions for Incompressible Flow 

Following the procedure outlined in the previous section, 
the analysis starts with the functions required in defining an 
i rrotat ional incompressible flow. For. this case, the sonic 
speed o Q tends to infinity, and the equations fob the veloc- 
ity potential cp 0 (x,y) and the. stream function ^(^y) all 

became harmonic: 

A<p 0 = 0 

^o = 0 

where A stands for the laplacian operator, 
the complex potential, it can be shown that 

W 0 (z) = <P 0 + 1 ^o (80) 

where 

z = x + i y 

If w denotes the complex velocity u— i v, it is connected 
with W q (z) by 

w = s w ( z ) (81) 

dz 

If w»(z)=^0, it always is possible to solve 1 for z in terms 
of w; namely, 



z = z 0 (w) 


(82) 
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In general, this solution is not s ingle— valued and will "be 
discussed later,, By introducing this relation into equation 
(80), the complex potential function in the hodograph plane 
can he obtained 


W (w) = cp <u,v) + i ^ (u,v) (83) 

y o o 

In case equation (82) is many— valued , this would correspond 
to one branch of the function. 

It is clear that in this case X Q (u,v) is also a 
harmonic function. Let cr Q ( u,v) be the conjugate function 
defined by 


Hence 


dX 0 _ _ dcr 0 
ou 6v 


(84) 


d*o _ 5cr 0 
c)v cu 


(85) 


A 0 (w) = X, - i 


( 86 ) 


where 


w * u — i v 


Thus A q (w) is an analytic function of w. Prom equation 

(13) the derivative of A 0 (w) with respect to w must be 
z. That is, 


aA 0 

dw 


z 0 (w) 


But z 0 (w) already has been found from equation (82). 
Therefore, 


A 0 (w) 



z 0 (w)dw 


+ 


constant 


(87) 
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The real part of A (w) gives x 0 ( u » v 5 as required, accord- 
ing t o ( 86 ) . 


7. Conformal Mapping of Incompressible Plow 
on the Holograph Plane 

Before the construction of solutions for the compressi- 
ble flow, the general character of the solutions in the 
holograph plane should be examined. This can be done easily 
by investigating the behavior of the transition function 
z (w) for. an incompressible fluid. To start with the 
simplest case first, consider a steady irrotational flow in 
an infinite, simply connected domain D bounded by a curve 
C in the z— plane, with a parallel flow at infinity (fig. l). 
At every point z of B there is one, and only one, velocity 
vector if. If the curve 0 is mapped into C and infinity 
corresponds to a point JP on the axis of reals of w within 
C.« then the domain D is mapped into D by a mapping func— 
t ion 

w ** w( s ) 

defined in (8l), where w(z) is an analytic function of z. 
The inverse function 


8 = Z Q ( w) 

will set up a continuous one— to— one correspondence between 
w— and z— plane, provided the mapping is conformal. This 
requires that w(z) is analytic, simple within B, and 
w 1 ( z) 0. 

However, for most problems these conditions cannot be 
satisfied throughout the field of flow. In the first place, 
the function w(z) is generally nonsimple, for example, in 
the case of a uniform flow, w( z ) = constant, thus w*(z) — 0 
and the whole z— plane corresponds only to a point in the 
w— plane. Furthermore, the complex velocity for a two- 
dimensional boundary— value problem generally can be put in 
the following form: 


W = Wjo +■ w*( z ) 



NACA TN No. 995 


35 


where Woo is a constant. Che "boundary condition requires 
that w*(z) = 0 and, as a consequence, w'*(z) = 0 as z 
"becomes infinite. Therefore, in all cases, the point P 
in the w—plaue, is a singular point. It is a "branch point 
at w OT if z ( w) is many— valued} or a pole, if otherwise. 
In practice, £here are two kinds of singularities that play 
a dominant yole in the problem of two-dimensional flow. 

These singularities will be investigated presently. 

Branch point of order 1 . It may be recalled that, 
when a closed body is present in a uniform flow, there 
always exist two stagnation points both of which correspond 
to the origin of the w— plane. If a streamline PS is fol- 
lowed, for ins tance, ( s ee fig. 2) from +» to S, the por- 
tion SMS' and then to — a curve PS_ in w— plane would 
be described twice. This indicates that the function z 0 (w) 
possesses two branches of Riemann surfaces Joining together 
about the branch point P* In order to make the domain D 
s ingle— valued , a cut Is put along the axis of reals from 
the branch point to + 00 . Then one portion of the z— plane 
is mapped into a definite branch of the Riemann surfaces 
in the w-plane, and this will be defined as the domain 3>, 

If the body is symmetrical with respect to the coordinate 
axes with parallel flow at infinity, then the domain 
DiRlz < 0 will be mapped conformally into D on one branch 
of the Riemann surfaces and D*:RJz >0 on the other, where 
the region within C is excluded. 


Logarithmic sing ul arity .— The flow over a wavy surface 
for instance, placed parallel to a uniform stream has a 
periodic nature. Por such flows there are infinitely many 
points in the physical plane that have the same velocity. 
Hence, there are an infinite number of branches in the w— 
plane, each of which corresponds to a definite portion of 
the z— plane. The function z Q (w) must have a term log 


(l - and 

If, however, 


the point P now is a logarithmic singularity, 
a cut is introduced from the branch point to 




+a> and — rr < arg Ql — ~J < tt, 
made s ingle— valued. 


then the domain D is again 


1 The function z Q (w) is said to have a branch point of 
order k at w = Wco if its inverse w(z) contains the 
part w* which has a zero of order k + X at z * °°. 
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8. Construct ion of a Solution about the Origin 


Prom the considerations of the last 


section, the domain within a circle with radiuB Iwt = PS =» U , 
where U ie the absolute value of w at infinity in z— plane, 
is in all oases s ingle— valued, If a function W 0 (w) is as- 

sociated with a definite flow in z— plane, from section 6 It 
is an analytic function of w and regular within the circle 
|wt = U. Consequently, the following Taylor expansion exists: 


C£> 

V 0 < w) * ^ A B w , Iw I < U (88) 

n= o 

j a 

where A^s are, in gefteral complex. Since w «= qe"" 
and by (80) the imaginary part of W 0 (w) is equal to in- 
compressible stream function \|> 0 , it can be written as 

CO 

^ 0 (q,*) • I» ^W 0 (w)j' ■ q* 1 ^A n cos np + A^ ^ sinn6V(89) 

n-o 


According to Chaplygin’s procedure, 
compressible solution can be obtained by 
the faction q n in equation (89) by the 

n ^ r ' 

q (t) as shown by (40), The second 


the corresponding 
simply replacing 
corresponding 

integral is ex— 


cluded by the regularity requirement at q = 0, However, 
in order to preserve the proper singularity at the point 
(U,0) in the hodograph plane, the compressible stream 
function \J/ is written as 


00 

q n £^ r ^(T) |a 

n=o 


(O 


(•> 


cos n£ + A Bis nj 
a n 


>} 


(90) 


where 


J< r \r) . i^lL ■= q < u oi) 

P(a .b;c iTi) 


n' n' »’ 



XI AC A TN No. 995 


37 


and T i = — 


U 


2.0 c, 


the value of corresponding to the 


It is seen that if 


then 


(t) -^l due to the nomadizing con- 


free— stream velocity U 
T = T, 0, and P ^ 

“'ll 

dition (30). Thus the solution Is reduced to the incompress- 
ible form. Furthermore, if q U the character of . the sol- 
ution is exactly like that of the Incompressible solution, ' 
Hence, all the specified condii* ions ' are satisfied. Of course, 
for the mixed subsonic and supersonic flow, the free— stream 
Mach number is always less than, unity. Thus < 1/20 + 1. 

For later analysis as given in part III, it is convenient 

\ I* ) 

to write if in a different form. Since F (T ) Is a nurely 
real quantity, a: complex function W(w;t) can be constructed 


as 


W(w;T) = ) A F ^ 1 t ) w n , twl < U 

‘ ■ n ”n 


(93) 


E.-0 


Then, similar to .the relation between eouations (88) and 
( 89 ) , vKq,6) can be taken as the imaginary part of the 
new function' W(w;t). Th^s , 


ty( q , ■& ) = Im 




.(93) 


Transfo rm ed ro t?:.i tial functi on,-- Similarly, it is 
possible to construct another function A(w:,t) defined by 


A( 


w,t, ..y 


~ ■ ~( r ) . . n 

-S. < T > l< > 


Q < U 


(94) 


Uss.o 


In this expression, the coefficients A are obtained from 
the expansion of A (w) for the mcompr ess ible flow (87): 


A. 


| i - 

,( w) = ^ A n ■'■w 31 , | w | < ir 


n=to 


t 


(95) 
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and 


Jr) 

It, (t) 



(96) 


Equation (96) is the result of equation (91) and the equation 
of compatibility given by equation (34). Then the function 
X(q«#) for the compressible flow is given by 


x( q»d) 



( 9 ?) 


The functions W(w;T) and A(w;t) are actually regular 
at the origin and satisfy the imposed conditions. However, 
the following question may be raised: Do the series (92) 

and (94) converge and represent the functions 't'(q^) and 
X(q,£) in the domain of validity ? To settle this question, 
it is necessary to prove the following theorem: 

Theorem (88). If the constants A„ and A are sriven in 
n . , n / ) 

equations (88) and (95), while Jr (t) and E (t) are 

— n n 

defined respectively by equations (9l) and (96), the series 

(92) and (94) are uniformly and absolutely convergent in the 
same domain as those of (88) and (95). The proof is given In 
appendix B, 


9. Analytic Continuation of the Solution 
Branch Point of Order 1 

Strea m function . — As proved in the appendix B, the series 
(92) is absolutely and uniformly convergent and does represent 
a regular function W(w;t) for every T in 0 — t — T- and 

i a 

on the circle of convergence it agrees with W 0 (Ue ), of 
which the Pourier expansion exists? 

CO 

W 0 (Uo“ 18 ) ■= £.A n U n e- ln9 
n=o 


( 98 ) 
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In the present section, it is proposed to continue 
the solution (92) analytically outside>the domain { w i 5. II 
with the initial value given "by equat ion ( 9 8 ) , The domain 
outside | w | < U is generally many— valued. To fix ideas, 
discuss first the case of a branch point of order 1, Gen- 
erally, the function W Q (w) has other singularities besides 
the one at w = II . However, such ' s ingular it ie s lie outside 

the region of interest and thus need not be investigated. 

Let the nearest singularity be given by | w I * T > U, Then, 

the domain to be considered outside }w} = U is an annulus 

with a cut joining the two singularities. The proper repre- 
sentation of' W Q (w) in such a. region which has a branch 
point of order 1 at w = U , is 

tf Q (w) = iw~£ W Q *( w) (99) 

t 

where ,rf 0 ’ |c ( w) is single— valued and regular within the open 
annulus U <| w| < V. Hence, in any closed domain 


U + 8 <|wi<V - 8, .8 


being a small number, there exists a uniformly and absolutely 
convergent series: 

00 

W q *(w) '*■ ^ w n + C n w“ n J (100) 

n=o 

which, on substituting in (99), will give the continuation 
of the Taylor series (88). 


1 For instance, in the problem of the .flow around an 
elliptic cylinder treated in part V, there are two singulari- 
ties of the W 0 function given by equation (280): namely, 
w = 1 and w = l/e 3 . The first singularity corresponds to 
the flow at infinity and is the singularity under discussion. 
The second singularity corresponds to a point inside the circle 
of the £ — plane, the plane of the circular section. Since only 
the flow outside the circle of the plane is of interest 
here, the singularity w= l/c 8, need not be investigated. In 
other words, it is necessary only to expand the W Q function 


in the annular region 1 < 



w 

U 
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The solution for a compressible fluid, which has the 
same character of singularities as' W 0 (w) and is valid 

in the annulus U < |w I < V, can be obtained from (100) 
.by introducing the proper hyper geometric functions corre-^ 
^ spending to each exponent of w. That is: 


CO 

w(°)(w;t)= i £ v .(T) w v+ c n * G v (t)w~ U 

n=o 


(ioi) 


which is the continuation of 


W ( w;t). 


Here V 



n being a positive integer; Z y ( T ) and q SV G^(t) are 

the first and second integrals of the hypergeometric equa- 
tion; and B * and C * are constants. It should be 
n n. 

noticed that the coefficients B n * and C n * in the outside 
series for the compressible flow are not equal to B n and 
C n in equation (100) for the outside series of the incom- 
pressible flow. The outside series of the incompressible 
flow is used only to give the proper form of W'- 0 '(w;t) 
for the desired branch point characteristics; while the 

exact determination of tf ° (w;t) has to be made by the 
conditions of continuity, which will be discussed presently. 


Since the partial differential equation considered here 
is of the second order, to ensure that w°^(w;t) is the 

analytic continuation of W^^(w;t), two conditions have 
to be imposed at the boundary of the respective regions of 
convergence; that is, the circle q = U, These two condi- 
tions are the following; 

-«/ ^ (Ue~ i6 ; t x ) =* (Ue"* 10 ; t x ) (102) 


ff W^)(w;t) 

» 

“ v/ °)( w;t) 

Ldq 

1 - ' 1 

T=T 

j 


( 103 ) 


On account of equations (103) and (103), there are two 
relations which have the imaginary parts? 
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4:1 


CO 

X [ b : !,!’:) ^ 

n=o u 

oo 

+ C£ a v (T l ) U~ v j cos^ = - X A^ U s in n 


0 < * < 2n 


CO 

r-i f 


Z_i L 

n=o 


) I (vr > (r)+ (T x )) 

OS d -v (-»o b (t i )+ 2 t x O v '.(T 1 ))] 


COB V 


= ~ \ A U ( n+ St, — ) 

L * V 1 Jj T ,) ' 


s in n' 


n=o 


— a' l 


Here the prime denotes differentiation with respect to T, 
Evidently, the coefficients on .the left-hand side can "be 
solved for in terms of the known constants A n . They are: 


jn= o 


U m ( -i- + -±-') ( 104} 
\m+ v n— v/ 


B* U u (vr v (T i )+ st i i u ‘(t 1 ))+ o* 


U V (“V0 v ( Tl ) + ST 1 G U '(T 1 )) 



^ m A m U S m (T l ) 

n = O 


(105) 


From these two equations, the constants B* and can 

he uniquely determined, provided the determinant ) 

does not vanish. These results are: 



42 


NACA TN No. 995 


Gh>( T i) 

B £ U * - 'Z~ T 7 " TP 


—u 3? 


00 

V A^u Y-J- ■+■■■■ * - Y m£ (t )-v| (t (106) 
2Un(l-T 1 ) H Z-< o V m+U m “ u A m 1 ~ v 1 / 

CO • 

liLg Y A u Y 1 ■■ » ■ A -. Y nl (t )~i>i (T )") (107) 
( 1— T 1 )P j£Z J 0 01 V n+V m 1 V 1 / 


as the Wronskian AdLyi^y) = — ■- q 2U ( 1— T ) =£ 0 and 

ly( t) is defined in (51). 

The solution is again formal. To prove that the 
function W(w;T) is a regular function in the annulus 
region, the truth of the following theorem must "be first 
demonstrated. (See appendix 0.) 


Theorem (9 8) . If the constants Bg and C* are determined 

according to (102) and (105) and if the series ( 100) con- 
verges uniformly and absolutely in a closed domain 

U, + 6 < 1 w I < V — ' 8, then the series ( 101) will converge 
uniformly and absolutely in the domain U+ 6 < I w I < V — 8, 
5' > 0. 


Transformed potential function .— By a similar procedure, 
the continuation of (94) is 



i 


X [ S n + «S 1„< T > « ”] 

n=o ■ 


(108) 


where ? y ( t) and G y ( T ) are the first and 
of equation (28) and the constants B£ and 
similarly determined; namel.y, 


second integrals 
0* can be 


ff* U p 

n 

T i) 

2ttv( 1-t i )~ P 


~ (r) 


I (Ti) 

(109) 


n=o 
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U 


v 


Ivjll l 

2tt1>( 1-T 1 )~& 



( 110 ) 


The solution determined so far ‘is understood to "be the 
principal branch of the function W(w»t), It was assumed 

that the flow at infinity is parallel tc the x— axis. If, in 
addition, the body is symmetrical with respect to the co- 
ordinate axes, the expression for the second branch of 

v/°\w;t) will be identical. However, in a more general 
case where asymmetry exists, the two branches will require 
separate consideration. 


10. Continuation — Logarithmic Singularity 

Stream function .— Consider now the second important 
type of singularity; it is assumed here that the only 
singularity possessed by the function W Q (w) in the finite 
part of the w— plane is a logarithmic branch point at w = U 
about which infinitely many Hiemann surfaces are joined. 

By analogy with (99), tf o (w) bow can be conveniently written 
as 


W Q (w) = Wj(w) + W 0 (w) 


( 111 ) 


where W*(w) is a regular function in the entire domain with 
possibly an essential singularity at infinity, and hence 
generally is given by a Taylor Series or a polynomial in w, 

and W Q (w) = 5 (q,d) + i^ (q,#) is an analytic function 
w o o 

which characterizes the singularity of W Q ( w) , Thus, aside 
from a constant factor, • 


W ( w ) = 
o 


7 l0 * ( l ~ S) 


( 112 ) 
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If a out is laid from 



is restricted 


+U 

in 


the function W Q ( w) will 
cut plane. 


to +® and the argument of 
- tt < arg ^1 - 
he s ingl e— ‘ valued in the whole 


w 

« ) < "• 


then 


The Question of constructing a solution for the com- 
pressible fluid consists, therefore, of two parts; W£(w) 
and W Q (w). However, the construction for W£(w) is, in 
principle, exactly the same as that of (92) and hence only 
(w) will he considered. First, let ^(w) he developed 

into power series in the respective domains of validity. 

The imaginary parts axe; 


~ (i) 

( q • d) 


00 n 

~ cos n q < U 

n = i 


(113) 


\j/ ^ \ q , S) 



q>U 


(114) 


The corresponding expression for 'i' ( q , i) , accordingly, will 
he; 

CO 

J (l) (q,^) = X AaIn(T) (u) 003 q < U (115) 


B j 

(1- 

t > I T U 2 °»-n (Tl 

(§)' 

~n 

cos n$, q> U 

T i 


n=i 


(116) 

where 3^(1*) stands for 

a J1 1 * c n > ^ ^ 

and 

G n (T) is 

defined hy (39). 





The function 

W ( w) 

may he regarded 

as 

the complex 

potential of a complex source situated at 

w 

= U. It is 

known that in this 

cas e 

the normal derivative 

of ^(q.d) 
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on |wj = U is a constant, except at w = U, where it 
becomes infinite. This "boundary value can "be expanded 
uniquely: 

CO 

— c o S n -6 = ”, d 0 ( 11*7 ) 

n*=i 

The corresponding problem in the case of compressible i^low 
can be put in an analogous manner. to find a function 
vKq,*) which is continuous together with continuous partial 
derivatives and the normal derivative of which on |wi=U 
is constant. Thus, the conditions (102) and (103) in con- 
duction with equation (11?) demand: 


F 


u <n> A n - s n (Tj) o n = 


(118) 


[n ^ ( r ^ ) + 8T, l. n 

+ fn 0 (T ) - 2 T, 0 * (t ) C =4B(l-T / (119) 
L — n l 1 n 1 J 1 


where the constant B can be determined when the normal 
derivative xj/^q.e) on | w | = U is assigned. By solving 

equations (118) and (119) and using the relation of the 
Wronskian of the two independent integrals of equation 
(27), there is obtained 

A = — B ft ( T ,) (120) 

n n n 1 


C 


n 





( 121 ) 


Thus the function ^( q , 8 ) is completely determined. 

Transformed potential function .— The associated function 

%(q,8) can be similarly constructed. As A Q (w) is derived 
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from (87) by integration of the inverse mapping function, 

it must involve a term ^1 — log which represents 

the singularity of the function A c (w)» As in equation (ill), 
A q (w) is split again into two part si 


A q ( w ) = A*( w ) + A q (w) 


(128) 


where A*(w) is an entire function and A 0 (w) is 


A 0 («) 


(x - §) 1.. (l - g) 

M) 


(123) 


Now the solution corresponding to log ( 1 — — ) can he 

determined in exactly the same manner except that the 

hypergeometr io functions Involved are F (t) and G- (t) 

“H n 

instead of J! ( T ) and G (t). The part that will require 


special consideration is the term — log 


o - «)■ 


he denoted hy \ o (w) 


o o 


Let it 


= 


- i * log 
i U 


( l "S. 


(124) 


This function is also multiple— valued. Let the argument of 
^1 — again he restricted in — tt < arg — ~J < rr * th< 


in the cut plane the result will he 

CO 

hi t 
O i 

n=i 


a i V i 

i A n VU, 


n+i 

v?; • 


r( o ) 


1 w 

i L~ U 


, w in , 
log - e + 


I ; ©■ 

n=i 


-n+i 


| w I <’ U 


. jwl > U 


(125) 

(126) 
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According to equation (-86), the function X 0 (q«d) is 
defined as the real part of A 0 (w). That part represented 
by equations (125) and (126) is then 

CO 

Xq^U.S) = - ^ (u) sin (127) 

n=2 

Xo°^(l»$) = § log - sin 4 - ^ (tt— #) cos * 

co 

+ ^ ~1_ C~^) sin n-d (128) 

n+1 VU/ 

n=i 

The particular solution corresponding to 


■2- log ^ ein d — ^ (ti— d) cos d 
U U U 


already has been given in equation (50). Hence the solution 
for the compressible flow is 


CO 

J (i) (q,4) = - ^ A n f Q (T) 0N sin n * (129) 

n=a 


X^°^(q.^) = ^ k(T) sin £ — 1 (it— d) cos 


) °n f-') n sin nd (ISO) 


n=i. 


where 


k (T)= — [(2P+D iog -(“fV f (1-T) P il](!3l) 

2(6 + 1) L T l V T V T J T J 
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The conditions (102) and (103) together with an expansion 


1 

2 


sin d + 


/ 1 


n=s v 


n+1 n—1 


sin n £ = ( rr— d ) cosd, 0 < ^ < 2tt 


require that: 


£ (T ) A + G ( T) C~ = -1- + -i- (132) 

— n 1 n — n 1 a n+1 n—1 


[n (v + 3 t 4 t i ) j r n 


4( T .)* ST > -a^V 


n n+i n—1 


+ ~-i n f 1 (133) 


and 


£i(T ) o x = i 


(134) 


[. “ Jl( T i) + 2 T i »I(T X )J 0, + 2T j k'( T x ) = n = 1 ( 136) 

By solving (132) and (133) for A and 0 , there is 

obtained: n 


t \ M v 

^n = T 1_n ^ ( O ) ^ ( T x ) (136) 

“ n 2 -l ■ \ • ~ T - n 1 / n 


C (1- n? (T )) F (tJ, n f 1 (137) 

n 3 _l n ’" I1 x 

by using the Wronskian of the independent integrals of equa- 
tion (28). With Cj given by (134), the constant Kj can 
be solved for from Cl35); it is 
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, K x = - (l-T ) P [l + pTj. + (p+l)T 8 - - --- - - 1 (138) 

L ii( T i) J 

M ** 

The solutions \]*( q , -6 ) and x(q.,d) in the whole domain 
under consideration are uniquely determined. Since the 
dominant properties cf the hyper ge ome tr ic functions discussed 
in section 4 hold, in general, the equation of convergence 
can be similarly settled. 


11. Transition to Physical Plane 

—In- th-e- prev^mLs _s.eciiorLa_» .i.t_hjas..J).eejn__Er_aJCfi-d_ that_ t ., for 
a given incompressible flow for which two associated functions 
'I'qC q , d) and X(q,d) are defined, there exist two associated 
functions ’Kq.d) and X(q,d) for the corresponding com- 
pressible flow, depending upon two parameters Y and T 1# 

The question is whether the associated functions i f (q,d) and 
X(q,d) belong to the same flow pattern in the physical plane. 
To answer this question it is necessary to fall back once more 
on the equation of compatibility (24); since when ty( a , d ) is 
given, cp(q,.$) is known by solving equations (20) and ( 21 ) . 
Hence, if X(q,d), satisfying equation (23) and approaching 
X Q as C 0 » — => is to be associated with "■l f (q,‘$) for the 
same flow, then it is necessary that the equation of corapat la- 
bility be satisfied. Sxoept in the case of logarithmic singu- 
larity in section 10 where the complete 'l 7 ( q , ^ ) function was 
not discussed, this condition has been properly considered. 

Onoe the relationship between \|/(q f d) and X£;l,£) is 
established, the next object is to cal"cdlEtre-- tirg "flow pattern 
*(*,y) = constant in the physical plane corresponding to 
"'Kq.d) and X(q,d). In the first place, the fact that the 
transformation defined by equations (9) and (10) is generally 
one— to— one must be recalled. Suppose that in the hodograph 
plane there is a line defined by 


^ ( q , d ) = constant = K ( 139 ) 


which will correspond to a definite streamline in the physical 
plane or a definite part of it. The streamline can be obtained 
by eliminating one of the two variables in x(q,9) and y(q,B), 
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To do this, first the equation (139) is solved for dj 
namely, 


d - d q ,K ) 


(140) 


provided ^(q,d) 0, Introducing this relation into 

equation (13) which, when transformed Into polar coordinates, 
are 


x 


y 


cos d 


sin $ 


6X s in# &X 
5q q d# 

d X cos # h X 

dq q d # 


(141) 


( 143) 


gives a parametric representation of this particular stream- 
line corresponding to ty(q,9) * K in the holograph plane. 
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PART III 


IMPROVEMENT OP THE CONVERGENCE OP SOLUTION 
BY THE ASYMPTOTIC PROPERTIES 
OP HYPERGEOMETRIC FUNCTIONS 
12. General Concepts 


The significance of the general solutions constructed 
in part II of the present report when viewed from the prac- 
tical point, rests in the fact that they constitute an ex- 
istence theorem. It shows that an irrotational isentropic 
flow about a "body can he obtained from the corresponding 
problem of an incompressible fluid, if the free-stream Mach 
number is not too high. However, the solution in the form 
of a slowly convergent infinite series cannot be conveniently 
used to obtain numerical values, as the labor of computation 
would be prohibitive. 


By examining the infinite series obtained in part II, 
the essential difference between the compressible flow solu- 
tion and the incompressible flow solution is seen to be as- 
sociated with the fact that, while in incompressible flow 
solution the individual terms of the series are of the forms 




COS V& 

sin 


, -V 


cos 

sin v# 


in compressible flow solution the individual terms of the 
> series are of the forms 


„ cos . cos 

I„(t) G u (t) 

sin u# sin 


If It were possible to write 


<L U J,(t) - [«(<,)] , *- v 9„(T). = [tu'q)] 


-15 


there would be no difference between the incompressible flow 
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solution and the compressible flow solution except the "dis- 
tortion of the speed" q "by the new scale Q. In fact, 
this possibility is realized under the special condition of 
Y * -1 as shown by von Karman (reference l) and Tsien (ref- 
erence 9), 

for the cas e of isentropic flow with the general expo- 
nent Y there is no such scale factor However, if v 

is assumed to be very large, then there 1 b such a function 
<%, at least to a first approximation. In other words, the 
leading term in the asymptotic representations of E^T) 
and G. V ( T ) does give the desired form. On the other hand, 
the use of asymptotic representation necessarily implies an 
approximation. But this defect is not difficult to remedy 
as the difference between an exact hypergeometric function 
and its asymptotic form can be added to. the approximate so- 
lution as a correction term. Sinoe there are an infinite 
number of terms in the series form' of the solution and each 
gives a correction term, the correction terms also constitute 
an infinite series. Therefore, the original infinite series 
is now transformed into a closed function plus another in- 
finite series of correction terms. At first' sight, such a 
transformation seems unable to give a result that will avoid 
the difficulty of prohibitive computational work, But actu- 
ally, owing to the good approximation given by the asymptotic 
representation even for moderate values of v, the correc- 
tion series converges very rapidly. A few terms seem to be 
all that are necessary. Thus, for all practical purposes, 
the original infinite series is now converted into a closed 
function with "speed distortion" plus a few correction terms. 
The fundamentally interesting point is that for the case of 
a general exponent Y, the simple method of speed distor- 
tion will not give an accurate enough solution, (Cf, ref- 
erence 8.) 

The change in type of the differential equation at the 
sonic speed also introduces a singularity in the speed dis- 
tortion function Q. However, by using the correction terms, 
the effeot of the singularity can be limited to a very nar- 
row range in the neighborhood of sonic speed, and! no practi- 
cal difficulty is experienced. This will be made clear by 
the numerical example given in part V of this report. 


13. Asymptotic Solutions of the Hypergeometric Equations 

Let U v (t) and V v (t) be two new dependent variables 
defined by 
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• \ A w . r . ' ■ ■ - . ,f •• 

•■Mr v (T)- «‘V-‘* a - (1 : tf U V (:T.)--‘S -. - ' •• ••'(•143 1 ) 


Xy(T) = r 


‘ . ~ » 

p+i 

3 ' (1 - T) T V (T) 


(144). 


The differential equations (27) and (28) reduce respective- - 
ly to 


U^'( T ) - cp ( T ) + p.g ( T ) . U u (t)., s ; i Q 


(145) 


V (t). 


p 2 q>( T ) + p (t) 
-P‘ 


V„(T) . 0 


.(146') 


where 


qp(r) = 


1 - (3P * l)T 

4 T 2 (l. - T) ' 


p (t) = Mpt ± 3 ) - r t - T) 3 . 
‘ ' 4T 3 (i « T ) 2 


r * 


Both equations (145) and (146) involve a constant param- 
eter v which is real and pogiti.ve but otherwise arbitrary 
for any fijoed constant ' 6. In the interval 0 < T < 1- in .. 
which the flow takes place, the functions cp(T”) and P±p(' T ) 

are finite and continuous except at the extremities T = 0 
and T ss 1 . To avoid the repetition, let equations (145). 
and (146) be replaced by 


u cc,p( T ) 


u 


<¥>(t) 


a 


(t) 


^a,v (T ) 


•(147) 


where U p ^ y ( T ) c U v (t) when a = 6; and U-p , p ( T ) = ,V ^ (t ) 

when a = In the interval 6 ^ T ^ - - 8 , 5 > 0, 

2p + 1 

cp(t) is bounded from zero and is positive. 3*. Horn 
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(reference 18) shoved that when p is a large positive num- 
ber, a pair of solutions of the following forms exist in the 


interval ^concerned* 

, 


u a^ T ) V * * 

fn( T ) A 

. V 

f 13 ( T ) j, 

"U 2 — ' ’ 

U (a) ( T ) w e -VK | 
U a,u v ' •- e • 

' _A 
cp 4 + 

+ 

u ' 

f ^ 3 ^ T ) + 

. Vfi * * 

' " J 

where 



‘ . “»* 


K(r) 



•(t) 4t, 


+ "£iiilAl (i48) 

US J 


lSJ. 


(t) 


u 1 




(149) 


0 < T < 


2P + 1 


(160) 


A constant in equation (150) was left out, as it can be ab- 
sorbed in the constant fact or- in -equations (148) and (149). 
This representation can be shown to be unique as long as v 
remains g^e^ter than a large positive number' N. By substl- 

u a,u( T ) and Uqui?( t ) in equation (147) and choosin 


tuting 
the coefficients 


r , s 


ing 

(t) (r s» 1 and 2; and s*l, 2, 3, 


to make the individual- . t'erins Vanish, equation (147) reduces 

to :‘:.- 


2K 1 


2K 


- t “■ r I 

£> ■ +• 2E 

,, 8 + 1 .. 


1,8+1 


' 3 , 8+1 


1 , 8+1 


"Pa f 3 


x 

. 4 , 


•Vf" 


8 . 




8 — . o , i , 2 , , r 


where f i •■•‘o (-CT ) -= •f s i j o^ T ) * <P‘ '* • ■ The . coefficients f 


r , s 


+ K« f « P a f ltS - t[ fB (151) 

(152) 
(T) : 

then are given successively by’ a first order : ordinary dif--; •• 
ferential equation and their determination does not involve’ 
any difficulty , The problem is thus formally solved,' 

Obviously, the solution' i's- of the exponential ty.pe when 
cp(t) is positive in the range concerned and of an •oscilla- 
tory' type ^hen ' tp(T ) is negative. Now in the inteival 

6 •£ T £ 1 - S , 6>0 where qp ( T ) ^ O'"- when ; T.^. — 4- — , Tpoth. 

* ■ • . ■ „ 2(3 t'l 

types of solution exist, It is "evident that in , the neighbor- 


hood of T 


a" change ' of- charaoter of .the solutions 


2(3+1 
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must take place, "but the manner in which, the transition oc-- 
curs cannot he deduced from equations (148 ) and (149) be- 
cause of the failure of the representation of .the solutions 


in the neighborhood ■ T 
the Stokes phenomenon. 


— . This is closely related to 

2{s + 1 


The method was extended by Jeffreys (reference 19) to 
include the case where cp(T) has a simple root in an inter- 
val under consideration and can be applied suitably to the 
first order of approximation. The general problem has been 
treated by Langer (reference .20-). in. -a series- of papers, .con- 
sidering both the case .where P and T are real and that 
where u and t are complex. Attention was given especially 
to the Stokes phenomenon, and a law of connection of the so- 
lution valid on each side of the critical point was explicitly 
stated. In the present case, hov^ever, only the. first approx- 
imation is used and Jeffreys' method is adopted for conveni- 
ence . 


It is seen from equations (148) and (149) that the 
first approximation depends only on <p(t), and the effect of 
PqA t ) is felt only by the higher order terms. Hence, for 
the first approximation only, equation (147 ) can be written 
as 

u”(t) - u s cp(-r) U v (t) * 0 (153) 


where Up^p = U_p 9 p - Up* Thus, when u > N, the dominant 
terms of the asymptotic solutions are 




(2 ) . 

u„ (t) 



VK 

e 



1 + 
1 ' + 



(154) 

(155) 


Here 0 in each case, denotes the fact that the term 

is uniformly of the order p -1 when p is sufficiently 

large in an interval 6 ^ T ^ — 1 - 5 , 6 > 0 and is a 

28+ 1 

function of p -1 . 
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On the other hand, in the interval 


+ 6 . 


26+1 

where ’cp(T') < 0 and 2 is a pure imaginary quantity ,iuj 
where tu is real, the dominant terms of the asymptotic so- 
lutions must he a linear combination of equations (148) and 
(149) ppd must he of t he forms! 


U. 


(l) 


V 


(t) «* fLi. oos C vcu + e„ ) 

A ' v 

D ! a!(T) ~ 4 sin (uu,+ £ i>) ! 


< T < 1 


(156) 
* \ 

(157) 


where Cj, c 3 , 


and 


■v 


are constants to he determined.' 


The question of determination of these constants is ac- 
tually the same as that of determining the mode of continua- 
tion of the asymptotic representation of the solutions in the 

range. ^ - - + {j^T^l-§. This can he done, according 

to Jeffreys, by considering the solutions valid in the neigh- 


Let I - T - 


When I is suf- 


horhood of T = _ 

36+1 26+1 

flciently small and U is large, equation (l£3) can he written 
approximately as 


U p "(t)+ V 3 (p‘(0)-| U V U) = 0 


provided 


cp 


<n) 


( 0 ) 


n l cp 1 ( 0 ) 

The • independent integrals are 

, £* H<°U) . £* H^U) ; 


(158) 

1. This is known as Stokes equation. 


with 


t - ~ Ixp'^(O) P (159) 


where' 


are the Hankel functions of 


and Hi £ £ ) 

TT '3' 

order ■ ± -,' Consider as two independent solutions the follow- 
3 

ing linear combinations : 

(i) 


U y U) 


£ H^(£) + I® H J) 

z 


(160) 
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5? 


' ■ : ' • ^ s) ce) =e^H| l} (t) ■ ■ (161) 

■ ‘ .3 • S 

As ^(£) and H^ a ^(£) are analytic Emotions in the 

whole" 5 ^ t— plane, thif immediately provides a means of iden- 
tifying the asymptotic forms that’ represent the same func- 
tion. 


Suppose first that for arg 1=0, the solutions are 
given in equations (ISO.), and. (l6E_L The_n ama solut i ons f op 


which 

arg t 

= rr 

and' 

. • ’3 

arg £ - — rr 

are 



D ' l)(£) 


TTi 

e 3 

(O 

H 1 • 1 
3 

ft 3Tri \ A 
e -a J+ i 

El (s ) / 
e a H t ( 

3 

f STTi \ 
y e 2 J 

(162) 

^ 3) (e) 

-e 

nl 
e 2 

3 

f. »?i\ A 

J. .,»)-£ 

EL (a ) / 

3 

f btt 1\ 
t e 3 j 

(163) 


Now ■ 


H 


U >./, ^ f.i\ ..U) 


r (t 


e * e 


) 


e 8 H 


I*'- (t 


nL\ 
e 3 ) 


H i 2) ( £ ®" i ) “ 3 008 | H i a) ( £ ’ e ^) * ^ ( £ e ^) 


ni 


and when t' is large and -tt < arg t e 3 < it, the dominant 


terms of the asymptotic expansions of (£ e 3 ’) 

(a) TTi 


TTi 


H 


T ( £ -*)' 


are 




and 



58 


NAQA TN- No. 995 


By substituting in equations (162) and (163) and neglecting 
the term of lower ord'er in £, there is obtained by expand- 
ing at the same time equations (160) and (161): 


2 1 


A 

*4 


• 0 - 5) 

r 4 Sin (t - 1) 


-t 

I e 

£ 

£ e 


(164) 

T 

(165) 


Here the arrow is used to indicate the transition of the as- 
ymptotic representation of the same function from the left- 

_A _A 

hand to the right-hand member. For small £, | 4 ~ cp ‘ 1 > 

and $ ~ — UUJ; . (156) and (157) finally become 


n. C .°(T) 


3 


V 


(a) 


q> 


( V 


COS ( UUJ - 


TJ_ 

4 


1+0 


ft)} 


V (T) ~ 006 ft" + f) ft + 0 ft)j 




< T < 1 


(166) 


(167) 


with Cj = 2, c 3 = -1, and e-y = * Under the hypothesis 

just made, the pair of expressions (164), (166) and (155), 

(167) actually represent respectively the dominant terms 
of the two asymptotic expansions of the solutions U J a )( t ) 
and U^ ; (t) for a v which may be any positive but large 
numb e r , 


^ ( a p « i 


V 


14, The Asymptotic Representation of 
and 

The dominant terms of the asymptotic expansion, of 


F(a v + p, by + Pi c v ; T ) 


T> 


Uy X ^(T) and U^ 3 ^(t) are given respectively by (154), 
(166) and (155), (167), By evaluating the simple integrals 

in (l54) and ( 1 6 6 ) , the explicit expressions for the first 


approximation of 




and 


u< s) (r) 


are 
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u^Ct) - *”*(T) 


u(cc-i ) ^ , '.-s-4 * '■" <-u+a 

U'“’ : Ct). r (at)—*— filizll ) 4 T~^ t- v *(t) 

u . , .• . U.- a 3 TJ 


(3) 


0< T < i 

2 6+1 


(168) 


(169) 


‘ ‘ ' * . , ■ * 

\ r ) ~ 2 f 4 ^ 1 ~ T ) - ) ^ cos (You - j"V ' " . , (17 0) 

p ,, i v a s T - l j \ 4 / ~ • 4 . . \ 


. J 


• 4 i!> (~)-~ cos (’ xd * 4) 


36+1 


-< T< 1 


. (171) 


where'-. 


T * ( T‘) = ^ 


a( 1 - t )' 2 ; + (1 - co 2 T 


a 


(1 - T ) ^ + (l.~ a s T)^ 


a 


'V+JL 

7-1 




(172) 


a) 


(t-) = a ten" 


a -T.'— j 
a. 3 ( 1 - T ) 


- tan. 


-i / a _!..r-l 
y 1 - t 


The values of ; the function u>(.T.5 *' are' -given, in figure 3 to- 
gether with the •' funclj lo.n |a(t), 'defined- "by cos p. = l/M. 

In the respective ranges ' of 'validity, each pair of expres- 
sions differs from the exact solution', only by a constant 
factor which can he determined to satisfy, the normalization 
conditions ( 36 ) and (36). , ; By substituting 'equation (168) 
into equation (143), these'^were found to be 


' ' -■ 

°-=75< se) e 


r 

{(l+ 2 a)®} 


Thus, the expressions for the desired asymptotic forms, when 
1 ? > Nj are, for the interval 0 S T < ^ , 
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T (t) ~ f ( t) T V (t) 


— V 


0 (t) ~ f ( t) T (t) 


(173) 

(174) 


where 


. £+i. 


‘ [a(l- t)^ + (1 - a S T)^] 


a 


f (t)= ( X -- T . t(t) a § : T r- 

. (l -r a 3 1) (l + a) a '• (1 -"t ) 2+ (i „ a ®T) * 


(175) 


i. 


’Tor the interval 


20+ 1 


< T < l, they are 


where 


F.(.T).~ f'( T ) T ^ ( T ) , c 0 s (Vo - 

0 (t) - f ( T ) T~ V ('T-). oos (vuj + 

1 ^ , • , ■. ■ i a • ’ ' f 


f(r .to-* 


(a T — 1 ) 


ui-df ys^T 


(176) 

(177) 

(178) 


i-he values' of t(t) are given (fig, 4) as a function of T 
together with the local Mach number M, 

Similarly, as from (153) U v ( f ) ,~ .V y '( t) , corresponding 
expressions, for + 0-, 'P | • c v » t) are: 


£ v (t) ~ -g(.T) : T U (t)- 


• .(179) 


0 <£ -T < 


-i. 


where 


& (t) - g (t ) t" U (t)’ - • > ; 2P+ 1 1 (180) 

“ v • ... . I 

: ' -iVi- • *•. 

<(T) - f 

(1 - a 3 T)t 


( 181 ) 
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and 


9 (t) ~ g ( t ) T V (t) cos (w> - &■') 
V \ 4 / 


5 (t) 1 g ( T ) T~ u (t) COS ( vu ) + 

u 2 \ 4 / 


1 < T < 1 

36+ 1 


where 


- 1+ £ 

2 4 


g(T) = 3 


(1- T), 
(a S T - 1) 




(182) 

(183) 

(184) 


Here F^(t) denotes invariably the first integral 

y(a Ut b v ; OyiT) while G, y ( T ) , when multiplied by q 2U , 

denotes the second integral S' u (t), defined by equation 

(37) when u is not an integer or by equation (39) when u 
is an integer, since the asymptotic expansions are valid for 
both integral and nonintegral values of u, provided v > N, 


In the domains of validity, the asymptotic expansions 
may be differentiated with respect to T with the same order 
of approximation. Hence, for u > N, it can be shown that 


for 0 <: t < 


36+ 1 


where 

h( t) a 3 (l-r J) 
1 

26+ 1 


I. v ^(t) - h ( t ) T V (t) + 0 J> 
^(t) - h ( t ) T” v (t) |i + 0 ^ 


a 


and for 


(1 

< T < 1 


(185) 

(186) 


- <x S T) 4 £<1 «. T j^+ (l - a 2 r)^j (187) 


^ ( T ) ~ h ( T ) 1 3 ( T ) cos ^VU) - m> - i) I 1 * ° (0 } (i8s > 
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6 ( T ) - i. h ( T ) T“ u ( T ) cos ^uu. > + n + jjL + ° ($) } (is9) 


where 


a c 


-* 


h(T)«4(l-T) (a 3 T-l)” ? (20T) , h(t)=cos 


■ 1/ 1 - T 

4 3PT 


(190) 


The values of the functions g(t) and h(”r) are given in 
figure 5. 


It is interesting to note that when V a t 1 the con- 
stant a vanishes and only the exponential type of solutions 
exist. In the case of ^(t) the solution is exact, namely, 


for 



i. ( T ) = 
u 



(191) 


F (T) a 

~-V 



(192) 


of which the first is in agreement with the result obtained 
by Tslen (reference 9), while for Xy( T ) the solutions which 

are not exact reduce to . • 
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-His may fee the cause that destroys the analogy between the 
coordinates of the corresponding compressible flows and the 
incompressible flows. 


For y 


1.406 and V 



n feeing a positive 


integer, the three groups of functions F ^(t); 

F (t); and IL ,(t), „ „( T ). together with their asymp- 

totic expressions were calculated for T varying from 0 
to 0,5 and n from 0 to 10. The results are presented 
in tables 2 to 13. The behavior of the approximation is il- 
lustrated in figures .6 to 11. It can fee observed that the 
degree of approximation of the functions increases, on the 
one hand, with v for any fixed t. On the other hand, for 
any fixed n, the approximation becomes worse as T ap- 


f 


proaches the critical point T 


26+1' 


corresponding to 


the local sonic speed. On the 

T = — 1 is not reached, the 

213+ 1 


whole, if the critical point 
agreement can generally be 


regarded as excellent, especially for larger values of n. 


15. Transformation of the Function *^(w;t) 

.Branch Point of Order 1 

The funotion tf(w;T) for a flow that possesses a branch 
point of order 1 was given in sections 8 and 9, The forms of 
representation, as can be seen, are not, in general, sxvitable 
for practical calculation. The difficulty is twofoldJ First, 
the series involves infinitely many hypergeometric functions 
which themselves are, in turn, defined as infinite series. 

The convergence of hypergeometric series generally decreases 
with an increase of the parameter This means that it is 

very difficult to compute the value of the later terms of the 
series for W(wjT). Secondly, the convergence of the power 
series defining the function W(w; t) itself is, as expected, 
very slow in the neighborhood of the sirgularity. To in- 
crease the convergence, the following method is used! 

Observe that the corresponding function for the incom- 
pressible flow that has the same character of singularity is 
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W. 


(i) 


w 

(w) = y A n w n , I w | < U 


n=o 


which, is absolutely and uniformly convergent in any closed 

~( I* ) . y 

domain in (w| < U. Now, if in (92) 3? v T ) is replaced 


~n 


*>y 


(t) s t n (t), o < t < — l — 

" 26+1 


f (T x ) 


(195) 


T (t) 

- 1 ( T x 5 

then it is clear that 


where t(r) = as by hypothesis, 0 < T x < < 


W ^ ^ (w; T ) ~ -iSZ 


CO 

$ l 


f(T x ) / n 
n=o 


A (tw) n , [ tw| < U 


(196) 


v;hich is also absolutely and uniformly convergent in the 

same domain as W (w) and, consequently, (196) will be de- 
f ( T ) 

noted by W Q (tw). In doing this, however, the re- 

f(T x ) " 

striction that (195) holds only when ■ n is greater than a 
large number N is violated. The error can be removed by 
adding to and subtracting from (91 ) the quantity given in 
(196); then it follows immediately that 


where 


tf (i) (w;T) = W X ( W ;T) + V^Uvn) 

f ( T ) 

w x Cw;T ) = JJJ-J W 0 (tw) 


(197) 

(198) 


Wa^Cwjt) * ^ A n G n (T) w n , jw[ < 


u 


with 


n=o 


(199) 
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& (t) =y r) (T) - jllL t n (T) 
n ^ f ( T x ) 


Here n is a positive integer. The function 'rf(w; t) then 
i 8 represented by the sum of two functions W 1 (w; T ), which 

(i) 

is of closed form, and W 3 ~ (w;t), which is the difference 

of two convergent power series and hence is also convergent. 
But, according to the theory of asymptotic expansion, G- n (T) 

tends to zero as n approaches infinity. In fact, ^ n ( T ) 

is of order n -1 ; therefore, the convergence of W(w; t) is 
increased "by the order of a -1 . This actually is the gist 
of the whole problem. 


As the form of the representation of the hypergeometric 
function given in equation (195) is valid for all T in 

0 <L T < Wj (v ;T j given by equation (198) holds auto.- 

2 p + 1 

matically even outside the circle |w| = U. For this reason, 
W 1 (wt) should be identical in form with that derived from 
equation (101). That this is the case can be seen from the 
following consideration. For, in addition to equation (195), 
if it is assumod that 


-V , 


0 (t) s f(t) T (t) 


( 200 ) 


it follows that 


1(f ) » - l ( Tl ) . p—S. h . 
V J 1 - T. 


(301) 


thep equations (106) and (lO?) yield, by equations (108) and 

(10^), 


+* Cf 


B* = Jp- T" U (t ), 0* S T U (T Z ) 

n f ( t x ) 1 n 1 


( 202 ) 


By using these sets of approximate coefficients and replacing 
2L and. Or (t) by their respective asymptotic expression, 


v 


the following relation is obtained with the aid of equation 

( 100 ) 
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(30S) 


(204) 


Xi ... 

functions ® (t) and G. (t) used in deriving ¥j(w;t), ar< 
-v V ( 1 ) . . 

approximate. Hence, if both are corrected, G^ (*r; and 

G< a >(r) should be of the forms 

* *g * 

g£ iJ (t) * dB* 3 ^(t) + AZ^(t) 

> (2051 

0< a) (T) =. &0* a^C.T) ♦ jVj T V (t 1 ) 42 v (t) 

J 

where 

dB* * b* - » 1 v (t) “ *(*0 * p (t) 

> ( 2 

AC n a °n r^—r ^(Tj) , dG (t) = G (t) - f (t) T -V (t) 

n f( Tl ) V V 


# * 

Here the differences dB n and d 0 n depend upon the condl- 

dition at infinity for any sets of constants B n and C n , 

while those of dJL ( T ) "and dG, (t) are functions of T 

V * w 

only and, for this reason, can be tabulated once for all. 

It also can be shown that the order of dB^ is at least of 

n- 1 and therefore the convergence of (204) is again increased 
by n~\ 

( % ) 

Consequently, if \j/(q,0) - ^i(q.S) + % (q, 6) where 

the superscript (l) denotes either (i ) or (o), and if 
the coefficients are real,, the stream function for the sub- 
sonic flow is according to (9s) .given by 
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•J'jtl.05 = fffjj * 0 (t4.8). 


0 < T < 


~ 2j3 + 1 


(207) 


T O 


(i) 


(q.,6) = - y A n & n^ T ^ S. 11 sin n0, a < U (208) 


n=o 


and in U < q. < V 


Co) 

V U.G) 


V -1 l ( 

= -A K 


(i>, 


(t) 


cl" 0 ♦ 4 a) (x) 


_p 


n= o 


cos u& (209) 


with 0 restricted "by 0 < G < 2'rr, This result is striking 
in that for T = T l( WU.b) = \t'i( T J,8) as & v (t 1 ) * 0; 
that is, the function ^^(q, 6) represents the correct sin- 
gularity of the exact function. Par away from the singular- 
ity the term U' 3 (q,9) (l = i or o) gradually comes 

into prominence, especially near 7 = — — r» hut the con- 
ap + 1 

mergence there is already so rapid that a small number of 
terms is enough to secure a high accuracy in >|/(q,9)« 

It is interesting to estimate the magnitude of the sec- 
ond part of the stream- functi on. By noting the fact that 
G n ( T 3 / = 0, G v (t 3 ) * 0, the expansions of the G n (T) and 

Op ( T ) are 

& n ( T ) = (t - Tj) + . . ,, 0 < t < t 3 


G v (t) = G v « Ct x ) <t - t 1 ) + . 


T < T < ± 

1 20 + 1 


Then from corollary (52), it is shown that for d =js 0 

*a l> < 1 . 6 ) - (T - T > + , . . 

In other words, the second part of the solution is of the 
order of magnitude of (t — t 3 ). However, the magnitude of 
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(T — T T ) depends essentially upon T x for a given incoa- 
pressible flow. If Tj. is not small, then when T = 0, 

I T - TjJ will he large. Thus for large free-stream Mach 
numbers, the second part of the solution cannot he neg- 

lected, This means that for high free-stream Mach numbers 
the correct solution for compressible flow is considerably 
more complicated than the usually assumed simple speed dis- 
tortion rule would lead ope to believe. Thus, any theory 
based upon such a simple rule cannot be accurate enough for 
transonic flows. 

On the other hand, if t is small, or t < < —It-—, 

1 1 2 fs + 1 

then the value of |T t T x I for T = 0 is small. Further- 
more, if the maximum velocity of the flow is well below the 
sonic velocity, then the maximum value of T also is small, 
thus | T - t x i for the whole flow i6 small, Then the sec- 
ond part of the solution is negligible. However, even 

then the solution for the compressible flow cannot be ex- 
pressed as the solution of the incompressible solution by a 
simple distortion of the velocity scale, as Is generally as- 
sumed in the so-called hodograph method. For this would be 
the case only if the multiplying factor* fCO/fCTj) is 
identically equal to 1, Since the multiplying factor is a 
funotion of the magnitude of velocity, the streamlines of 
the compressible flow and the streamlines of the incompress- 
ible flow cannot be made to correspond to each other. On 
the other hand, equation (30?) shows that if is zero, 

then tyj, is also zero. Thus there is this one streamline, 
the zero streamline, in both flows satisfying the require- 
ment of direct mapping. Since the zero streamline generally 
is chosen to' represent the contour of the body* on the sur- 
face of the body in purely subsonic flows, the velocity of 
the compressible flow can be calculated from the incompress- 
ible flow by a simple '/correction formula." This formula is 
given by equating the velocity q of the incompressible 
fluid to the velocity function tq of the compressible flow. 
Thus 



where the subscript o denotes the quantity for incompress- 
ible flow and T(t) ie given by equation (175), This for- 
mula is the same as that suggested by G. Temple and J. Yarwood 
(reference 11), This coincidence of Temple's theory with the 
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present investigation can "be considered as a further substan- 
tiation of the method. 

for the supersonic regions, ^(t) and G. v (t) in (lOl) 
should he replaced by b 


F (t) S f (t) T u (t) cos fwi) - 1L ') 
*“V \ 4 / 




(310) 


& (t) 

“U 


^ f ( T ) T- U (t) cop (v U) + 
& ' 


I) 


> 


5£ + 1 


< T < 1 


.-( 211 ) 


where f(T), T(t) and u>(t) are given in (l?8) and (l?2); 
then by writing 


F (t) 

v 


s i. 

2 


f (t) {. 


i(pcu~2-) - i ( ULU - ) 


+ e 


} 


and substituting as before in eauation (lOl), it leads again 

(o') 

to the sxxm of W 1 (w;T) and W a (w;t), where 


N 1 (w;t) » if { t^T [® 1 Y, { Bn (tweitU)V + C a (twela)) | 

n=o 

v -1 ( , -iuj s V . -icu 

+ e 4 i ) <B n (twe ) + C n (twe ) J 

n= o 

and 

« S (»!T) . i Y »” * 4 s ’( t) w- u ). 

n= o 

According to equation (lOO) , Wj(w;T) also can be summed! 


( 212 ) 
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Furthermore, from (17£$) it can he seen that | tw j = ^ 

"being a constant given by 

a 

\ = J - — > 1 (313) 

(1 + a) a ( 2^T 1 p T(Tj ) 


which is a function of the Mach number and the characteristic 
constant 6 of the gas but independent of the shape of the 
boundary. The value of this function ^ is given in table 
14 and figure 12 for Y = 1.405. As a consequence, the func- 
tions constituting the stream function for the supersonic 
flow are 



f(t) 

f(T a ) 


\J/ Q (0 + u;) + \j/ Q ( d *• U)) + $ 0 ( d + u>) - ij» 0 (d-u)) 


(214) 


■b - to > 0 

'Mn,*) 


= Y ( G p l)(T) + 4 a) < T ) <T V } cos v*. U < q < 7 
n=o 


(216) 


Here the funotlons and $ 0 are defined, on account of 

(213) , by 

i|f Q (d ± i.u) s \t/ o (^u, b * U)) , $ 0 (b ± w) s cp 0 (M7, d ± to) (216) 


where <p 0 ‘ and \|/ Q are the velocity potential and the stream 
function, respectively, 1 of the corresponding incompressible' 


flow. The functions Gp (t) and G-p (t) are the same as 
defined in (205) except that 'dF (t)‘ and hG.^(T) now are 
given by 
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^(t) = Fy (t) - ^ COB (vu) - j) 

> (217) 

AG (t) = G (t) - T" U cos (vuj + j') 

“u 2 \ 4 / 

Unlike the previous calculations, G^ 1 ^!) in (211) is 
not of the order of p” 1 due to the presence of l/2 in front 

of f(t) T^ co$ (uu) - §•)• This, however, does not offer a 
serious objection, as the series in which it appears already 
converges with (tq) v , t "being less than unity, 

It is worth noting, moreover, that in the hyperbolic 
domain the function ^(q.d) depends, aside from a factor 
f(i), only on the two independent families of characteristics 
defined by 

| = d + uj, n = v - ui (218) 

This result is most striking, as it shows that the main part 
of the solution satisfies the simple wave equation and thus 
clearly demonstrates its hyperbolic character. With both 
the incompressible stream function and the incompress- 

ible potential function cp 0 appearing in the solution, it 
i« impossible to establish a simple relation between the in- 
compressible streamlines and the compressible streamlines. 
Since such a simple relation is the foundation of the so- 
called speed correction formula for a quick estimation of 
velocity distribution in compressible flow from that of in- 
compressible flow over the same body, this idea cannot be 
extended to supersonic regions. On the other hand, this 
also indicates that although the differential equation for 
^(q,d) is hyperbolic in the supersonic range, it cannot be 
reduced to the simple wave equation by a mere distortion of 
the speed scale as given by the function u»(t). For if this 
were the case, then (q, v) would constitute an exact so- 
lution without the additional ^3 ^(q,5). This fact Is all 
the more important as the additional ^2 (q,d) is not 

small in comparison with ^i(q,$) for the mixed subsonic 
and supersonic flows, especially for the transitional region 
near sonic velocity. However, in the case of pure supersonic 

flow, tyi°^(q.,d) might be small; then ^(a, d) alone may 
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"be used as a satisfactory approximation. Of course, when 
V * -1, then, as in the corresponding case in subsonic flow, 
the exact differential equation for the stream function can 
he reduced to the simple wave equation. In this case, the 
appropriate form for the speed function U) is 



(219) 


where the subscript l denotes the conditions at the point 
of tangency of 'the true isentropic curve and the approximat- 
ing tangent. This agrees .with the result obtained by N. 
Coburn . (See reference 21.) 


16. Continuation i Logarithmic Singularity 


In the case of the logarithmic singularity the function 
W(w; t) was broken up into two parts of which only the one 
that characterizes the singularity was’ given in equations 
(115) and (116).; As an example, it is proposed to show that 
thi6 problem can be treated by the same method. If the same 
approximation is introduced as In equations (195) and (20l), 
then the coefficients defined in equations (121) and (122) 
become approximately! 


£ 0 B l €shl - 

n f(-T 1 ) * n n f (t j, ) 


( 220 ) 


with B f (tj ) ~ 2 * s0 chosen that the form of equation 

(207) is again preserved. With these coefficients and if 
there is written for the function inside the circle q * U: 


Jt) 

V 


(q. d) 


$i(q,0) + 


~(i ) . 

Vs (q, 


Equation (115)- reduces to the sum of 


** f ( T ) ~ 

Vi(q.d) = yyy-y V 0 (tq, 


0 < T < 


2B+ 1 


(2 21 ) 
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~ ( i ) 


n 

(a, 3) - Y — ® (t) (’i.') oos n$, q < U 
/ , n u mJ / 


(222) 


where 


with 


„ hF (t) 

<* n <T) = I (T) A& (Ti) + — — =^-5- — - 
n -n —a f ( T x ) (t, ) 


AF (t) = F (t) - f ( T ) T (t) 
n ~ n. 


l 1 9~ y^ ( / b 1 > k "II / k 

Ti = f%ry ■ f Ti) n 


(223) 


(224) 


Similarly, in the case of equation (116) it reduces to 


*~( o ) , . " . "-( 0 ) 

\|/ (q, d) = ^ (q, $) + \\f a (a, 


~( 0 ) 


Here %(a,d) is again the same as (331); while \i/ s (q,d) 
is 


_i 


2f 3 (Tj ) 


» tSt '•* 




cos 


(325) 


where 


with 


G^°Yt) = G (T) AF (T x ) + f" 1 (T 1 ) T n (T 1 ) AG (t) (226) 

It n n x * * r\ 


AF n ( T 1 ) = -f-r-' V ~ T , AG n (T) = G (T) ~ f (T) T“ n (T) (237) 
n f 3 (t j ) f (t j ) n n 
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Unlike the previous case, $(q,d) - ty 0 (q,d) when, and only 
when, c 0 tends to infinity. Because of (221), however, 
the singularity of \j/(q,d) remains unchanged. 

Again, if in (llS) 

G. (T ) = i f (t ) T~ n (t ) cos (ntu + N-'') 
n 2 V 4/ 


is substituted for 0 (T) It can similarly be shown that 


'i'l (q . i} a 3 ^ ^ o (d+Oj) + $ 0 ( -6— u>> - $ 0 (d+tu) + $ 0 (d-(ju) (228) 


^3 0) (q,^) = / (1-T) 11+ 2 

af*(T x ) d 


T 

P (i T ) P djj. f.(r) 

J t 

i 


v - U) > 0 


(log A - U)) 


f(r x ) 

y ~y~- (|) cos nd (229) 


co g(°) -n 


n+ i 


where ty 0 (^ ± 10 ) and ^ 10 ) are defined analogously to 

(216), and AG. (T) in G^°^ (T) is now given by 
n n . 


AG (t) * G (t) - 1 f(T) f' n (T) cos (ntu + *r) (230) 

n n 2 \ 4 / 


This seems to indicate that the results obtained so far 
for %(q, d) are quite general: It. may differ for different 

oases, at most, by a constant factor. The general property, 
however, is not shared by \l/ 3 (q,d), the character of which 
changes radically with the nature of the singularity and the 
shape of the boundary. Hence, its importance in the present 
problem is evident. 
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17. The Coordinate Functions x(q,b) and y(q,b) 

Whenever the function X(q,b) for a boundary problem 
is determined, the coordinate functions x(q,b) and y(q,b) 
can he calculated according to equations (141 ) and (143). 
Suppose, for instance, a boundary is assigned with the prop- 
erty that the function A (w;T ) is truly described by (94) 
and (llO), of which the real part X(q,b), defined within 

i s 


A n F^ ; (T) q n cos nb, q < U (231) 


A n in (94) are again real and are re- 
F^ r) (T) = £ n (T)/F ti (T 1 ). 

As the series is absolutely and uniformly convergent in 
q < U, it can be differentiated partially term by term with 
respect to q and d . When the differential coefficients 
X^ta.b) and X^Ca.b) are calculated and are substituted in 

equations (l4l) and (142), there results J 


the circle ! w I = U, 

CO 

, X V 1 

X ( q , & ) = ) 

n=o 

where the constants 
garded as known, and 


(q,b) = V n A n F^ r ^ q n ‘" 1 cos (n - l) b 


n=i 


CO 

V 


~(r ) 


-Pt ; n Ay, F (t) q 11- cos nb cos b 

A_, n + 1 n * x 

n=i 

q < U 

CP 

y(q,b) = ~ Y n A^ ^ q. n_1 Sin (n - l) b 
n=i 

y n (T) 

z_. n n + l^.i 


(232) 


~PT 


cos nb sin b (233) 


n= i 


where 


^(r) ^ , _ F(a n + £ + 1 , b n + g + 1 ; c n + 1; t) 
~n,i 7 “ ^(a n , b n ; c n ; t 2 ) 


( 234 ) 
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Nov , • s ince 

CD 

' x 0 (q.*<>) = n A n q. n_1 cos (n - l) d 


and 


n=i 


CO 

y 0 ( a -.d) ■ - ^ n <i n ” 1 


sin (n - l) 


n= i 




CO 

A n q 11 sin ni 

n=i 


■by introducing the approximation given by equations (179) 
and (185), that is 



~ n 


(r) £ &( r y. t n 

' ' f(T,) ‘ 


(T) 


0 < T < A 

2 P + 1 


by defining 


t r l < T > 

"T1 . 1 . 


h ( T ) 
f (Tj, ) 


t n (T) 


x(q,d) = x 1 (q,^) + (a, d) 


y(q,d) = yi(q,d) + yg 1 ^ (a.d) 


( 235 ) 
( 336 ) 


it can be shown by the same manner that 

x x ( q,d) = -frr T t (t) *o<t q,d) - ^ ~- T J Q 0 (tq,d) cos d (237) 
f( T l^ q f(T,) 


0 < T < 


2P + 1 


y i (q, d ) = • £ -- - T J- t (t ) y 0 (tq , d) - — iiilJL ft 0 (tq,d ) sin d ( 238 ) 
f(T a ) q f(T x ) 
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and 


00 

(q,d) * ^ n A n G- n (T) q 11 1 cos (n - l) £ 


n=i 


00 


- PT ^ el A n G nfl (T) q 31 ' 1 c o 8 nd cos £ (S39) 

q < U 


n=i 


y< i} 


00 


(q, d ) = - » A n 


5 n ( T ) q n 1 sin (n - l) 


n= l 

OO 


n A n S njl (T) q cos n£ sin £ 


(240) 


n=i 


where 


& r T ) = * ^,.PjlLJ 2 - -aixL t n (r) (241) 

F(a n , b n ; o n ; T x ) f(T 1 ) 


5 (t ) » a_=_l i g l a n 4- .fi' t ~ 1 • . ^ot£~ 1 ' 1 ; °n +1 i ll - t P (T) (242) 

n ' J n+ 1 F(a n , b n ; c n ; T x ) f(T x ) 




(243) 


On the other hand, the expression for X(q,$) valid 
outside the circle of convergence is 


OO 

X(q,a) = Y [ S £ 2Lu ( t) ^ “ °n ^u (T) *- V ] sin ^ (244) 


n=o 


provided the coefficients B* and 0^ in (llO) are real. 
The functions x(q,d) and y(q,d) corresponding to (244) 
can he found similarly. These are! 
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co 

v - ' r 


x(q,£) = y |ubJ1 u ( t ) q P_1 Binfu-D^+uoJo^CT) n y_1 s in( v+1 )$j 
n=o 

CO 

- {” 5i (T) " u ' 1 
n=o 

+ uc n S.v,i ( T ) sin v& cos 0 (245 ) 


U< q < V 


CD 

y(q,$) = q U-1 C0S (y-l)i^ - u>Cn?p(T ) q u ' 1 cos ( u+1 X’j- 


n=o 


03 

-I 

I 

j *** 

n=o 


uC n 

P + 1 

v - 1 


.y-i 


(246) 


Here the constants B n and C n satisfy the relations 
( 109 ) and (110) and can he reduced to 


B C 

~ * n _ _u . . ~ * n 

3* = — T U (tj), C* S 


f( Tl ) 


f ( T i ) 


T ^ ( T x ) (247) 


provided the same approximation is made as in (202). Further- 
more, 


,* 0 (q, if) " U “ 1 


■ Z - t 


' uB n q u 1 sin (p - l)& + i;C n q~ y_1 sin (v + 1W 


n=o 
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y 0 (q, d) = Y juB Q cos (l 5 " 1 ) d-vG n ^ 1 008 ( y+1 H 


n=o 


} 


#w 

and if H v (t) and By 1 ( T ) . for the high- order terms are 
substituted by the asymptotic forms! namely, 


£^(T) = g(T) T V (T), i 


_-U 


(T) S h(T) T (T); 0 < T < 


3P+ 1 


then in like manner (245) and (246) can be transformed and 
can each be represented by the sum of two functions Xi(qi^), 
yi(q.d), and x 2 (q.,$), ys(q,d). where x 1 and y a are 
the same as (257) and (258); while x z and y 3 are! 

Xg°\cL,^)-^ U -jjjy 1 ^ (T)q U X sin( U-l )d + Gy 3 ^ ( T)q U 1 sin(p+l) dj- 
n=o 

op 

- PT U ^Gy X i(T)q 1? *’ 1 + Gy^J (T)q~ U-1 | sin ”d coed (248) 
n = o 


T, < T < 1 

1 2£+ 1 


7s \q,d) = ^ v^Gy 1 Tjq^^cos (v-1 )d - Gy 2 ^ (T)q U X cos (p + l ) dj> 

n= o 

y^^:i( T )^-l:l( T )^} 8 - ^ sin« (249) 

1=^3 


- 'pT 

h=b 

n(o-) p(a) 


~( a y "(a) 

where G p and G u>x are defined by! 

5 i l5 <T) = 4B; I„(t) + af^T) 

4 S) (t) = Gy(T) + ^( Tl ) AGy 


/ ( 250 ) 
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1 t ( T ) = AB* iLrJL 5 (T) + 3 - n — T~ U (t ) a£ 

n U+ 1 f(T ) nj,l 


(t ) 


v - 1 £ 


Gy S ?(r) = AC 

U > 1 n U+ 1 v,l f (T a ) 


G + 


n m v 


T (Tj) AG 


V . l 


with 


AP (t) = F (t) - h ( T ) T V (t) 

U , 1 U + 1 ^,1 


v+ 1 G ( T ) - h ( T ) t“ v (t) 


AG. (t) 

v , i u - 1 v , a 


> (351) 


> (252) 


T 

while AB n and ^(t) are defined just the same as those given 
in equation (206), 

Similarly, if the hypergeometric functions involved in 
the high-order terns are substituted by 

f (T)Sg(T) T U oos p (T)Sh(T) T U cos ("uuj - ^ - S.') 

~P \ 4 / *”i> , a \ 4 / 

G (T)Si.g(T) t" U cos (Vo+nA G (T)Slh(T) T ~ V cos ( vu »\ i .+%) 
and by resolving the products of the trigonometric functions 


ins tance , 




cos 

(vh) - = sin 

(u- 1 ) ( d + U)) + ((JO- 3-^) j 

8 in 

£(v- l) ($ - to) • 

■ (»-* 

)] 


cos 

(wu + X sin 

( P + 1 ) ( it + UU ) - ( U) ~^ 

sin 

£ (v + 1 ) ($ - 0)} + 

i 
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a 'brief reduction gives when 


2p + 1 


< T < 1 , 


: x («l.*) = f T X °^ + + X o(^ - ®) cos “ “) 


-[ 


Y 0 (£ + *u) - Y ( 


, ( # - u> ) sin ^ *- u> 


Pt h ( t ) 
4q f ( 


~ ^ ^ Q q ( £ + U) ) + Q 0 (^-.U))JC0S ^ ^ 


- Jj 5 Q (£ + u>) - @ Q (£ - 0))J sin (^ + u) j*cos & (253) 


t (t) g( T) I - . . .1 /rr \ 

Yx (q . =“4 f ^ T ) |^L Y o ^ + “>) + ? 0 (^ - ^ ) J cos (5 - W J 
+ [x 0 ($ + oj) - X 0 (£ -,iu)J sin ^ - cu) j 

“ + ®) + " «>)] 003 + j) 

-[e G (* + ®) " 8 0 ^ " «>)J sin - <■«) j* sin ^ ( 254 ) 

by the fact that qt = XU in the interval under consideration. 
Here 

X l> ± uj) = x (XU, £ ± uj), Y n (£ ± to) = y n (XU, £ ± w) 


0 Q ( £ ±a>) = © 0 (XU, & ± a>) , Q n (il ± to) = Q n ( XU, S- ± a>) 
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where 


®o ( «' 4) - jr 


and 

(o) 


C q , ^ ) 


=» V j^Gy 1 ^ ( T ) q U X sin (l>~ 1 ) d + Gy ^ ( t) q~ V ~ X a in (u + 1) d*l 

n=o i 

-pT y u (T)q v_1 + G^(T)q” u ~ 1 | sin vd cob d (265) 

n=o * i 


( 0 ) . 

y 3 (q.,d) 


28+ 1 


< T < 1 


o° i f f \ f s » 

** y v (f)q U 1 003 (v — 1 ) - Gy (t ) q -u "~ X cos (p+l)dj 

n=o 

-P T V ^ 1 (T)q U '" 1 + 1 (T)q" u “ 1 ^ sin ud sin d 


(256) 


~(a) . ~(a) , . 

where Gy (T) and Gy X (T) retain the definitions given 

in (250) and (25l) except that A? (t) A? (t) AG (t), 

~ . H; u , i v 

and AG^ ^T) are replaced by 

^(T) « ly^ T ) “ g g(T) ^ C0S ( UUJ “ j) 

j) 


V, 1 


£ _ hill rp U 

U + 1 v, l 2 


AjT (T ) = Z — £ - iLU-L rp" C 0S / VOJ - n 


AG (T) = G (T ) - 
V "~V 


S ^ ^ _ l _i_ TT 


T “ COS ( UU> + — ^ 
2 \ 4 / 


>( 257 ) 


ag (t ) = 5 - «r“ u cos fvtu+M.+ n.'N 

u,i u - 1 “v , l 2 \ 4/ 
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respectively. It must "be noted again that the orders of 
3-y 1 (t) and are the same as those of AF (t) 

and AF^ ^(7), respectively, because of the way they are de- 
fined in (257). For the same reason as stated' in section 15, 
this again cannot jeopardize the basic assumption of conver- 
gence of the series. 
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PART 17 


CRITERIA FOR THE UPPER CRITICAL MACH NUMBER 
18. Limiting Line and the Breakdown 
of Isentropic Plow 


The solutions constructed in the previous sections are 
known to he regular in the hodograph plane except at a few 
singular points. It is also known that for the limiting 
case of infinite sonic speed, or c 0 — > aCl , the solution 
will give the desired flow pattern in the physical plane. 

When the sonic speed is finite or when the Mach number of 
the free stream is different from zero, there is no guar- 
antee as to the behavior of the solution in the physical 
plane except the probable continuity of the flow pattern 
with respect to the fr ee-stream Mach number. It is found 
that such continuity in the flow pattern actually exists 
up to a certain Mach number. In other words, the pattern 
of the compressible flow is only slightly different from 
that of the incompressible flow up to a certain Mach number 
at which the so-called limiting lines appear. At the limit- 
ing line, the acceleration of the flow is infinite and the 
flow is reversed. It was shown by Tollmsin (reference 1?) 
and Tsien (reference 2) that, without considering viscosity, 
the flow cannot be continued across the limiting lines, and 
a forbidden region is created in the space where no fluid 
can enter. In other words, continuity of flow pattern exists 
up to a critical Mach number beyond which no isentropic flow 
i 8 possible with the imposed physical boundary conditions. 

The breakdown of isentropic flow, or the compressibility 
burble, can be effected in two ways. First of all, the ac- 
celeration in the neighborhood of the limiting line Is very 
large. Thus each one of the following factors gives appreci- 
able alterations In the dynamic relations? 

(a) Viscous stress due to ordinary viscosity of the 
fluid (reference 22) 

(b) Stress due to expansion or compression of the fluid, 
or viscous stress due to the second viscosity coefficient 
(reference 23, pp . 351 and 358) 
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(c) Small "but appreciable relaxation time required for 
the vibrational modes of the molecules to reach equilibrium 
state (reference 34) 

(d) Heat conduction from fluid element to fluid element 

Secondly, the isentroplc flow also can break down through 
the appearance of shock waves. The breakdown of isentroplc 
flow is associated with the introduction of vorticity to 
the flow. Thus the flow becomes rotational with part' of 
the mechanical energy of the fluid converted into heat 
energy. All these factors tend to increase the entropy of 
the fluid and finally to increase the drag of the body. 

Thus the critical Mach number so defined is of great physical 
importance to the aerodynamic characteristics of the body 
concerned. 


Of course, the isentroplc flow might break dewn due to 
the instability of flow fluid with the final appearance of 
shock waves. Furthermore, the action of boundary layer and 
possible condensation of one component of the fluid 1 on the 
flow might lead also to the premature destruction of the 
lsentropic flow. On the other hand, shock waves can appear 
only in supersonic flow; thus, if the speed of th* fluid is 
everywhere subsonic, there is no danger of the compressibility 
burble. Hence, the free— stream Mach number for the first ap- 
pearance of sonic speed in ths field is called the “lower 
critical Mach number"; while the free— stream Mach number for 
the first appearance of limiting lines is called the "upper 
critical Mach number." (See reference 2.) The latter is 
always higher than the former, due to the fact that limiting 
lines appear only in supersonic flow. The actual critical 
Mach number for the compressibility burble must lie between 
these two limits and depends, among other parameters, upon 
the Reynolds number of the flow. 


19, The Condition for the Limiting Line 


At 

limit ing 
that 


the limiting 
line, it was 


hodograph, or the 
shown (references 


hodograph 
1 > 2 , 12 , 


of the 
and 13) 


d( x,y ) 
d( u,v ) 


c 


Po\ S Ts /I 1\ , a“ 
) \|) - ( -a 2 ) ^ 

PV L q \ c * / « J 


= 0 


( 258) 


1 The phenomenon of condensation shocks due to water 
vapor in the air flow around an airfoil was first brought to 
the attention of the authors by Kate Liepmann, who observed 
them in wind-tunnel experiments. 


< 
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Since the factor “before the term is positive for 

supersonic regions only, c < q. , where p 0, the 
limiting line can appear only when the local speed exceeds 
that of sound. It should he noted that the vanishing of 
the Jacobian is the condition for the failure of the hodo— 
graph method, as the transformation (9) and (10) would no 
longer he one-to-one and continuous. Thus, the appearance 
of the limiting lines is then the physical counterpart of 
the singularity of the transformation. 


As ty(T,$) is known, equation (258) defines two lines 
in the holograph plane: 



( 259 ) 
( 260) 


Geometrically, this expresses the fact that the streamline 
\J/ ( q ,£ ) = constant and a characteristic curve belonging to 
either family has a common tangont (reference l). The 
problem can then be formulated based on this property? 
the necessary and sufficient condition for the existence 
of a limiting line is that there exists a solution between 
the two simultaneous equations 


2T 


fJL »-T -v 
La®T~l_ 




\|f = 0 


or 


2t 




* T + v 0 


\if ■ 0 


(261) 

(262) 

( 263) 
( 264) 


where \|/(T S 6) is ' a • def in i t e branch associated with, the 
largest possible t for a given boundary and a freo-stream 
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Mach number. The zero streamline is chosen, as it generally 
gives the highest velocity ana is the place for the earliest 
appearance of the limiting line. 

Generally, these equations may not possess a solution 
for a known function \l/(T,d) when the parameter is 

assigned. This means that there will "be a system of bounda- 
ries corresponding to a sequence of values of M lf for 
which the limiting line does not occur. The first Mach 
number for which equations ( 36l) and (262) have a solution 
will be defined as the upper critical Mach number and the 
corresponding boundary as the critical boundary. 

The actual solution of the equation is, in general, 
difficult owing to the fact that is, in most cases, 

represented by an infinite series. However, if the stream- 
lines are determined in the hodograph plane for the calcu- 
lation of the shape of the body, a simple graphical test 
of whether there is a point of tangency between the zero 
streamline and the characteristic can be easily made. On 
the other hand, if the form (214) and (315), for instance, 
is used, an approximate analytic solution can be obtained 
without involving- much labor. 


20. The Approximate Determination of the 


Upper Critical Mach Number 


(o) 


As can tfe seen from section 


( t, &) relative to ^x( T ,^) 
a 

recedes from the critical circle 


15, the importance of 
will decrease as T 

T = — - — toward the 
28+1 


supersonic region. Tor the first appearance of the 
limiting line, T is. almost always high, especially when 
the boundary is a slender closed body. Let this be the 


case; then \1> S ( T, d) can be ne-glected in comparison with 

*i( T, *) and a great simplification is possible. The zero 


streamline then can be represented approximately by 


'K T, e) = ^i( T,«) * 0 
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furthermore, a s lmple r-educt ion -shows that the two pairs .of 
equations, (26l), (262) and (263), (264). reduce respee- > 
tively to . .. 

TJ > + ijfj(u) = 0 

® 0 u) + ^ ^ 0 (ti) • 

or . 

>;a> + *;<n - ° (se?) 

• ' ® 0 (tj) - $ 0 (h) = 4> 0 ( 0 + t 0 ( i) ( 268) 

where | and n are the characteristic parameters defined 
in equation (.218). This reduction is made possible ."by the ' 

fact that f( T) never vanishes in the interval — - — < T< 1. 

2P+1 

Whenever the stream function and the potential 

function cp 0 of the incompressible flow are given, the 
functions $ 0 and $ 0 can be easily obtained by substi- 
tuting XU for q according to equation (216). Then, 
since X decreases with an increase in the free— stream 
Mach number M* as shown in table 14 and figure 12, the' 
upper oritical Mach number will be given by the largest 
value of X that gives a solution either of equations 
(265) and (266) or equations (267) and (268). An analytical 
solution can be made, ‘as the functions and $> 0 are 

quite simple, ■ • 

. There is, however, an interesting direct geometrical 
interpretation of these sets of equations in the physical, 
plane of the incompressible flow as shown by figure 13,' 
According to equations (216), the functions and $ 0 

are the stream function and the potential function c p Q 

at the constant value of the speed,- XU, Since X > 1, for 
the body shown in figure 13, the constant speed XU curve 
0 ^ forms a loop symmetrical with respect to the y— axis. 


( 265 ) 
(266) 
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The variables are really the angle of inclination of the 
incompressible velocity vector* Along the constant speed 
curve Cj^ from the point S 2 to P, the angle of inclina- 
tion of the velocity vector is monot onically decreasing. 
Therefore, the parameter of the angle of inclination can be 
replaced by the distances along the curve C^. let equation 

(267) be satisfied at the point S » S g ; then 

$£(S s ) = - ^(S 8 ) (269) 


This means that, at the point S = S 2 » the rate of change 
of the potential function <p o along is equal to the 

negative of the rate of change of the stream function \i/ Q , 

Since potential lines and streamlines in incompressible 
flow form an infinitesimal square mesh, this condition 
requires that the angle between the tangent to the curve 
at S = S 2 be 45° , as shown in figure 13* This is 
easily seen by remembering that from S 2 to P, the value 
of the stream fvinction increases while the value of the 
potential function decreases, because of the indicated flow 
direction. Thus the point S 2 can be easily determined by 
this graphical condition. Squat ion ( 268 ) can then be written 
as 


* 0 (b) - ¥ 0 <S) = $ (S a ) + 4 ( S 3 ) (270) 

u u o 0 


If this condition is satisfied at a point S lf then the 
condition for a limiting line is completely satisfied. A 
similar graphical interpretation for the equations (265) 
and (266) can be worked out for the side of the constant 
speed curve lying to the right of the y— axis , Prom these 
considerations, it is clear that the upper critical Mach 
number is the lowest free— stream Mach number which gives 
a constant speed C^_ containing two points, S x and S Sf 
defined by equations (269) and (270). 
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' ' PAST 7 

APPLICATION - ELLIPTIC CYLINDERS 
21, Preliminary Discussions 


This part of the report is devoted to the application 
of the general method, developed in part III - , to the study 
of the flow of a compressible fluid, around an elliptic 
cylinder. According to sections S and 9, if a solution 
were constructed about the stagnation point, the continua- 
tion of this solution would require that conditions (102) 
and (103) and, hence, (106) and (10?) be satisfied. These 
equations involve two sets of hypergeometric functions with 
parameters m and m + l/2, as well as their derivatives. 

To shorten the lengthy calculation’s, in view of the limit-od 
amount of time available', the following approximate procedure 
was adopted, 

Given the domain £>, the solution valid in the annulus 
region, rather than that about the stagnation point, was 
first constructed. The constants which determine the Laurent 
expansion of the solution, B* and C * , for example, are 
now assigned and, consequently, the set of hyper geometr ic 
functions with Integral parameters is not immediately re- 
quired. The difficulty, however, is the question of whether 
it is possible to continue the solution within the circle- of 
convergence. This continuation may not be possible owing to 
the stringent continuity conditions given by equations (102) 
and (103), and to the requirement that the function must be 
regular within the circle q = U. 

This, however, does not offer a serious objection from 
the practical point of view. In the first place, the summed 
function ( 'I', ( q , #) , for instance) actually holds even within 
the circle of convergence q < T J , and the correction function 
'l'sU.3). is generally small compared with ^i(q,$) due to 
the close asymptotic approximation of the hypergeometr ic 
functions in the elliptic domain. In other words, although 
the solution within the circle of convergence strictly repre- 
sents a different flow, numerically it approximates very 
closely that defined in the annulus region. In the second 
place, since this region q < U is relatively unimportant 
In the case of mixed flow, where T x is very much less than 

— that is, for free— stream Mach number considerably less 


2p + l 
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than unity - the inaccuracy of the solution is limited to a 
small region in the hodograph plane. Furthermore, the most 
interesting phenomena of such a flow, such as the appearance 
of limiting lines, always take place in the annulus region* 
Therefore, this modified procedure, although unsatisfactory 
from the general view point, is an expedient capable of 
yielding an interesting result and furnishing a test of the 
practicability of the proposed solution. 

The situation also may be considered from another angle. 
The procedure used in this section can be derived by replac- 
ing the functions ^(t) and C^Ct) with the approximate 

values given in equation (201) in the expressions for the 
coefficients involved in the solution within the annulus 
region, that is, (106) and (107). Thus the procedure may be 
regarded as an appropriate method of approximation. The 

error introduced is generally negligible if Tj < < ■■■ ■ ■ , 

2 p+1 

This is indicated by the fact that the correction function 
for instance, is very small in comparison with 
1|f 1 (q,^) when q < U. 

Another simplification is made by using the elementary 
integral q~ 2U (t) instead of q“ 2U ^(t) in the con- 
tinued solution, as, in this case, JE_ v (t) is a well- 

defined function. In doing so, the asymptotic behavior of 
the second solution remains unchanged because the first term 
in G^Ct) is always small in comparison with the second. 

If, however, all the required hypergeometr ic functions 
are computed, there is no difficulty in carrying out the exact 
method developed in part III of the report for any accurate 
study of two-dimensional flow. For this reason, the expres- 
sions for the hydrodynamic functions derived for both the 
exact and approximate procedures for the problem at hand are 
given. 


In the numerical example, detailed calculations are made 
for the flow of air about a cylindrical body derived from the 
incompressible flow about an elliptic section with a ratio of 
the minor and major axes equal to 0,6. The calculations were 
carried out for two differ'ent free-stream Mach numbers, 0.6 
and 0.7 . 


22. The Functions z 0 (w), W 0 (w) and A 0 (w) 

The irrotational flow of an incompressible fluid about 
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an elliptic cylinder with the velocity at Infinity pscrall el to 
the major axis is represented by the complex potential 

*„( l 0 ) » t + T ( 271) 

with 

s o = ^ + X ( 272 > 

For convenience in practical calculation, all the physical 
quantities z Q , q, and p, will be normalized consistently 

throughout the present part. The major and minor axes of the 
section are respectively 1 + e® and 1 — e 3 » where 
e < 1; q « 1 at infinity and p = 1 when o = 0. This will 
automatically render the hydrodynamic functions dimensionless 
and the constants U and p 0 will be eliminated from the 
formulas in the succeeding sections. 

By differentiating (27l) with respect to z Q , the di- 
mensionless complex velocity of the flow is 



Thus 

| i _ e 2 w ] 4= 0 ( 273 ) 

This function is two— valued with two branch points at w* 1 

and w = e“ 2 . In order to make z q (w) a s ingl e— valued 

function of w, the expression (273) is supposed to be the 
principal value so that I ar g( 1 — w ) 1 <tt and 1 < [w 1 < e -2 . 

The condition ! e s w 1 <1 must be- satisfied, for w = e 
corresponds to t, = 0, which is another singularity. With 
the principal value so defined, if the negative sign in (273) 
iB taken, then the domain H corresponds to the half plane 

Hi t < 0 and 1 t i > 1 . On the other hand, since the'trans— 
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f or mat ion ( 27 2 ) is on e— t o— -one when t£i > 1 , then the 
domain D, which is Rt < 0 with the region inside 
the section excluded, corresponds uniquely to D. 


Consequently, the inverse mapping function z Q (w) is 


z 





(274) 


which will he s ingle— valued , provided a cut is introduced 
to Join the branch points in such a way that the argument 
of ( 1— w ) is restricted to — tt < arg(l-w) < rr and 
1 6 w | < 1 . On separating into real and imaginary parts, 
it is found that as 0 < d < 2tt 


X Q (q, d) = 




|l( q » ■&) + J( q » *) V 

+ ^I £ (q,<) + J _1 (q,^)| j (275) 

y O ( q » ^ = Jjlj- I(q,d)+J(q,d)^ 


~ e a< r I £ (fl3).+ j‘* 1 (q,^)j' ] (276) 


-H 


with w * q e , where the functions l(q,d), I e (q,d), 
and J(q,d) stand for: 


I(q,d) = 


1— (l+e )qcosd+e c 
1 — 8 q cos d + q s 


1 ,(q t d) = f S } ? s - t _ t c . 5 q 3 


2s s q cos d + S 4 q a 


( 277) 


1 


(278) 
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J(q.*) = 


1-2 Z* q co8 £ + e 4 q 2 
L 1 - 2 q coa & + q 3 J 


,11/2 


(279) 


On the other hand, substituting equation (273) in equation 
( 271)1 the function W Q (z 0 ) is carried over into D; namely, 


W n (w) = - 


• l-£ B W ~^ 
l-w 


l-w 


Ll-€ 2 vd J 


( 280 ) 


Now W (w) * cp 0 (q,$) + i\Jr 0 ( q, •& ) , and similarly 


<p B (q.*) = ~x 

2*a 


Y 0 U»*) = “ 

2a 


Li 


A* . J 


l(q,$) + J(q,*)> + <1 (q,£) + J- l (a,#), 


I' 


A 1 


/ j V M « /j J 

t' i — T 


( 281 ) 


i-l(n,<0 + J(q. «)]■■*- J-I (<.,#)+ J- l (q, 

Ll J l £ 




r 
J J 


( 282 ) 


By integrating z 0 (w), according to ( 87 ), the trans- 
formed potential function A 0 (w), aside from a constant, 
takes the form: 


A„(w) = 2 ( l-w )'i (l-€ 3 w)& 


( 283 ) 


The principal value of this function is again defined 
by restricting the argument of (l-v r ) to -tt < ax g (l-w) < tt 

and iwl < e~ 3 . Within this domain 32, the real and imaginary 
part s are: 

•X 0 (q,*) = ^ [K(q.,d) + L(q,*)3& ( 2 84 ) 

OS & < 2tt 

cr c (q,&) = - 2 's C-K(q,£) ■+ L(q,^)]^ (285) 

as A 0 (w) = X 0 (q,^) - icr 0 ( q ) , where the • f unct i ons K(a,£) 
and L(q,d) are defined by: 



NAGA TN No. 995 


95 


K(q,£) = 1 - (l+e 2 )q cos $ + e 2 q 2 cos 2 £ (286) 

L(q,i>) = [1-2 q cos d- + q 2 ji [l-2e 2 q cos £ +e 4 q 2 ]^ (287) 

23- Expansions of W 0 (w) and A 0 (vr) 


The function W p (w) defined in (2SQ) is single-valued 
and regular everywhere "in twi < 1 and, hence, possesses the 
following expansion: 



n=o 


where the coefficients A n are real and given "by 


A n = (1 + £S ) n>l 


(288) 


( 289 ) 


with 


A o = 




n 

1 'X' r ( n- m+-g' ) T ( m-t-g- ) 3m 

77 Zl r( n- m+ 1 ) r ( as+ 1 ) 
m=o 


Hewever , in the yegign outside iwl < 1 the function 
W Q (w) is douhl ed-valued; and when a cut is put between the 
"branch points w = 1 and w = e ~ a , the principal value is 
discontinuous along the positive axis of reals within the 
annulus region. To obtain the desired expansion, the function 
is written in the following form 


, N i 2 - . ( 1+e 2 )w 

^(w) = 2, i i 

t (l-w' 1 )* (1-.€ 2 wF 


(290) 
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tfow ( 1- e 3 w r 2 is s i ngle- valued and continuous 

within the annulus region; its Laurent expansion is 


(l-w 1 )"2 (l-e 3 w)'"^ 



where 


CO 

8<°> [e 3n w n + *-]. 

n=i 



1 < I wi < £ 



CD — 

1 -"V r (n+m+i)r (m+ t) sm 
17 ^ ^ P ( n+m+1 )V (m+1 ) 


2 


m=o 

Substituting (l-w)** 5 ( 1- e 3 w )“ & from (29l) i» (290) 
expansion for W 0 (w) in the annulus region is 


( 291 ) 


( 292 ) 


the 


W 0 (w) 

CO 

n £ 3n W u + 

0 n w U ], 1 < 

J w t 


n=o 




when the 

constants 

®n 

and the exponent 

V 

fined by: 

B n = 

2 £3 S n ll 

- (l+£ 2 ) Sn 0) 



C n = 

2 si o) - 

(1+e 3 ) S 

> 


V = 

n + -1- 




< £“ 2 ( 293 ) 

are de- 


( 29 14 ) 


Similarly, the transformed potential function A Q (w) 
can he expanded and is: 


A 0 ( v ) = 2 ) A n w n , I w I < 1 

n=o 


when the constants A n are 


A n = s£ i} - (1+e 3 ) S^l + £ 3 ] 

Si = - i (1+e 2 ), A 0 = 1 


and S 


(i) 


n 


( 295 ) 


(296) 


is given in (289) 
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On the other hand, in the annulus region the expansion is 


AqC w) = — 21 V |^B n 


L, 

n=o 


w n j, 1 < | w 


e Sn w n + C w“ n |, 1 < | w | < €“ s ( 297) 

n 


with the constants and defined as 


S n " s n-i - (1+£S) ^ + £ ° S l+l' n - 1 


B ' = 2 e 2 S ( 1 ° ) 


/ a . ( o ) 

(l +£ ) s o 


( 29 8 ) 


C = S ( °^ - (l+c a ) + e 8 s ( !> 

n n n+ i n+ s 


where ^ ( e 8 ) is defined in (292), 


24. The Stream Function 'Kq,$) 


The relationship between the domain D and D is 
thus fully established and the functions corresponding to 
such domains are also given. From the general scheme de- 
veloped in sections 8 and 9 the solutions for the similar 
motion of a compressible fluid can be constructed. First 
®f all, the stream function \|/(q,6) governing the subsonic 
flow is the sum of and ty s ( a , #) . According to 


( 207 ) , ( 208 ) , and ( 209 ) , f or 0 < T < 

cp + J. ' 


'I'l ( a , $ ) = 


1 f(T) 


2^ f( 


~ {[-Ktq 
tO LL 


+ J( tq , & ) 


- [ - I e (tq,*) + J- l (tq,«) 2 | 


( 299 ) 
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where the functions l(tq,$), Ig(tq,$), and J(ta_,d) 
are obtained from I, I £ , and J in (272) to (279) by 

replacing q by tq, t being defined in (195). For 
q <1, the function is 


U), 

( q.r 


00 

= ^ A n & n (T) q 11 sin n i 

n = o 


(300) 


where A n is defined in ( 289) and & n (*r) in (199). For 
q > 1 and in subsonic region the function 'i^ 0 ^ ( q ,$ ) : 


,( o) rf (i), ( an v 

♦ a (q.^) = > I (T ) e ,q 

n^o'~ 

( 3 ) — 1 

+ & (t) q | cos v $ , 0 < * < 2n 

(i) (a) 

where (t ) and (t ) are defined by ( 205) with 

the constants B n and C n defined in (294). 

When the motion becomes supersonic, the continuation 
of ^i(q,^) defined in ( 899) gives 


(301) 


^l( q,d ) = 


- Kx.e ) + j(x,e)]^-j-i € (A,i)+j^ 1 (x' i {)j^ 






* 


23+1 


< T <1 


+ |^-I(\,n)+J(A,n)J 8 ~^-I e (X,q) + J -1 (\,q) 

- [+km> +j (m )] - [-y*.e ) + j“ a (x,e )J 

+ j^i(\,T l )+j(\,n) ~ + ^i € (\,n) + J" l (\,n^ j- (302) 

according to (214). Here 5 and T| are the characteristic 
parameters defined in (218). She upper sign in the last two 
terms corresponds t© ri > 0 while the lower one, to rj < 0. 
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The accompanying function '4^ \q,^) is 


Co) f 

(q.^) = 2_. 


r .i°(o «•* *’ 


n=o 


+ &i S) ( T ) q V 


cos v 4, 


2£ + l 


< T < 1 


(303) 


(l) 


Here the functions 


(t) and 


(O 


G- (t) are defined "by 


(205) in conjunction with (217) in such a way that (303) 
will "be the continuation of (30l). It also should he 

1 

< T < 1 


noticed that the variable is restricted to 
1 


instead of 


< T < T l€ - 4 , 


as 


T x e 


-4 


2(3 + 1 
is generally 


2(3 + 1 

greater than unity, which is impossible for the actual gas. 

It should be remembered that 4'g (q,$) is always 

negligible compared with 4»i(q»'&) within and on the unit 

circle q = l when T, is small in conioarison with — - — : 

1 2(3 + 1* 

^(q,^) can be approximately represented by , 'l' 1 (q.,£) alone 
throughout the interior of the unit circle. As a consequence, 
the calculation can be simplified cons iderably by constructing 
first a solution for the annulus region by using £ ( T ) in- 

stead of &^( t) and making an approximate connection across 

the unit circle. In that event, the stream function will be 
reduced to 


'K q , $ ) * q > 3 ) 


(304) 


when 0 < q < 1; here x ( q ) is again defined in ( 299 ). 

On the other hand, when T, < t < — - — , 

1 2g + l 

(o) 

vl/(q,$) = ^i(q.'S) + ( q * ^ ) 


(305) 
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where the function tyg°^(q»^) which is small on q =» 1 
is given by 


't4° ) (q,«) = Y, [ B n & u (T) 

n=o 

Here the functions &u(t) 


€ Sn q V + 


cos u $ (306) 


and G_y( T ) can "be shown to he 


V T) = 4 r) ‘ T >-~ *”• y * ” (so 1 ?) 


and the coefficients B n and C n are defined in (294). 


The continuation of ^(q,^) is naturally the expres- 
sion given in (302) while that of (306) differs only in the 
definition of G v (t ) and G_ v (t) which are 


O v (T)-'4 r) (T) 


(r) 

G^(T) = l k '(T) 
v — u 


1 f(T ) 

2 f(Ti') 


1 

2 


f(T) 

f(Ti) 


cos 


(,. - ;) 


2P+1 


<T< 1 


(308) 


— v 

t cos 


(~* 5 ) 


25. The Coordinate Functions x(q,0) and y(q,d) 

With the functions z (w) and A (w) defined in 

° o 

sections 22 and 23, the corresponding functions A(w;T) 
and consequently z(w;T) for the motion of a compressible 
fluid can he constructed, : The coordinate functions derived 
from A( w ; t ) are given, respectively hy the sum of two 
functions *i(qi6) and yi(q»$) which, according to equa~ 
tions (237) to (238), 


ar e 
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*i( q , ,*) • -sr~— ^ — ”■ y^^Ctq,*) + J(tq.,.*)j' 


+ e s [l £ ( tq , d) + J 1 ( tq. 


(3T h( T ) t sin 24 


) t sin 24 f 

— J — 1 + 4 

2 f(T x ) CT q ( tq ,* ) \ . . 

+ J(tq,«) + e S J"" 1 (tq,*) y 

yi(q,*) - •' ^ T' - r I — *(tq,*) + J( tq , * ) j 
2h f(T x ) \L J 


€ t q cos * — 6 


& 


— £ 


2^ f( 

X c ( tq.,-0 ) + J -1 (tq,*) ^ 

h( T) t sin 3 * 


Pt 


{- 


1+ 4 6 t q cos * — t 


f(T x ) O o (tq,«) 

+ J(tq,*) + € J (tq,*)^" 


where cr (tq,*) is obtained from a (q,*) in ( 285) 
replacing q by tq. The functions x s (q,*) and 
y g (q,*), according to equations ( 239 ) and ( 240 ) , 

00 

x 8 ( q, *) * 2 n A n ^n^ T ^ 1 cos ^ n “^ ^ 

n=i 

CO ... 

~ T n A^ & n t i( T ) q 11 1 cos n * cos * 

n= 1 


a 

(309 ) 


(310) 

^y 

ar e 

(311) 


q < 1 
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y U> 
y a 


(q,d) = - 2 y n A n G n (T) q 11 ** 1 sin (a - 1 ) 


n=i 


CD 

- 28t y n A n 5u. f i(T) q. n ~ x cos n d sin $ (312) 

n-i 

Here the functions O n (T) and G n|1 (T) are defined by equa- 
tions (241) and (242) and the constants A n by (296), 

The same funotions valid in the annulus region are 

again represented by the sums x^q,#) + x a ^°^(q,^) and 

yi(q.»£) + yi°^(q.*^)» Where Xj.(q,S> and yi(q,$) are 
defined by equations (309) and (310), respectively. When 

Tj < T < — 1 — , x|°^(q,d), and ya°^(q,$) are 

26+1 

c 

*3°^ (q,^ ) = - 2 ^ 


U ^^(T) € an q^" 1 sin (v-1) 0 
+ Si a) (T)q“ v “ 1 Bin (v+l)*]+28T Y u (T)e an q U '’ 1 


n-o 


L-, 

n=o 


p ^ \ “P-i 

a u,i (T) * 


sin v « oos 


(313) 


yi°^(q,*) — 2 y u [g^ 1 yT) e an q v - x oos (v-l)d 

i 1 L, 

n=o 


+ Op ^( T ) q u "* 1 cos £v + l) 0 

' A *) 


26 T y p X (T) e an q u ~ l + &p* i( T )q"" U ~ 1 sin sin $ (314) 


n=o 


The functions G^ a ^(T), 0p®i( T ) are defined in equations 

(250) and (251) together with equations (352) with the con- 
stants i3 n and defined in equations (298). 

On the other hand, when — 1 — < t < 1, the continued 

2 8+1 

expressions of x 1 (q,^), yi(q»^) across the critical circle 

T = — — are, acoording to equations (263) and (264), 

28+1 . 
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=i( o.» 


t(T) s ~{[i(x,n + j(x,e il 

Tj.) II 


i* 7 * f( 


% r Tl 

+ e 2 ' f i € ( x , | ) + J“ 1 (x,!)j + [ l(A,n) + j(x,n)j 


h 


+ € 


1 

[l € (X,ri) + j“ 1 (X,n) j* COS 


l 

- 4U K^.e> + j(x,e)l - c s \ ~i e -(x,e) 

2 5 8 f(Ti) LL ' . j l 


+ J 


' 1 ( X , |)1 - -I(x,n) + J ( x , n )] "* + gS [“ 


-i, C\*\ / rr \ P T h(T) cos 4 ffX sin * 

J ( x >J } • Ci -) - 4 , f(V7) - U^ttt 

1 + 4 € 2 X cos | — £ S + J(X,f)+e J (X,£)^ 

— — 3 -- n — ^ ( — 1+4 c 2 X ccs T)~c S + J(X,ri) 
o A X , n ) \ 


+ e 8 J” : 




in £ 

IT 


-C 1 - 


4 € X cos | 


+ £ 2 +j(X.f)+£ J X ( X,|)^ 


— x y - in . ( 1 — 4 £ S X cos n + C S + J ( X | T) ) 

x 0 ( x ,n ) V 

8 J -1 (X,ti)^ s in j" (315) 
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/ .N _ t(T) g(T) f ' 
y iU»*) = ~T7b — r 1 “ 

n5 B '( T,) v. L 


I(A,t) + J(X,t) 


2 “ f( T x ) 

x 1* r 

J“ (x, I) + l-.KA.r,) + J(\,n) 

£ 


- e S [-l e (xJ) 


2 


~ € 


[- 


i e U,n) 


+ J'^U.rO 


+ J(x,|) 

i(x,n) + J(x f n) 


\ ( t \ \ . t(T) g(T ) r [ T/% M 

J C0S 7 ^ 7 oTj‘t“L I(x, ° 

i ir 


— 6 


I,(X,?) + J 


a r 


+ € | I e ( X ,T) ) + J 


■‘a. t>] 


„ _ . , - , PT h(T) sin^/Tx sin i 

x 8ln ^ - V - 4 1 f(T- 1 )-U^TT) 

X X + 46 ^ cos E ~ ( + tX(\,E) +£ J ( X , £ 

i .'AP-? /-! * 4 5 e X CO. T, - e 3 + J(x.'n)') , 


€ S J” 1 (X,tl)) 

' tt\ f 

cos + y - [ 

X s in | 
.Xo ( X , i ) 

e + J ( x , £ ) + 

C 3 J-^X.e)) - 

X gin ri 

X 0 (X ,T)) 

£ S + J( X ,T) ) + 

e s J _ 1 (X t Ti)^ 

s in ■ 


COB 


4 X COS T( 


( 316 ) 
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While x a ( q , ■&) and y^q,#) remain to he defined "by equa- 
tions ( 313 ) an & ( 31 ^) except the functio ns ^ an & 

* \ ^ 

&y^( T ) are replaced "by those given in equations (? 50 ). (?5-0 

together with equations ( 257 ). 

By the same argument as that used for the stream function, 
the practical calculation of x(q,q) and y(q,$) can he sim- 


plified hy 
namely, 

neglecting x^^Cq,#) and y£ 1 ^ ( q , i5- ) when 

q < 1; 


x(q,$) = XiCq,^) 

(317) 


y(q,^) = yxCq,^), 0 < 0 < 1 

(31S) 


where x^q.^) and yiCq, 1 *) are defined in equations (309) 
and (3IO) ; and in the. anauA.U.s...r.egi on .. 


x(q, b) = x x (q, £) 

+ x i°^ U. 1 *) 

( 319 ) 


-1 

H 

A 

A 

H 


y(q, *) = y 1 (q,^) 

+ yi o) (q.^) 

( 320 ) 


Here x 1 (q,i9-) and y 1 (q,^) are either given by equations 
(309)* (310) or (315)* (316). The terms x^°^(q,&) and 

ya°^(q f £), on the other hand, become 




5 /T) = F W '(t) - lili- t - u ( 323 ) 

V — u ry f t .. \ 


p, m t<t 

CD O M'"'* 

w> P O 

H. w 

P 

CD -t t° 

pj (-■ *• 

*> 

ED jA ^ 

at 


-i x I II 



*+■ Vi i- 4 Pi 

J3- Vi 

a h o 


o -l 

Hj ,0 01 

P I m 

P Vi W 

O I Vi P 



o Vi 

P at , h 1 I 

_a> h- h 


P r}. 

I I O 

Pi vi o 

vi <3» + at 

* s 

-4 (0 CO ^ 

H- Ti> ^ 


CJ° », . 

Pi H ^ 

vi o L J 



-1 WJ 

— P 

Pi 

vi 

P 

P H* 

a> w 

CO H 

CO ^ 

•— - m 

w 

p 


M 

M 

w ^ 

O 

o 

<5> 

o> 


X 

+ 

n 

(O 

O } 

1 

vi 

p 

w 

1 

p l 


f* 

+ 

1 

vi 



-1 


Oi 

' — ^ 


p 

lO 

i — i 

PI 

1 

bdi 

vi 

vi 

p 

• 

1 

pi 

M 

w 

vi 



S — s 

’ -l 

CO 

-» 

' — t 

H- 

' ■■ 


P 

m 

tQ 


M 

1 


P 

vi 

Vi 


1 



M 

+ 

fO 

’ 1 


vi 


M 

1 

CO 


w 


<&> 


0 

1 1 

CO 



H- 


+ 

P 

«** 


■ 


CO 


o 


Vi 

o 

Tn 


OT 

-1 

l 

3* 

t? 1 ^ ! 
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CONCLUS IONS 


As an example, the motion of air- past a ayl indr ieal body 

was considered "by taking e =. The flow patterns in the 

T,^- plane for two free— stream Mach numbers M x = 0,6 and 0.7 
have been calculated and were given in figures 14 and 15, It 
should be noticed that there is considerable distortion in 
the shape of the bodies in the compressible flow from that in 
the incompressible flow, If the compressible flow around a 
given bodj»- is desired, a series of computations should be made 
with various geometric parameters 6, so that the desired 
body shape at a definite Mach number M x could be picked out. 

These computations definitely demonstrate the practica- 
bility of the proposed method. They also show that, in the 
case of two-dimensional motion of a compressible fluid, the 
mixed subsonic and supersonic flows exist within the field of 
an irrotational isentropic flow about a suitable body, and the 
transition from one to the other is continuous and reversible. 
Furthermore, the breakdown of the irrotational isentropic flow 
depends soltlyupon the occurrence of limiting lines which, in 
turn, is determined by the condition at infinity or the- shape 
of the boundary,, while the magnitude of the local speed at- 
tained is immaterial. In the case of M x = 0,6, the irrota- 
tional supersonic flow continues to exist up to the local Mach 
number M = 1.25; whereas for M x = 0.7 it breaks down- as 
soon as M =» 1,22 is reached. The singular behavior of the 
streamline is marked by the point of tangency of ^ = 0 with 
a characteristic at M = 1.22, 

The calculation of the flow pattern in the physical plane 
is yet to be completed. When this is done', the pressure dis- 
tribution can be compared with that over the same body of the 
incompressible flow. 


Guggenheim Aeronautical laboratory, 

California Institute of Technology, 

Pasadena, Calif., April 17, 1945, 
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APPENDIX A 


PROOF OF THE OREM (52) 


To facilitate the discussion, equation (7l) is first 
written in the. form: 

® C !i) (t) + a u i s . o 


where 


^ (t) a + ^I'T) + V r) 


e (t) « e v (l) ( T ) + 


JBl. 


v(l - t) 


- Y v (t) 


and 


Vp(T) xx f - 1 - ( 2S . . * 1 l T + 


P®T 3 


, 4(i , T) 8 ^ T 


1 - T 


V 


(1 - t)‘ 


t u hS(t)J 


when u is large, the character of the functions £, and £ a 
can he easily studied in the' T , i ' 1 ' -plane (fig, 16) hy neg- 
lecting the third term under the radical sign. This can he 
justified in the following manner; Consider the case visa v> 
is positive and large hut not an integer. In the interval 

0 < T < 1 j X V < T ) << because f^T) - T“ V F -V (T) bjr 

2& + 1 v _u 

equations (35) and (55). Then rs R„(t) - t a T. . There- 

t? 1/ 

fore, T-g Ry(T) » 1 -when u is large. But both F v ( t) 
and Fp(T) are oont inuous, with respect to v ; so the fore- 
going result applies equally to the case of integral v. 

Hence, the third term in the Radical for Y (t) can he neg- 
lected for large u. . . 

Owing to the manner in which Y^ is defined, correspond- 
ing to each v> there is a line T = T 0 > i . such that 

2p + 1 
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' "Yy (t ) ^ 0 when t < t^-.. ' Ad a consequence £ x and. £ 2 are 

real or .pompiei' con jugate according as T ^ T Q . In 

0 < t <i 0l £i = 0 and £ a = 0 will- give two 1-parameter 

families radiating, frpin (-0, -l) and (0,l), respectively, and 
joining t.ogethpr kt . a point where ^^2 = 0. If 0 < T < T 0 , 
the product' £j,£ a may he negative or positive according to 
whether the point lies, to the left 1 or”' the ' right of the curve 
£i a 0 and £ 2 = 0-. , On the other hand, if” T > t 0 , £i £ s 

1 8 always f positive . 


Now 


= £, while the initial, slope of 


is p.^1. ‘ the integral curve, must lie' above ^ = 0, 

and.helow £3 =* 0. 'If'it were not,' the integral curve would 
croa,s the curve « O', ■= 0, where ^(t) = 0, and- 

■ would he negative somewhere in 0 < T < 1 

■(I)-' - - 3P + 1 

This is not possible, for ~ by an argument simi- 

lar to that used for determining the ..magnitude of t" 2 R p (T) 

1 ) T ^ k- . 


and according to (55) T f 1 + 1 

. . • v • 1 — T 

< T 

(l) 


> r > -ci - t)* 

-U 

0 < T < ■ — _i . Eence | 1 ( 1 ) > 0 in 0 < T < — and 


in 


2{3 +' 1 V , ; , ~ - 2P + 1 

% , * 

p continues. to increase until it intersects with ■ = 0. 


(1) 


v 


<■ 0 and never 


After it crosses the curve £1 fe~ 0, | 

changes sign as £ s £ a > 0 -in. ’t 0 < T ■< 1. Consequently, 


(t) is monotonic and decreasing in the interval 
When 1? is sufficiently large, T ( 

1 ' . _ 1 


T q < T <,1. 


proach very rapidly to 
becomes infinite. 2& + 1 


and 


will ap- 
when p 


2g + 1 

Proof of theorem (52) .- Form the following identity: 


(1) 


(O 


xitn;*') ^ (n;*' - + 01; 1 ' - + % << 

+ e B (l) )] * ? t„ a u - T) 2 VV* ( ^ l) W 1>) V 8inh v 

_ £ v (l) ) £ 0 (Al) 


(i) 


(l)xf 


(O 
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It can be shown $hat the differential expression possesses an 
integration factor 



It wj.ll be noticed that the sigh of (A2) is determined by the 
first factor 1 ^ ) • only.. On multiplying (Al) by 

(A2) and integrating the resulting total differential from t 0 

to T, with a suitably chosen initial value 
it is found t.hat 1 : 


- T)-' P * t* S« 0 (T 0 )B= v (T 0 j 


coflh v J < T1 1 5 - 1 1 - 1 

T n . . 


(l) 


> 0 


which is positive 'if and only i f T) J- 1 * - £ $ 0 every- 

where in T 0 < T < 1 . Since both i v ( 1 ) and T) u ( 1 ) are con- 
tinuous and moi^otpnic, the condition is both necessary and 
sufficient. Furthermore, it should be noticed that the condi- 
tion = Iu^ 1 ^(t 0 ) is purely a convenience. If 

Tlv (l) ( T o) ^ the validity of the theorem is not in' 

the least impaired. 1 
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APPENDIX B 

■PROOX OX THEOREM (22) 


Consider the first series: Multiplying throughout the 

ineauality (5&) f namely. 


U'r) > 


' 1-(2P+1)t 
1 - T 


0 < T < 


2P+1 


hy 2 ^ and integrating hoth sides from T to T x shows 
that 

X^(t) < *i(t) 

T x (t) 

where t-^r) = ~ — (~r~) ~ ^hen it follows that 


A n 


A n (t 1 w) 


n 


Now 


cc 1 ■ 

A n (tjw) 11 conTsrges when 1 1 x vr j < U due to 


n=o 

equation (28). By Weirstrass’s theorem, the series (92) is 
uniformly and absolutely convergent if {t^j = t x q < U. 

Now = 1; thus t x q is eaual to U when q = IT and 

= 1 . The term _ t ! q is zero if a = 0 and remains positive 
for 0 < q < U. By the definition of T^t) given hy eauation 
(56), it can he easily shown that 


dq 


t x q 


> 0 


for 0 < T < T » Thus t x q increases mono t oni cally from 
zero to U in the interval 0 5 T <. t x . Therefore, the 
series (92) is uniformly and absolutely convergent in any 
closed domain in |wl < U. 

Similarly, the convergence of the series (9^0 can be 
established. 
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APPENDIX 0 


PROOF OF THEOREM (98) 


It is observed 
among the constants 


that the following identities exist 
involved in (98) and (99): 


B U V — V A U m f-i- + -i-\ m + v) 

n 2vrr Z_ i m ' yn+v m —v/ 


m=o 


Cr. 


c U V = JL y A U m ( -i- + — ^ ( m— v ) 

n 2uvr j m \m+u m— v / 


®*=o 


Now, by the 

S v ( T i>. L V 

all v=f 0, 

Ml*,*’*) is 


inequalities (58-) and (59), the functions 

(t ) can be bounded both above and below for 

1 

when 0 < T < » And if a smaller value of 

2(3+1 

taken, it can be deduced that 


B 


* i 

a ! 


< M x 


_J_B n t 


* i 

c n ; < m s 

** t 



where Mj_ and M s are constants independent of n. On the 
other hand, from the Inequality (58) 


I ( T) < ( 1- 


T) P , 


0 < T < 


1 

2(3 + 1 
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it follows that 


I„CO 

V T !> 


V 

t g( T ) , 


< T 


< 


1 

2^ + 1 


Consequently, the first part of (101) can be dominated: 

| 3 n T > * V 1 < j B n (t ^ w)V ! 

T-( T) 

where t 3 (*r) = — — . The continuation of this inequality 

Te< T i) 

for t > — - — can be easily done by defining a new t 2 (T). 

2^+l 

By hypothesis , T.\ B n (t s w) | converges if i t 2 w j < V. 

n=o 

Since t a (T) < tg(T x ) for T, < T < 1, the inequality 
| t 8 w I < V is uniformly bounded. 

Similarly, it can be shown that 

jo* & (t) wM < jC (tjw) V j 

1 n v \ i a i 

OO 

But i C^ (t]_w) i converges if | t x w 1 > U. Since on 

I w I = U °t,(T.) = 1 and — log i t x w! >0 when 0 < t <—J - — 

d <3. 2f + l 

or — | tiv ( = 0 when — = — < t < 1. the condition 
dq 2^+l 

I t x w i > U holds for all T in T x < T < 1 . Hence, by 

Weiers tras a 1 s theorem the series ( 101) converges uniformly 
and absolutely in U+6<!wi<V-— 6. 
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NACA TN No, 996 


CABLES 0E CHE EYPERGEOKETR IC FUNCTIONS 


The values of the hyper geometr ic functions given 
in tables 1 to 5 are calculated from power series for 

y = 1.405, The function Jj ( t) in table 6 is con— 

— v , i 

nected with j£ y ( t) , Z_ V ( T )» and ]T y ^ t ) through the 

following equation: 

t f (t) i ( T ) 

2( i>— 1 ) v -v , i 


* S (t) S' (t ) 

— U V 


P( v-l) 

2( »+l) T 


f ( T ) | (T) 

, 1 ~v 


(1-t) 




This is s imply the Wr oaskian of the two independent inte- 
grals of the hypergeometric equation and it holds every- 
where except at the s ingular it i es T « 0 and T s l , 

Tables 7 to 12 contain the corresponding approximate func-" 
tlons as indicated. 

The numbers in these tables are expressed in terms of 
appropriate powers of 10. However, a notation was devised 
in which only the powers are given while the base 11 1 0 11 is 

omitted. Thus, 3,14159 X IC 12 & 3*14159, m. Here m nay 
be either a positive or negative integer, or zero. Unless 
indicated by the sign t on the heading, accidental errors 
were detected and eliminated by the difference method. 



NACA TN No. 995 


117 


TABLE 1.- CORRESPONDING- PARTICULAR INTEGRALS FOR THE SOLUTIONS 
OF COMPRESSIBLE FLOW AND INCOMPRESSIBLE FLOW 



The functions |„(t), q _2U G; v (T > and F u (t), q~ 2U Gu(T) 


are respectively the two independent integrals of equa- 
tions (27) and (28)» 

TABLE 14 


















£ 


TABLI 2 


T 

V T) 

w r) 

* 5 / a (T > 

W ” 

^ 9 / 3 ^ 

- 11 / 3 ^ 


- lB / 3 ^ 

E 17 / B ^ 

* 19 / 2 ^ 

0 

1 . 00000 , 

0 

1 . 00000 , 

0 

1 . 00000 , 

0 

l . ooooo . 

0 

l . QOOOO , 

0 

1 , 00000 , 

0 

1 , 00009 , 

0 

1 . 00000 , 

0 

1 . 00000 , 

0 

1 . 00000 , 

•0 

.19 

*. 40592 , 

-1 

B . 26748 , 

-1 

7 . 25608 , 

-1 

6 * 61720 , 

-1 

6 , 60840 , 

-1 

4 , 79680 , 

-1 

4 . 17817 , 

-1 

5 . 63614 , 

-1 

5 . 16 # 59 , 

-1 

2 . 76137 , 

-i 

.12 

9 . 29261 , 

-1 

7 . 94 S 94 , 

-1 

6 . T 5636 , 

-1 

9 . 72109 , 

-1 

4 . 834 B 7 , 

-1 

4 , 08016 , 

-1 

3 * 43976 , 

-1 

2 . 69778 , 

-1 

2 . * 3900 , 

-1 

2 . 06332 , 

-1 

.14 

9 . 19113 . 

-i 

7 , 83846 , 

-1 

6 . 29660 , 

-1 

6 * 16695 , 


4 . 22613 , 

-1 

5 . 44761 , 

-1 

£. 80676 , 

-1 

£. 26561 , 

-1 

1 . 66621 , 

-1 

1 * 60907 , 

-1 

.14 

9 . 12606 , 

-1 

7 . 40607 , 

-1 

6 . 07675 , 

-1 

4 . 90498 , 

-1 

5 . 94290 , 

-1 

5 . 181 E 0 , 

-1 

£. 32960 , 

-i 

2 , 03124 , 

-1 

1 . 61 S 32 , 

-1 

1 , 28841 , 

-1 

.16 

9 . 07143 , 

-1 

7 . 55892 , 

-1 

6 . 66410 , 

-1 

4 . 652 S 2 , 

-i 

5 . 67498 , 

-1 

2 , 09541 , 

-1 

2 . 27271 , 

-1 

1 , 76196 , 

-1 

1 . 39616 , 

-1 

1 . 09110 , 

-1 

.166 

9 . 04429 , 

-1 

7 . 26616 , 

-1 

6 . TB 01 O , 

-1 

440035 , 

-l 

3 , 34821 , 

-1 

2 . 76621 , 

-1 

2 . 16219 , 

-1 

1 . 07109 , 

-1 

1 . 3 * 58 *, 

-1 

1 . 00303 , 

-1 

.17 

9 . 01726 , 

-1 

7 . 19196 , 

-1 

6 . 65605 , 

-1 

4 . 41018 , 

-1 

5 * 42067 , 

-i 

2 . 64535 , 

-1 

2 . 05672 , 

-1 

1 . 60561 , 

-1 

1 . E 0165 , 

-1 

*. 20856 , 

-2 

•ire 

6 . 99036 , 

-1 

7 . 11930 , 

-1 

6 . BB 249 , 

-1 

4 . 29256 , 

-1 

3 , 29831 , 

-i 

2 . 62400 , 

-1 

1 . 93811 , 

-1 

1 . 40665 , 

-i 

1 . 13535 , 

-1 

8 . 44375 . 

-2 

.18 

8 . 96556 , 

-1 

7 . 04719 , 

-1 

6 . 46130 , 

-1 

4 . 17882 , 

-1 

5 . 10004 , 

-l 

2 . 41012 , 

-i 

1 , 82027 , 

-i 

1 . 37103 , 

-i 

1 . 03036 , 

-1 

7 , 72 * 80 , 

-2 

• 186 

8 . 93684 , 

-1 

6 . 97663 , 

-l 

6 . 36 140 , 

-1 

4 . 06356 , 

-1 

3 . 06440 , 

-1 

2 . 2 * 968 , 

-i 

1 . 71898 , 

-1 

1 . 28111 , 

-i 

9 . 52406 , 

-2 

7.06635 

-2 

.19 

0 . 91 * 17 , 

-1 

6 . 9 iH 81 , 

-1 

6 . 25201 , 

-1 

$. 98263 , 

-1 

2 . 06200 , 

-1 

2 . 19288 , 

-1 

1 . 62215 , 

-1 

1 . 19582 , 

-1 

8 . 76286 , 

-2 

6 . 44991 , 

-2 

.196 

8 . 8 * 580 , 

-1 

$. 83411 , 

-1 

6 * 15565 , 

-1 

5 . 84570 , 

-i 

2 . 84274 , 

-i 

2 . 09005 , 

-1 

1 . 92965 , 

-1 

1 . 11525 , 

-1 

8 . 10047 , 

-2 

6 . 90499 , - 

-2 

.20 

8 , 86743 , 

-1 

6 . 76417 , 

-1 

9 . 03 * 82 , 

-1 

3 . 73707 , 

-1 

2 . 75642 , 

*1 

1 . 860 * 3 , 

-1 

1 . 44109 ,' 

-i 

1 . 030 BO , 

-1 

T . 46306 , 

-2 

6 . M 693 , 

-£ 

• 21 . 

B . B 0 T 1 B , 

-1 

6 . 626 BS , 

-l 

4 * 87054 , 

-1 

3 . 33020 , 

-1 

£. 85258 , 

-i 

1 . 00316 , 

-1 

1 . 37901 , 

-1 

8 . 98454 , 

-2 

0 . 29806 , 

-2 

4 . 40048 , 


.22 

0 . 76312 , 

-1 

6 . 48974 , 

-1 

4 . 60017 , 

-1 

3 . 53172 , 

-1 

£. 34028 , 

-1 

1 , 62911 , 

-l 

1 . 12662 , 

-1 

7 . 75225 , 

-2 

5 . 2822 *, 

-2 

3 . 69819 , 

-2 

<23 

8 . 70191 , 

-1 

B . 36671 , 

-1 

4 * 60996 , 

-1 

3 . 14141 , 

-i 

2 . 15882 , 

-1 

1 . 46839 , 

-l 

1 . 0836 E , 

-1 

6 . 618 * 5 , 

-2 

4 . 3 * 813 , 

-2 

2 . 90644 , 

-8 

<34 

6 . 66094 , 

-1 

6 . 22579 , 

-1 

4 . 33882 , 

-1 

3 ** 8901 , 

-1 

1 . 90702 , 

-1 

1 . 51 ( 93 , 

-1 

8 . 18938 , 

-2 

6 . 63346 , 

-2 

3 . 63254 , 

-2 

2 . 326 £ 4 , 

-Z 

.26 

8 . 69991 , 

-1 

6 . 0 * SSB , 

-1 

4 . 18007 , 

-1 

2 . 78498 , 

-1 

1 * 82708 , 

-1 

1 . 18146 , 

-1 

T . 84363 , 

- 2 . 

4 . 76140 , 

-2 

2 . 97302 , 

-2 

1 . 05760 , 

-2 

.23 

8 . 64973 , 

-1 

6 . 96018 , 

- 1 . 

4 . 00426 , 

-1 

2 . 61721 , 

-1 

1 . 67591 , 

-i 

1 . 09480 , 

-1 

8 . 938 ®, 

-2 

5 . 99519 , 

-2 

2 , 40756 , 

-2 

1 * 43192 , 

-£ 

.27 

8 . 4 SSB 9 , 

-1 

6 . 64048 , 

-1 

5 * 64607 , 

-1 

£.* 5734 , 

•1 

1 . 03401 , 

-1 

9 . 30418 , 

-2 

5 . 6 KE 4 , 

-2 

5 . 333 * 9 , 

-2 

1 . 92 MZ , 

-2 

1 . 09662 , 

-2 

. 28 . 

8 . 46066 , 

-1 

S . 718 B 2 , 

-1 

5 . 69040 , 

-1 

£. 30467 , 

-1 

1 . 40094 , 

-1 

9 , 31 $ 00 , 

-2 

4 . 02837 , 

-3 

2 . 73684 , 

-3 

1 . 61785 , 

-2 

8 , 20153 , 

-5 


8 * 40177 , 

-X 

8 . 6 * 519 , 

-1 

5 . 54054 , 

-1 

£. 16009 , 

-1 

1 . ETB 30 , 

-i 

7 . 55811 , 

-2 

4 . 10444 , 

-2 

2 . 33087 , 

-2 

1 . 17400 , 

-2 

5 * 96261 , 

-5 

.30 

8 . 36332 , 

-1 

B . 47387 , 

-1 

2 U 39465 , 

-1 

2 . 01961 , 

-1 

1 . 13 * 69 , 

-1 

6 . 44493 , 

-2 

5 . 48251 , 

-2 

1 , 79225 , 

-2 

9 . 875 **. 

-5 

4 . 1 * 079 , 

-5 

.32 

8 . 26771 , 

-1 

6 . 24232 , 

'1 

5 * 11619 , 

-1 

1 . 76056 , 

-1 

9 . 49764 , 

-2 

4 * 09140 , 

—2 

2 . 3 B 164 , 

-2 

1 . 09548 , 

-2 

4 . 5694 T , 

-5 

1 * 65491 , 

-3 

.34 

9 . 19381 . 

-1 

6 . 01802 , 


2 . 88491 , 

-1 

1 * 62658 , 

-1 , 

7 . 67759 , 

-2 

5 . 61699 , 

-2 

1 , 56592 , 

-2 

5 . VT 06 B , 

—5 

1 * 75030 , 

-5 

1 * 01766 , 

-4 

• 36 

8 . 07192 , 

-1 

4 . 79026 , 

-1 

2 . 60899 , 

-1 

1 * 31370 , 

-1 

6 . 11164 , 

-2 

2 * 68167 , 

-2 

9 . 45060 , 

-3 

2 * 55652 , 

-5 

5 * 66667 , 

-e 

- 5 . 94075 , 

-4 

.36 

7 . 96111 , 

-1 

4 . 66 B 25 , 

-1 

2 . 37900 , 

-1 

1 . 12320 , 

-1 

4 . 77666 , 

-t 

1 . 75608 , 

-2 

4 . 66806 , 

-3 

2 . 90912 , 

-4 

- 9 . 06462 , 

-4 

- 8 * 00554 , 

-4 

*40 

7 . 88 * 20 , 

-1 

4 . 56070 , 

-1 

2 . 18418 , 

-1 

9 . 62000 , 

-2 

5 . 64802 , 

-2 

1 , 10041 , 

-2 

1 , 62290 , 

-2 

- 1 , 00371 , 

-5 

- 1 * 61917 , 

—6 

- B . M 30 T , 

-4 

.42 

7 . 80612 , 

-1 

4 . 19614 , 

-1 

1 . 963 S 8 , 

-1 

8 . 00764 , 

-2 

E , 70013 , 

-2 

6 . 06927 , 

-3 

* 4378 , 

—4 

- 1 . TT 029 , 

-5 

- 1 * 58955 , 

-5 

- 7 * 6251 ?, 

-4 

.44 

7 . 71962 , 

-i 

4 . 00562 , 

-1 

1 * 77706 , 

-1 

6 * 66119 , 

-2 

1 . 911 T 8 , 

-2 

£. 36821 , 

-3 

- 1 . BQB 72 , 

—3 

- 1 . 02384 , 

-5 

- 1 . 25888 , 

-5 

- 5 , 83137 , 

-4 

.46 

7 . 63677 , 

-i 

3 . 82294 , 

-1 

1 . 80388 , 

,-1 

6 . 47482 , 

-2 

1 . 27950 , 

-2 

1 . 30069 , 

-4 

- 2 . 63 * 49 , 

-3 

-IV 7881 , 

-5 

- 1 . 02127 , 

•8 

-3 , 58237 , 

-4 

.48 

7 . 36566 , 

-1 

3 . 64809 , 

-1 

1 . 44504 , 

-1 

4 . 45664 , 

-2 

7 , 87404 , 

-a i 

' - 2 . 12640 , 

-3 

- 2 . B 02 B 4 , ' 

-5 

- 1 . 74616 , 

-5 

- 7 , 8 * 880 , 

-4 

— 1 * 7555 ? , 

-4 

~ 6 Q 

7 . 47296 , 

-l 

3 . 48202 , 

-1 

1 . 27464 , 

-1 

5 . 65204 , 

-2 

3 . 66270 , 

-3 

- 4 * 02362 , 

-3 

- 2 . 80466 , 

-3 

- 1 . 42691 , 

—5 

- 4 . 95182 , 

-4 

- 2,925 ?, 

-5 



1 . 00000 , 
2 *$ 9 e ^, 

; 

1 . 82678 , 
1 . 02406 . 
B . 6246 T , 
7 . 76607 . 
7 . 04764 . 
6 . 59286 , 
6 . 78904 , 
fi . 25290 , 
4 . 7212 *, 

4 . 26122 .. 
2 . 61990 , 

2 .06207, 
2.42151, 
1.90756, 
1.47590, 
1.12607, 
B. 68040, 

B .09524, 
4.2B202, 
2.86751, 
1.76066, 
5.05966, 

-3.44649, 
-6.16166, 
-6.26445, 
- 6 . 21186 , 
-5.B205S, 
-2.54660, 
-1.12716, 
-2.46121 , 
+*. 98821 , 


0 

-1 

-1 

-1 

-1 

-2 

—2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 


► 

Q 

► 


S 


• -i • \ 


VV V V ? 7 ? ? •? T T t T 1 T 1 TT V 



V 


I 

TABLE 3 


W T) 


*V T> 


~-5/3^ 


*V T> 


%/3 lT) 


■2-11/3^ T ' 


• 10 
*12 
•14 
.15 
-lfi 
.166 | 

• IT 

• 1TB | 
.18 
.186| 
•16 
.165 | 
.£0 
.21 

.SB 

.23 

•24 

•25 

•26 

•2T 

•28 

.29 

•80 

•82 

•84 

.80 

•88 

.40 

.48 

•44 

.46 

•48 

.60 


l.oodco, 

1 .,06389, 
1.0692B, 
1,07998, 
1.08610, 
1.06081, 
1.068 BB, 
1.08648, 
1.09796, 
1.10048, 
1.10898, 
1.10646, 
1.10793, 
1.11066, 
1.11626, 
1.12006, 
1*12481, 
1.1W60, 
1.16414, 
1.13871, 
1.14828, 
1.14770, 
1.15211, 
1.15640, 
1.16500, 
1.1T588, 
1.18146, 
1.1M56, 
1*16700, 
1.20456, 
1.21187, 
1.21899, 
1.22591, 
1.28255, 


1.00000, 
1.12065, 
1.13121, 
1.15814, 
1.14060, 
1.14168, 
1.14201, 
1.14218, 
1.14216, 
1.14196, 
1.14101, 
1.14106, 
1. 14080, 
1.15649, 
1.18724, 
1.15489, 
1.18095, 
1.12695, 
1.12244, 
1.11745, 
1.11196, 
1.10605, 
1.06978, 
1.09804, 
1.07801, 
1.00296, 
1.04080, 
1.02681, 
1.01064, 
6.91660, 
9.72948, 
9.58707, 
9.84406, 
9.10148, 


a 

o 

a 

o 

o 

o 

o 

o 

o 

Q 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

-1 

-1 

-1 


1.00000, 
1 .41771, 
1.49211, 
1 *66566, 
1.57687, 
1.86695, 
1.00889, 
1.61878, 
1.0184B, 
1.62910, 
1.68886, 
1.08087, 
1.68608, 
1.64000, 
1.68842, 
1.65210, 
1.62166, 
1.GO680, 
1.68605, 
1.66010, 
1.58082, 
1 .49786, 
1 .40973, 
1*41012, 
1.82801, 
1.21606, 
1.006B2, 
8.6T676, 

H. 81023, 
0.69672, 
0.44047, 
8.96089, 
2.68280, 

I. 10727, 


1.00000, 0 

1.T20O8, 0 

1.90884, 0 

2,07768, 0 

£.16207, Q 

2.21020, 0 

2.24003, 0 

2.87677, 0 

2.19904, 0 

1.82087, O 

2.88871, 0 

2.86820, 0 

1.36480, 0 

1.87172, 0 

2.87626, 0 


2.30298* 


2.58411, 0 

2.26610, - 0 
2.22477, 0 

2.14878, 0 

2.04004, 0 

1.82B64, 0 

1.76545, 0 

1.04685, 0 

1.29700, 0 

8.98074, -4. 

4.40206, -1 

-6.00266, -2 

-6.65467, -1 

-1.09306, 0 

-1.61928, 0 

-2.18027, 0 

*2.61298, 0 

-3.06472, 0 


1.00000, 0 

2.06116, 0 

2.56596, 0 

2.72468, 0 

2.86460, 0 

8.08204, 0 

3.10056, 0 

8.16514, 0 

3.22029, 0 

8.27119, 0 

3.81624, 0 

8.30162, 0 

3.86082, 0 

3 .40014, 0 

3.41159, 0 

8.86194, 0 

3.80745, 0 

5.16493, 0 

5.01161, 0 

2.78019, 0 

2.60098, 0 

2.17366, 0 

1.78068, 0 

1.54769, 0 

3.21F07, -1 

-6.76444, -1 

-2.21429, 0 

-5.68755, 0 

-6.06175, 0 

-6.61410, 0 

-7.84E0T, 0 

-9.01412, 0 

-9.96614, 0 

-1.06646, 1 


1.00000, 
2.4OST0, 
2.98085, 
8 .02660, 
8.01970, 
4.10001, 
4.28287, 
4.36026, 
4.47061, 
4.67866, 
4.04869, 
4.73692, 
4.79766, 
4.08764, 
4.06625, 
4.76045, 
4.56764, 
4.29654, 
8.06665, 
8.80162, 
2*39526, 
1.74670, 
6.53468, 
-8.77564, 

-5.02436, 

-0.09000, 

-9,40146, 

-1.20941, 

-X.09C38, 

- 1 . 91580 , 

-2.14604, 

-2.29447, 

— 234069 , 

-2.27281, 



t 

•2-13/3 

(t) 

t 

—16/3 

(T) 

t 

--n/a 

(T) 

t 

—19/2 

(t) 

--21/2 

(t) 

0 

1.00000, 

0 

1.00000, 

0 

1.00000, 

0 

1.00000, 

0 

1.00000, 

0 

0 

2.79840, 

0 

3.28017, 

0 

3.72368, 

0 

4.23471, 

0 

4.98498, 

0 

0 

3*60332, 

0 

4.88940, 

0 

5.22992, 

0 

6.26099, 

0 

7*52361, 

0 

0 

4.51962, 

0 

E.T8468, 

0 

7.26997, 

0 

9.19794, 

0 

1*16606, 

1 

0 

6.01470, 

0 

6,64263, 

0 

8-49624, 

0 

1.09096, 

1 

1.41735, 

1 

0 

5.60301, 

0 

T.34426, 

0 

9.76040, 

0 

1.29267, 

l 

1*70790, 

1 

0 

5.73721, 

0 

T.T3T04, 

0 

1.08932, 

1 

1.39181, 

1 

1,06944, 

1 

0 

6.96026, 

0 

B. 11670, 

0 

1.10104, 

1 

1.46961, 

1 

*,01070, 

1 

0 

6.16B44, 

0 

0.47290, 

0 

1.15990, 

1 

1.38567, 

1 

2*15772, 

1 

0 

0.86762, 

0 

0.60060, 

0 

1.21426, 

l 

1.67116, 

1 

2.29559, 

1 

0 

0,62420, 

0 

9.0S9B6, 

0 

1.28236, 

1 

1.74878, 

1 

2.41773, 

1 

0 

4.06345, 

0 

9.33005, 

0 

1.30229, 

1 

1.01209, 

1 

2.61778, 

1 

0 

6.7T089, 

0 

9.61862, 

0 

1.33197, 

1 

1.06906, 

1 

2.68765, 

1 

0 

6.B4198, 

0 

9.62201, 

0 

1.84916, 

1 

1*00307, 

1 

2.61861, 

1 

0 

6,85711. 

0 

9.60926, 
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-l 

4.92479, 

-1 

2.09491, 

-1 

1.87110, 

-l 

1,10237, 

-1 

*.31633, 

•10 

7.26975, 

-l 

4.6 T 960, 

-1 

2,98279, 

-1 

1 .7 X 084, 

-l 

9.60291, 

-2 

*.42840, 

.» 

8.70161, 

-l 

4.25917, 

-1 

2.92964, 

-1 

1*4388®, 

-1 

7.6646 B , 

4 

3 .8*380, 

.64 

6.20992, 

-i 

5.86159, 

-1 

2.20917, 

-1 

1,18510, 


6*92761, 

-2 

2.73131, 

•60 

5.76927, 

-l 

5.50441, 

-1 

1.92896, 

-1 

9.73488, 

-% 

4*47178, 

-2 

1.(2367, 

.16 

5.56904, 

-l 

5.18978, 

-1 

1.67620, 

-1 

7,90176, 

-2 

5.27109, 

-2 

1.11127, 

•40 

6.00 X 66, 

-l 

2.90193, 

-1 

1-49389, 

-1 

6,52097, 

-2 

2.39030, 

-2 

6.72187, 

•42 

4*68159, 

-l 

2.65899, 

-1 

1.29202, 

-1 

4,96805, 

-2 

1.49655, 

-2 

I . TT 128, 

•44 

4.54424,' 

-l 

2.59199, 

-1 

1.07046, 

-1 

5.90819, 

-3 

. B . 89751, 

«3 

.1,00216, 

•46 

4,04898. 

-i 

2.16475, 

•1 

9.07280, 

-e 

2,81887, 

-2 

5.81381, 

-0 

.2.82746, 

•46 

5.78255, 

-i 

1.96561, 

-1 

6.74099, 

-2 

1.90003, 

-3 

1.52100, 

-4 

-1.90216, 

•60 

5.60250, 

-i 

1.76684, 

-1 

6.42218, 

-2 

1.29782, 

-2 

-2.55621, 

4 

-4.39783, 





nmmm i ^hnmnnnn 



r 


>- 

o 

► 


TABU 9.- /^(T) - f(T) I~ D {T)—>lf( T ) l“ V (T) 00. (w» + 


T 



4V> 

#/ 8 < t > 

(0) 

£-9/3 

(T) 

—11/3^ 

>) 

£-13/3 

(T) 

r {°) 

£-15/3 

(T) 

,(o) 

£-17/3 

(7) 

-<0> 

£-ib/8 

(7) 

r (o) ( T ) 

--21/3 

0.02 

1.0145b, 

0 

1.04031, 

0 

3.06603, 

0 

1.09425, 

0 

1.12225, 

0 

1.16004, 

0 

l.iaoso. 

0 

1.21060, 

0 

1.24167, 

0 

1.27334, 

0 

1,30598, 

0 

0.04 

1,03534, 

0 

1.06657, 

0 

1.14645, 

0 

1.20756, 

0 

1.27165, 

0 

1.35010, 

0 

1.41027, 

0 

1.48523, 

0 

1.66416, 

0 

1.84T29, 

0 

1.734 64, 

0 

O.Ofl 

1.06986, 

0 

1.14763, 

0 

1.24310, 

0 

1.34*2*, 

0 

1.45802, 

0 

1.67004, 

0 

1.71010, 

0 

1.85606, 

0 

2.00677, 

0 

2.17225, 

0 

2.36156, 

0 

0.03 

1.09659, 

0 

1.22210, 

0 

1.56443, 

0 

1.52217, 

0 

1.69809, 

0 

1.80436, 

0 

2*11530, 

0 

2.35766, 

0 

2.63003, 

0 

2.93400, 

0 

3.27310, 

0 

0.10 

1.14952, 

0 

1.52362, 

0 

1,52410, 

0 

1.76493, 

0 

2.02073, 

0 

2.32879, 

0 

2.67920, 

0 

1.06499, 

0 

3.55224, 

0 

4.09025, 

0 

4.70075, 

0 

0.12 

1.25522, 

0 

1.48225, 

0 

1,76068, 

0 

2.06699, 

0 

2.40546, 

0 

2.06146, 

0 

3.62662, 

0 

4.20441, 

0 

5.00962, 

0 

5.96906, 

0 

7.11221, 

0 

0.14 

1.58760, 

0 

1.71669, 

0 

2.12375, 

0 

. 2.62744, 

0 

3.26062, 

0 

4.02181, 

0 

4.97646, 

0 

6.16638, 

0 

7.51553, 

0 

0 .42179, 

0 

1.16686, 

1 

0.16 

l.ftf666. 

0 

1.95944, 

0 

2.44959, 

p 

3.00543, 

0 

3.90680, 

0 

4.95446, 

0 

8.23139, 

0 

7.8608ft, 

0 

9,98913, 

0 

1.25526, 

1 

1.68550, 

1 

0.16 

1.86559, 

o' 

2.79764, 

0 

5.09669, 

0 

3.69966, 

0 

6.16665, 

0 

6.67170, 

0 

8.61692, 

0 

1.11203, 

1 

1*43741, 

1 

1.8564ft, 

1 

2.39777, 

1 

0.166 

Z.7J5115, 

0 

3.05463, 

0 

5.96607, 

0 

6.18822, 

0 

6.78566, 

0 

6.86042, 

0 

1.15967, 

1 

1.51526, 

1 

1.982+9, 

1 

2.69209, 

1 

3.38012, 

1 

0.17 

1.07704, 

0 

2.62022, 

0 

5.47696, 

0 

4.60257, 

0 

6.00966, 

0 

8.0B400, 

0 

1.07140, 

i 

1.41996, 

1 

1.86169, 

1 

2.49412, 

1 

3.30E51, * 

I 

D.17S 

1.59531, 

0 

1.56658, 

0 

2.46904, 

0 

3.53870, 

0 

*.47831, 

0 

6.00686, 

0 

8.06606, 

0 

1. 0806ft, 

1 

1.44044, 

1 

1.04403, 

1 

2.60737, 

1 

0.18 

1.19687, 

0 

1.62226, 

0 

2.19664, 

0 

2.97967, 

0 

4.03762, 

0 

6.47065, 

0 

7.41186, 

0 

1.00*10, 

1 

1.35014, 

1 

1.34229, 

1 

2,49608, 

1 

0.186 

1.06951, 

0 

1,40040, 

0 

2.05866, 

0 

2.78761, 

0 

3.61068, 

0 

6.20810, 

0 

T. 11616, 

0 

9.72066, 

0 

1.327/18, 

1 

1.81236, 

1 

*2.47302, * 

1 

0.15 

1.P1602, 

0 

1.50082, 

0 

1.92961, 

0 

2.66922, 

0 

3.66014, 

0 

6.03680, 

0 

6.92715, 

0 

9,52118, 

0 

1.50782, 

1 

1.79621, 

1 

2.46255, 

1 

0.1S5 

9.58226, 

-1 

1,33082, 

0 

1.84661, 

0 

2.66947, 

0 

3.64596, 

0 

4.00170, 

0 

6.77165, 

0 

9.34346, 

0 

1.28762, 

1 

1.77174, 

1 

2.43446, 

1 

0.20 

9.12566, 

-1 

1.27494, 

0 

1.77669, 

0 

2.47727, 

0 

3.44596, 

0 

4.77857, 

0 

8.51868, 

0 

9.14090, 

0 

1.2507ft, 

1 

1.73152, 

1 

2,37306, 

1 

0.21 

9.40756, 

-1 

1.18858, 

0 

1.66788, 

0 

2.35714, 

0 

3.26928, 

0 

4.50251, 

0 

8,24248, 

0 

6.66138, 

0 

1.18589, 

1 

1.57371, 

1 

2.10311, 

1 

o.zz 

7.87016, 

-1 

1.11600, 

0 

1.67419, 

0 

2,80307, 

0 

3.06776, 

0 

4.19690, 

0 

6.58606, 

0 

7.65023, 

0 

0*90*10, 

0 

1.25728, 

1 

1.63635, 

1 

0.2* 

7.29586, 

-l 

1.06451, 

0 

1.48679, 

0 

2.06469, 

0 

2.81962, 

0 

3.78601, 

0 

4.90319, 

0 

6.06090, 

0 

7.064TJ5, 

0 

7.50396, 

0 

T. 107+0, 1 

0 

O.Z* 

7.00080, 

-1 

9.96920, 

-1 

1.50893, 

0 

1.91522, 

0 

Z. 686 30, 

0 

3.21210, 

0 

3.79363, 

0 

3.03031, 

0 

2.05415, 

0 

-+,99959, 

-1 

-2.81*31, 

Q 

0.26 

6.66402, 

-1 

9.472 BO, 

-1 

1.51099, 

0 

1.74512, 

0 

2. 19851, 

0 

2.52055, 

0 

Z. 38206, 

0 

1.10702, 

0 

-2.51476, 

0 

-1.12276, 

1 

-2*91660, 

1 

O.20 

6.24156, 

-X 

Q. 96065, 

-1 

1.22046, 

0 

1.58129, 

0 

1.80446, 

0 

1.66068, 

0 

6.23766, 

0 

“8.48611, 

0 

-9.77446, 

0 

-2.01542, 

1 

-5.38ft TO, 

1 

0.2 1 

6.0BP84, 

-1 

8.60627, 

-1 

1.1287b. 

0 

1.56793, 

0 

1.55267, 

0 

6.89015, 

-1 

-1.4P096, 

0 

-6.86430, 

0 

-1.86760, 

1 

-4.23646, 

1 

-0.847*8, 

1 

0.2ft 

8.78147, 

-1 

8.06995, 

-1 

1.02956, 

0 

1.15606, 

0 

8.42660, 

-1 

-4.02905, 

-1 

-3.96362, 

0 

-1.20262, 

1 

-2.89045, 

1 

-6.27815, 

1 

-1,27228, 

2 

0.28 

6.64476, 

-1 

7.66697, 

-1 

9.28626, 

-1 

0.06666, 

-1 

Z. 78121, 

-1 

-1.78800, 

0 

-6.77226, 

0 

-1.79149, 

1 

-4.08776, 

1 

-8.6B920, 

1 

-1.70778, 

? 

0.20 

6.21244, 

-1 

7.15734, 

-1 

0.25600, 

-1 

6.58876, 

-1 

-3.42850, 

-1 

-3.18572, 

0 

-9.90920, 

0 

-2.44750, 

1 

-8.4065ft, 

l 

-1.11232, 

2 

-2.17616, 

2 

0.22 

4.68404, 

-1 

6.26910, 

-1 

8.0*706, 

-1 

7.707OC, 

-2 

-1.72*09, 

0 

-*. 37915, 

0 

-1.60446, 

1 

-3.906*9, 

1 

-8.28015, 

1 

-1.6*060, 

2 

-3.09169, 

7 

0.34 

4.49410, 

-1 

6.56627, 

-1 

3.74912, 

-1 

-6.54046, 

-1 

-3.26621, 

0 

-9.03145, 

0 

-2.34785, 

1 

-6.4889ft, 

1 

-1.10463, 

2 

-2,09*62, 

2 

-8.73929, 

2 

0.26 

4.1381P, 

-1 

4.64018, 

-1 

1.58204, 

-1 

-1.17927, 

0 

-4.871C7, 

0 

-1.36061, 

1 

-3.21568, 

1 

-6.80790, 

1 

-1.31890, 

2 

-2.33893, 

2 

-3.73154, 

Z 

0.26 

3.80168, 

-1 

3.71165, 

-1 

-1.01642, 

-1 

-1.830*6, 

0 

-0.49930, 

0 

-1.71478, 

1 

-3.8775ft, 

1 

-7.80525, 

1 

-1.40365, 

2 

-2.19883, 

2 

-2.72673, 

2 

0.40 

5.4 9107, 

-1 

2.90607, 

-1 

-3.40012, 

-1 

-2 .40511, 

0 

-8.08630, 

0 

-2.02670, 

1 

-4.36049, 

1 

-ft.l8R40, 

1 

-1.29898, 

2 

-1.54454, 

2 

-4.99521, 

1 

0.42 

3.19960, 

-1 

2.15287, 

-1 

-6.72770, 

-1 

-3.110*2, 

0 

-9.40830, 

0 

-2.26740, 

1 

-4.87657, 

1 

-6,48765, 

1 

-B, Pi 26 30, 

1 

-3. 14716, 

i 

2.931*0. 

2 

0.44 

2.92620, 

-1 

1.59091, 

-1 

-7.96*20, 

-1 

-5.6BWR, 

0 

-1.06211, 

1 

-2.40969, 

1 

-4 .46007, 

1 

-6*18020, 

1 

-2.5823ft, 

0 

-1.43021, 

2 

7.1843ft, 

2 

0.46 

2.66912, 

-1 

6.87472, 

-2 

-1.00362, 

0 

-4.20472, 

0 

-1.14689, 

1 

-2.43071, 

1 

-3.93644, 

1 

-3,64254, 

1 

4.5*826, 

1 

3.52680, 

2 

1.14148, 

3 

0.46 

2.42713, 

-1 

2.69076, 

-5 

-1.65868, 

0 

-4 .62620, 

0 

-1.19376, 

1 

-2,31634, 

1 

-2,08528, 

1 

-1.66651, 

0 

1.43301, 

2 

6.81980, 

2 

1.48ft 80, 

3 

0,60 

2.19986, 

-1 

-6.86492, 

-2 

-1.67610, 

0 

-4.94566, 

0 

-1.1984C, 

1 

-2.06754, 

1 

-1.72386, 

1 

-4. 00740, 

1 

2.48137, 

£ 

7.29125, 

2 

1. 8705ft, 

8 
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TABU 10.- f£ 0, (t) - s(t) ^(t) — ►sCt) i”(t) <x» (vm - J) 


S (o) r i 

13 / 8 ( 7 ) 

;( o ). , 
Zs / s ( T > 

;( o ), 

l 7/ a( T 

) 

* 0 / 8 ( 7 ) 

£ ll /2( T ) 

; ( o ) 

hs / s ^ 

) 


£ i 7/ b( t 

) 

=(°) 1 1 
Eig / aff ) 

f ( o ) (T 
- a i / s ' 

) 

1.0 SHU * 

0 

1.01389, 

c 

9.88402, 

-1 

9.63742, 

►1 

9.19697, 

-1 

9.16231, 

-1 

8,93390, 

-1 

S . T 110 X , 

-1 

8.493*6. 

-1 

*.29176,. 

-1 

5.07311, 

-1 

1.06660, 

0 

1.03078, 

0 

9. T 870 X , 

-1 

9.29566, 

-1 

8.82469, 

-1 

8.87985. 

-1 

T . 95848, 

-1 

7.66496, 

-1 

%. 17570, 

-1 

0.81168, 

-1 

S .46793, 

-1 

1.14018, 

0 

1.06279, 

0 

9.72109, 

-1 

$.97606, 


B .2 B 813, 

-1 

7.85294, 

-1 

7,00842* 

-1 

6,5 X 488, 

-1 

6 .02479, 

-l 

5.16906, 

-1 

6,15669. 

-1 

1,20746. 

0 

1.08257, 

0 

9.70230, 

-1 

8*98712, 

-1 

7.79607, 

-1 

6.96838 , 

-1 

6,26456, 

-1 

5*61535, 

-1 

5*03359, 

-1 

4.31209, 

-1 

4.044*3, 

-1 

1,29496, 

0 

1.12461, 

0 

9. TB 88 B , 

-1 

8*48 X 17, 

-1 

7.33046, 

-1 

6.59700, 

•1 

6.66600, 

-1 

4.8 X 518, 

-1 

4.19060, 

-1 

5.039 X 9, 

-1 

3.18060, 

-1 

1.41902, 

0 

1.19084. 

0 

9.99612. 

-l 

8.38888, 

-1 

T . 04026, 

-1 

E . ID 8 K , 

-1 

4.90*93, 

-1 

4.16138, 

-1 

5.49292, 

-1 

X . 95180, 

-1 

2,46031. 

-1 

i .« tr &4. 

0 

1,61663, 

0 

1.06555, 

0 

8.89476, 

*1 

6.94705, 

-1 

Sr .1522, 

•1 

4,65872, 

-1 

3.66060, 

-1 

£.96629, 

-1 

2.59681, 

-1 

1*95732, 

-1 

1.81530, 

0 

1.43828, 

0 

1.13883, 

0 

9.01060, 

-1 

7.1 B 959, 

-1 

5.85301, 

•1 

4.47606, 

-1 

5.54418, 

-1 

2,50630, 

-1 

2.22204, 

-1 

1.75043, 

-1 

Z *21066, 

0 

1,71181, 

0 

1.52640, 

0 

1.0 X 608, 

0 

7*94446, 

-1 

• 05694, 

-1 

4. T 8255, 

-1 

5.68745, 

-1 

£.05505, 

-1 

2. £1063, 

fl 

1,70998, 

-1 

Z . T 7097, 

0 

£.11930, 

0 

1.82089, 

0 

1.23970, 

0 

9.48160, 

-1 

7.281*8, 

-1 

5.34623, 

-1 

4.24191, 

-1 

3,24452, 

-1 

2.48135, 

-1 

1,89778, 

-1 

4,72857. 

0 

3*66034, 

0 

2*69061, 

0 

2.03030, 

-0 

1. 85193, 

-0 

1.183*7, 

0 

8.72144, 

-1 

6-88217, 

-1 

4.98853, 

-1 

5.74629, 

**1 

X .8 E 683, 

-1 

5.51 C 4 T , 

0 

2 .47233, 

0 

1.84300, 

0 

1.37586, 

0 

1.02415, 

-O 

7*83405, 

-l 

8.69069, 

-1 

4.06090, 

‘1 

3.18185, 

-1 

2.36670, 

-1 

1,768*9, 

*1 

£.88255. 

0 

2.12636, 

0 

1*66840, 

0 

1.1567 T , 

0 

B .03116, 

-1 

6.29112, 

-1 

4.65889, 

-1 

5.42050, 

-1 

2.52177, 

-1 

1.86906, 

-1 

UBTOtl . 

-.1 

£.83693, 

0 

1.92812, 

0 

1.40646, 

0 

1.0 X 865, 

0 

7*51078, 

-1 

•*4&505, 

-1 

4.00186, 

-1 

E . 92026, 

-1 

£.19061* 

-1 

1.65415, 

-1 

1.11 MB , 

-1 

2.47544, 

0 

1.78311, 

0 

1.29637, 

0 

9.39669, 

-1 

6.79759, 

-1 

4.91554. 

-1 

8.332*3. 

•1 

2.5686 T , 

-1 

1.85581, 

-1 

1.33056, 

-1 

9.** *40, 

«-£ 

2.56654, 

0 

1.69 X 67, 

0 

1.21001, 

0 

8.7 X 516, 

-1 

6.25456, 

-1 

4*48075, 

-1 

5.20796* 

-1. 

2.2946 B , 

-1 

1.64078, 

-1 

1.17 X 2 S , 

-1 

8.38970, 

-2 

2.25219, 

0 

1.61415, 

0. 

1.16030, 

0 

8.18491, 

•1 

8.91745, 

-1 

*.129*6, 

-1 

£.92841, 

-1 

2.074 X 9, 

-1 

1.46774, 

-1 

1.03749, 

-1 

7*32*90, 

-£ 

2.125*9, 

0 

1.49637, 

0 

1.06068, 

0 

T .35 646, 

-1 

5.15841, 

•1 

3-58072, 

-1 

£.46968, 

-1 

1.72742, 

-1 

1.19610, 

-1 

8.2 C 695, 

-X 

6*70166. 

-2 

ZJto 900, 

0 

1.41171, 

0 

8*77366, 

-1 

6.73567, 

-1 

4.83276, 


3.16014, 

-1 

2.16251, 

-1 

1.460 T 6, 

-1 

9.88151, 

-1 

6,66 X 88, 

-2 

4*47862, 

-X 

1.96765, 

0 

1.54739, 

0 

0.20208, 

-1 

6.2419 C , 

*1 

4.20366, 

-1 

2 .6 1732, 

-1 

1.87089, 

-1 

1.243*4, 

-1 

8.201 TO , 

-2 

6.58186, 

-2 

*.*1003, 

-2 

1.90304, 

0 

1.29631, 

0 

8*73889, 

-1 

6*62184, 

*1 

3.88226, 

-1 

2.52566, 

-1 

1,64 X 94, ■ 

-1 

1.06088, 

-1 

«. S 0 l 40, 

-X 

4*33992, 

-2 

2.73733, 

-X 

1*86048, 

0 

1.26642. 

0 

8.34283, 

-1 

5.47005, 

-1 

3.54363, 

-1 

2.2 T 014, 

-1 

1.43920, 

-1 

. 9.03263, 

-2 

8-61301, 

-2 

3-46341, 

-2 

2.10302, 

-z 

1.82895, 

0 

1.22137, 

0 

8.00298, 

-1 

5.15565, 

-1 

3.28748, 

-1 

2.04201, 

-1 

1.35678, 

-1 

7,6*623, 

-fi 

4,694 £7, 

-Z 

Z . 71871, 

-X 

1.58618, 

-r 

1.80166, 

0 

1.18142, 

0 

7.70731, 

-1 

4.87004, 

-1 

3.01787, 

-1 

1.83630, 

-1 

1.007*4, 

-1 

6*44880, 

-E 

5. T 190 E , 

-£ 

2.10285, 

-t 

1.18202, 

-z 

1.77972. 

0 

1.16852, 

0 

T .432 T 3, 

-1 

4,60400, 

-1 

£.78408, 

-1 

1.04390. 

-1 

9-51444, 

-2 

6.37404, 

-2 

2.96016, 

-X 

1.58448, 

-X 

8.19 X 66, 

-5 

1.70572, 

0 

1.10828, 

0 

6.94- lfT , 

-1 

4.21076, 

-1 

2.48080, 

-1 

1*42144, 

-1 

7,91670, 

-2 

4.27697, 

-X 

2,23181, 

-X 

1.11490, 

-X 

5. ZB 886, 

-3 

1.761*4, 

b 

1.12896, 

0 

6 *96719, 

-1 

4.1 X 747, 

-1 

2.30499, 

•1 

1.50964, 

-1 

6.99837, 

-X 

3.58860, 

-2 

1.76073, 

-2 

7.97632, 

-5 

3.2*045, 

-5 

1.75875, 

0 

1* 09860, 

0 

6.5*014, 

•1 

5.69786, 

•1 

1.99132, 

-1 

1.01940, 

-1 

4,920 TB , 

-2 

2.19926, 

-X 

S . T 2 TST , 

-3 

2.71854, 

-5 

X . 80349, 

-5 

1*73244, 

0 

1.07462, 

0 

6.16036, 

-1 

5.50061, 

-1 

1.65588, 

-1 

7.87615, 

-£ 

3.22470, 

-2 

1.15647, 

-2 

2. 673 64, 

-5 

-6.74978, 

-5 

-1.32606, 

-3 

1.73858, 

0 

1.05485, 

0 

6.90438, 

-1 

2.92719, 

-1 

1.54599, 

-1 

B -49526, 

-2 

1.83*19, 

-2 

3.4*102, 

-5 

-1.4*361. 

-5 

-2.407*3. 

-5 

>1.99136, 

-3 

1.74737, 

0 

1.05842, 

0 

6*48362, 

-1 

2.57239, 

-1 

1.08489, 

-1 

3.8 1771, 

-2 

7.42911, 

-3 

-2,11660, 

-8 

-3.91146, 

-5 

-5.19604, 

-3 

-2.0*435, 

-ft 

1,7011, 

0 

1.0 X 438, 

0 

8.16202, 

-1 

£.23507, 

-1 

8.06441, 

-3 

2.01980, 

-a 

- l . ooeos . 

-3 

*6,01232, 

-5 

-5.06977, 

-5 

-3,28674. 

-3 

>1.75*4*. 

-ft 

1.76914, 

° 

1.01*13, 

.0 

• 4.80363* 

-1 

1*90705,1 

-1 

5.75092, 

-2 

6. T 985 T * 

-s 

-7.24486, 

-3* 

-2,62393* 


-5.54367, 

-ft 

-2. SM 21, 

-5 

-1.811 TC * 

> -ft 

1.81853, 

0 

1.00119* 

0 

4,48061, 

-1 

1.89312, 

“1 

3.64048, 

-2 

•4. UN 1, 

-3 

-1.11778* 

-2 

—8.86106, 

-ft 

-6.96392, 

-ft 

-2.25619, 

-5 

*8,87482* 

-4 

1.66306, 

0 

9.91200, 

-l 

4,15491, 

-1 

1.29083, 

-1 

1*74711, 

-2 

•1.29466* 

-Z 

-1,41383, 

-2 

-8.60360* 

-5 

-4.16712, 

■X 

-1.63572, 

*8 

*3, *9673, 

-4 

1.89097, 

0 

9 .81 546, 

-1 

3.58806, 

-1 

8.99196, 

-2 

6.65945, 

-4 

-1.9*124, 

-X 

-1.67499, 

-2 

-8.05296, 

-5 

-3.14*23, 

-5 

-8.31431, 

-4 

-1.25787* 

**6 

1.93927, 

0 

9.72788, 

-1 

5.0078 T , 

-1 

7.1*620, 

-2 

-1 .40578, 

-£ 

-£.43514, 

-2 

-1.54161, 

-X 

’-6*78851, 

-5 

-2.08775, 

-a 

-2.23031* 

■* 

-2.319421, 

-4 


I* 


4 
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table 11 .- f^T) - g(T) r p (T)— ►|a(T) T _ *(t) 009 (to + D) 


T 

' 



?< 0) fT 
x -7/a l 

) 

2&<’> 

;(o) 

(T) 

-(a) 

~1S/B 

(T) 

f (0) 

— -16/3 

(t) 

? (o) 

l -n /s 

(t) 

o 

T) 

f(o) 

■£-21/3 

(T) 

.02 

1.06625, 

0 

1.09362, 

0 

1.12160, 

6 

1.16020, 

0 

1.17965, 

0 

1.20881, 

(J 

1.24077, 

0 

1.27262, 

0 

1-306DB, 

a 

1 .33847, 

6 

J. 37273, 

a 

•04 

1.14626, 

0 

1.20401, 

0 

J. 26801, 

0 

1. 33640, 

0 

1.40686, 

0 

1.48111, 

0 

1.56983, 

0 

1.64274, 

0 

1.70004, 

0 

1.82199, 

0 

1.91882, 

0 

•06 

1.23401, 

0 

1,55750,' 

0 

1.44850, 

0 

1.66651, 

0 

1.89870, 

0 

1.63969, 

0 

1.99239, 

0 

£-16777, 

0 

2.35585, 

0 

2.65062, 

0 

2.74080, 

0 

.06 

1.54702, 

0 

1.60270, 

0 

1.67634, 

0 

1.87014, 

0 

2. 03627, 

0 

2.32740, 

0 

2.59640, 

0 

£ .89649, 

0 

5.25126, 

0 

3.60471, 

0 

4.02133, 

0 

.10 

1.49108, 

0 

1.71691, 

0 

1.97696, 

0 

2.27837, 

c 

2.52115, 

0 

6.01015, 

0 

®.*762T, 

0 

4.oniop, 

0 

4.60772, 

0 

6.50669, 

0 

6 .10918, 

0 

.12 

1.69092, 

0 

2**01466, 

0 

e.^oo*?. 

0 

2.96014, 

0 

3.40700, 

0 

4.06068, 

0 

4.83828, 

0 

3.70*8®, 

0 

S.B6869, 

0 

8*18440, 

0 

9.76184, 

0 

.14 

2,01567, 

0 

2.49113, 

0 

3.08200, 

0 

3.91LB9, 

0 

4.71736, 

0 

6.88624, 

0 

7.22049, 

0 

0,93506, 

0 

1.10613, 

l 

1.36731, 

1 

1.89151, 

1 

• 16 

2.29386, 

0 

2.69702, 

0 

3.668TB, 

0 

4.B207B, 

0 

6.63674, 

0 

T.3T080, 

0 

9,30606, 

0 

1.176*5, 

1 

1.48465, 

1 

1.87602, 

1 

2.38803, 

1 

.16 

Z.B661B, 

0 

S.4B7&4, 

0 

4,75260, 

0 

6.15143, 

n 

7.04492, 

0 

1.02618, 

1 

1.32631, 

1 

1.71172, 

1 

2.21078, 

i 

£.85632, 

1 

5.66784, 

1 

.165 

3.62501, 

0 

4.73706, 

0 

6 .19364, 

0 

9.0BB14, 

0 

1,06 BBS, 

1 

1.38440, 

1 

1.8100©, 

1 

2.3606'), 

1 

5.09441, 

1 

4.04692, 

1 

C. 26999, 

1 

.LT 

3,15201, 

0 

4.16093, 

0 

5.50817, 

0 

7.29103, 

0 

9.66300, 

0 

1.28066, 

1 

1.69729, 

1 

2.24940, 

1 

2.98127, 

1 

5.96118, 

1 

6.23776, 

1 

.175 

2.22438, 

0 

2.98579, 

0 

4.00240, 

0 

5.56666, 

0 

7.20117, 

0 

9.66907, 

0 

1.29667, 

1 

1.73771, 

1 

£.53071, 

l 

3.12604, 

l 

4,19256, 

1 

•16 

1.95587, 

0 

2.64806, 

0 

3.bBB8B, 

0 

4.B63C1, 

0 

6,59062, 

0 

8,029*6, 

0 

1.20086, 

1 

1.83901, 

i 

2.22013, 

l 

3.00721, 

1 

4.07277, 

1 

,166 

1.80649, 

0 

2. *8999, 

0 

5.37822, 

0 

4.51955, 

0 

6.51491, 

0 

6,63066, 

0 

1,17020, 

1 

1.61087, 

l 

2.10906, 

1 

5.00337, 

1 

4.09916, 

1 

•IS 

1.70780, 

0 

2.55524, 

0 

3.24646, 

0 

4.47263, 

0 

6.15B30, 

0 

8,47457, 

0 

1.16661, 

1 

1.50196, 

l 

2.20046, 

1 

3.02049, 

1 

4.14328, 

1 

•19& 

1.63709, 

0 

2.273 64, . 

0 

3.1&WJ7, 

0 

4.37275, 

0 

6.05468, 

0 

8.37431, 

0 

1.16890, 

1 

1.69628, 

l 

£.19955, 

l 

3.02663, 

1 

4.15914, 

1 

.20 

1.68294, 

0 

2.21194, 

0 

5.00092, 

0 

4.29791, 

0 

6.97006, 

0 

8,29063, 

0 

1.14793, 

1 

1.6R6B9, 

l 

2.18681, 

1 

5.00408, 

1 

4.11712, 

1 

•21 

1.60460, 

0 

2.12385, 

0 

2.9B4BE, 

0 

4.13271, 

0 

6.65304, 

0 

B.06766, 

0 

1.11719, 


1.65270, 

i 

2.06820, 

i 

2.81642, 

l 

3.76387, 

1 

.22 

1.46362, 

0 

2.06056, 

0 

2.90731, 

0 

4.07042, 

0 

6.64724, 

0 

7.76109, 

0 

1.06013, 

1 

1.S&99T, 

l 

1.82910, 

l 

£.32195, 

I 

2.83383, 

1 

.2* 

1.406TB, 

0 

2.01020, 

0 

2.63368, 

0 

5.93646, 

0 

6.37584, 

0 

7.18046, 

0 

6.51053, 

0 

1.16732, 

l 

1.34700, 

l 

1.39261, 

1 

1.35513, 

1 

.24 

1.37862, 

0 

1,96 TOT, 

0 

2.7S477, 

0 

3.70751, 

0 

4.29251, 

0 

6.82642, 

0 

7,47022, 

0 

7.T6T29, 

0 

6.81731, 

0 

-9.84841, 


-5.54253, 

0 

.26 

1 .36266, 

0 

1.92738, 

0 

2.08741, 

0 

5.55274, 

0 

4.47279, 

0 

5,12803, 

0 

4.84086, 

0 

2.26240, 

0 

-6,52010, 

0 

-2.28444, 

1 

-5.95406, 

1 

.26 

1.35*76, 

0 

1.68884, 

0 

2.5669* , 

0 

3.2*372, 

0 

3.T98ZH, 

0 

3.5M84, 

0 

1.51193, 

0 

-6.22684, 

0 

-£.05678, 

1 

-6.29048, 

1 

-1.17059, 

2 

.27 

1.31767, 

0 

1.85100, 

0 

2.46219, 

0 

2.96652, 

0 

2.94385, 

0 

1.60188, 

0 

-S.2MB4, 

0 

-1,49300, 

+1 

-4.04266, 

1 

-9.21996, 

l 

-1.92560, 

2 

.26 

1.50525, 

0 

1,81384, 

0 

2.31064, 

0 

&.56S45, 

0 

1.89654, 

0 

-1.01922, 

0 

-0.9196B, 

0 

-2.70640, 

l 

-6.62494, 

1 

-1.41217, 

2 

-2.85510, 

Z 

.29 

1.29167, 

0 

1.76740, 

0 

£ ■ l6^2£| 

Q 

£.03880, 

0 

6.470B3, 

-1 

-4.06081, 

0 

-1.67760, 

1 

-4.1T329, 

1 

-9.62262, 

1 

-2.00086, 

2 

-3.97822, 

2 

.50 

1.26167, 

0 

1.72193, 

0 

1.98746, 

0 

1.54132, 

0 

-8.?7091, 

-1 

-7*84462, 

0 

-2.59063, 

1 

-5.90469, 

1 

-1.30406, 

2 

-2.68362, 

2 

•*.24766, 

2 

.52 

1.26572, 

0 

1.62457, 

0 

1.66757, 

0 

2.02055, 

-1 

-4.46807, 

0 

-1.66319, 

1 

-4.39125, 

1 

-1.01213, 

2 

-2.14066, 

2 

-4.28186, 

2 


2 

.34 

1.25376, 

0 

1.BC325, 

0 

1.04595, 

0 

-1.40988, 

0 

-0.08118, 

0 

-2.77088, 

1 

-6.64998, 

1 

-1.61680, 

2 

-3.08170, 

Z 

-6.94329, 

2 

-1.04319, 

3 

.56 

1.24468, 

c 

1.55560, 

0 

4.16996, 

-1 

-3,54961, 

0 

-1.46620, 

1 

-4.09645, 

1 

-©•67291, 

1 

-2.04918, 

2 

-3.6T01T, 

2 

-7.04020, 

2 

-1.12320, 

3 

.52 

1.25773, 

0 

1.20838, 

0 

-3.30892, 

-1 

-6.90851, 

0 

-2.1156®, 

1 

-5.68241, 

1 

-1.25235, 

z 

-2.64098, 

2 

—4.67013, 

2 

-7.16660, 

2 

—8.86530, 

2 

.40 

1.23236, 

0 

1.02686, 

0 

-1.20024, 

0 

-8.PC069, 

0 

-2.8438T, 

1 

-7,15425, 

1 

-1.53925, 

2 

-2.68346, 

2 

-4.67837, 

2 

-6-45223, 

2 

-1.75978, 

£ 

.62 

1.22516, 

0 

3.18886, 

-1 

-2.19644, 

0 

-1.10750, 

1 

-3,63030, 

1 

-8.70278, 

1 

-1*78624, 

2 

*-2.48243, 

2 

-3.06640, 

2 

-1.20798, 

2 

1.12514, 

3 

.44 

1,22464, 

0 

6.62101, 

-1 

—3.32932, 

0 

-1.64729, 

1 

-4.44858, 

l 

-1,00866, 

*i 

-1.58253, 

2 

-2,88679, 

2 

-1.49106, 

1 

6.02 397, 

2 

2.99872, 

3 

•46 

1.22216, 

0 

3. 147B3, 

-1 

696 36, 

0 

-1.02628, 

1 

-6.25146, 

1 

-1.11209, 

2 

-1.80168, 

2 

-1,06787, 

2 

£.12700, 

2 

1.81487, 

3 

5.22666, 

3 

.48 

1.21989, 

0 

1.36241, 

-2 

-7.73300, 

0 

-2.52465, 

1 

-0.99988, 

1 

-1.18370, 

z 

-1.04856, 

2 

-8.32672, 

0 

7.20692, 

2 

2.62444, 

3 

7.53202, 

3 

.50 

1.21792, 

0 

-3.24763, 

-1 

-9.27525, 

0 

-2.7375T, 

1 

-6.63670, 

1 

-1.15815, 

2 

-9.54618, 

1 

-2.21096, 

2 

1.34626, 

5 

4.03726, 

3 

9. 29462, 

5 
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TABLE 12.- F ^( t ) *= b ( T ) T p ( t ) — ► Mt ) T p ( t ) cob (W - p , - 0 


^l/2.1 tT) 




4/3,l< T > 


S < o ) 
* 7 / 3 . l ' 

T ) 

Ao ) 

* 9 / 8. 1 

( T ) 

=( o ) 

I11 / 2 , 

l ( T ) 

m 




? ( o ) . 

* 17 / 8 , r 

t ) 

=( o )* 

- 19 / 8,1 

O ) 

1 . 0095 , 

0 

6 . 8455 , 

-1 

9 * 5677 , 

"1 

9 . 3683 , 

-l 

9 . 1248 , 

-1 

8 * 8973 ., 

-1 

8 . 6761 , 

-1 

1 . 0238 , 

0 

9.7212 

-1 

9 * 2506 , 

-1 

8 . 7648 , 

-1 

6 . 5224 , 

-1 

7 . 9024 , 

-1 

7 . 5036 , 

“I 

1 * 0450 , 

0 

9 . 6495 , 

-1 

8 * 9098 , 

-1 

8 . 2269 , 

-1 

7 * 6964 , 

-1 

7 . 014 *, 

-1 

6 . 4767 , 

-1 

1 . 0771 , 

0 

9 . 6649 , 

•1 

8 . 6491 , 

-1 

7 . 7680 , 

-1 

6 . 0642 . 

-1 

6 . 2538 , 

-1 

6 . 6879 , 

“I 

1 . 1276 , 

0 

9 . 7919 , 

-1 

8 . 6059 , 

-1 

7 * 3854 , 

-1 

5 . 4139 , 

-1 

6 . 670 *, 

-1 

4 , 8378 , 

-1 

1 . 2152 , 

0 

1 . 0182 , 

0 

8 . 5461 , 

-1 

7 . 1718 , 

-1 

5 . 0190 , 

-1 

6 - 0616 , 

-1 

4 . £ 896 , 

-1 

1 . 3865 , 

0 

1 . 1197 , 

0 

9 . 0607 , 

-1 

7 . 5166 , 

-1 

6 . 9131 , 

-1 

4 . 7796 , 

-1 

5 . 8652 , 

-1 

1 . 5517 , 

0 

l . essG , 

0 

8 , 7910 , 

-1 

7 . 71 * 5 , 

-1 

6 * 1546 , 

-1 

4 . 8606 , 

-1 

3 . 8488 , 

-1 

1 . 9644 , 

0 

1 . 6209 , 

0 

1 * 1749 , 

0 

9 * 1176 , , 

-1 

7 * 0595 , 

-1 

6 . 4058 , 

-1 

4 . 2318 , 

-1 

2 . 6698 , 

0 

1.9665 

0 

1 . 8082 , 

0 

1 * 1497 , 

0 

0 . 7952 , 

-i 

6 . 7252 , 

-1 

5 . 1436 , 

-1 

4 - 3205 , 

0 

3 . 2800 , 

0 

2 . 4600 , 

0 

1 . 6562 , 

0 

1 . 4607 , 

0 

1 * 0659 , 

0 

7 . 9 TBG , 

-1 

2 . 4996 , 

0 

1.6880 

0 

1 . 8912 , 

0 

1 , 0599 , 

0 

7 * 7430 * 
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